V.D, ZAKHAROV 


GRAVITATIONAL WAVES 
IN EINSTEIN'S THEORY ~~ 


This book is an up-to-date review of work on gravitational 
waves in general relativity. Its central theme is the exposi- 
tion of rigorous approaches to the problem, principally the 
definitions and criteria for distinguishing wavelike gravita- 
tional fields from other solutions of Einstein’s equations. 
The introductory chapter (Chapter 1) contains a review of 
approximation methods. The mathematical apparatus 
required for setting up the problem of rigorous (generally 
covariant) wave criteria —Cauchy’s problem for gravitation- 
al equations and the Petrov classification of gravitational 
fields — is given in Chapters 2 and 3. Chapters 4—8 describe 
the well-known generally covariant criteria for gravitational 
waves of Pirani, Bell, Lichnerowicz, Zel’manov et al. Relat- 
ed to these is Chapter 12, which is devoted to the chrono- 
metrically invariant analysis of gravitational-inertial waves. 
Chapter 9 deals with the theory of propagation of gravita- 
tional waves and their classification according to the nature 
of the wave front (plane and spherical waves). Chapter 10 
discusses the special case of spaces with plane gravitational 
waves, more specifically spaces that admit of an absolutely 
parallel vector field. Chapter | 1 reviews work on the asymp- 
totic behavior of wavelike gravitational fields generated by 
insular distributions of sources, while the last chapter gives 
a brief review of the problem of experimental detection of 
gravitational waves and of the basic results arrived at exper- 
imentally to date. 
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FOREWORD 


The modern theory of gravitation — Einstein's general theory of 
relativity — has received practically no experimental confirmation 
outside of the amazingly accurate predictions about the advance of 
the perihelion of Mercury and the deflection of light rays in the 
Sun's gravitational field. Moreover, in the course of its develop- 
ment the theory itself has generated a whole series of fundamental 
problems. The latter include the problem of gravitational energy, 
related problems of gravitational waves and quantization of gravita- 
tion, and a number of other unsolved ones, 

This situation is due first and foremost to the absence of the 
basis required for the construction of a theory representing reality, 
i.e., absence of experiments which can be duplicated and which admit 
of variation of the parameters. Experiments have been inhibited 
chiefly by the state of the experimental technology. Recently, how- 
ever, a concrete program of experiments on the detection of gravi- 
tational waves has emerged in the works of Weber, Schiff, Bra- 
ginskii and others. These investigations are certain to be followed 
by further experimental work, and the theory of gravitation is sure 
to acquire a more secure experimental base, 

Whenever a branch of science, although definitely relevant today, 
lacks a sufficiently solid experimental foundation, elements of 
formal theorizing begin to creep into its development. To some 
extent at least, the modern theory of gravitation is no exception: 
suffice it to mention the many variants of a "unified theory," by 
no means all of which will withstand the test of time. 

Of all the unsolved problems relating to gravitation, that of 
gravitational waves has attracted the most attention among physicists, 
theoretical and experimental. This is because it is closely related 
to other unsolved problems in this field (the energy problem, the 
problem of the construction of a quantum gravidynamics, and so on), 
and its solution on the theoretical and experimental planes would be 
a stimulus to the study of many other problems of gravitation, 

A vast corpus of work has been devoted to gravitational waves, 
including both large monographs and isolated articles dealing with 
theoretical approaches to the problem. Authors differ in their 
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approaches: one can find the approximate approach, based on analogy 
with Einstein's method; definitions based on geometric or physical 
considerations; or lastly theories based on the analogy with electro- 
magnetic theory. 

To construct a theory of gravitational waves it is necessary, of 
course, to start from some definition. Yet in Einstein's theory of 
gravitation, two fundamental difficulties arise at precisely this 
point. The first is that the geometrical object which determines 
the energy and momentum of the field is not tensorial in character, 
and the second is that the field of gravitation itself becomes iden- 
tical witn the space-time continuum, complicating the problem very 
considerably. One might mention that there exists at present no 
physical theory (discounting 'unified theories") in which such a 
situation prevails, 

As we saw, the number of approaches to gravitational waves is 
very large, and the first concern of the present author, V.D.Zakha- 
rov, has been to arrange them according to some basic principle. 
Having chosen the definition of the concept of gravitational waves as 
guiding principle, the author arrives at seven groups of theories. 
Of course, these groups are not very clearly demarcated, but at 
any rate the author has succeeded in putting some sort of order in 
the truly enormous material and making it more susceptible of 
review. 

The significance of the various theories considered in this book 
is worth discussing at some little length. It is evident that until 
experiments have yielded decisive results, no criterion can be 
given by which to estimate the adequacy of any given theory; many 
of the approaches described by the author will eventually seem to 
be of historical interest alone if one considers them on the whole 
as closed theories of gravitational radiation. Still, each one 
of them contains something of objective value —be it a construc- 
tive invention, a fine point of physics, a true physical idea, a 
mathematical device or something else. And although some of the 
approaches analyzed by the author are likely to fail the test of time, 
many elements which enter into their makeup will remain in the 
final theory, proving useful and necessary, . 

Thus at present writing gravitation research is still at a stage 
where one must speak of a theory in the making, rather than of a 
completed theory, Partly for this reason the book has the following 
structure: analysis, interpretation, arrangement of the vast 
theoretical material now available; the author's own results; and, 
to a lesser extent, comparison with new experimental ‘studies. 

This particular structure is easily explained in the light of the 
present situation in gravitation. In addition, of course, the book's 
contents depend on the personal views of its author. As the author 


points out in his introduction, the stress here is on works which 
give a rigorous invariant definition of the concept of gravitational 
waves, in the sense of necessary and sufficient conditions imposed 
upon the quantities characterizing the space-time manifold, in 
order that it describe a wavelike gravitational field. Other 
possible approaches to the problem are described in less detail. 
This appears to be a feature of all authors. 

And yet surely there is one very important rule which authors 
writing on any physical problem must obey: relate theory to experi- 
ment (and observation). Zakharov seeks to adhere to this rule, 
and a special chapter has been devoted to the experimental aspects 
of the problem. The trouble is that for the present there exist 
very few such comparisons with experiment (disregarding classical 
astronomical observations and first attempts at detection of gravi- 
tational waves by Weber, Braginskii and a few others). 

Five to seven years from now, many of the theories considered 
in this book will have been sifted out by the experimental sieve, to 
lie upon the shelves of history. But nothing of value in them will 
go to waste; rather, it will become incorporated organically in the 
fabric of the future theory of gravitation. Here Zakharov's mono- 
graph will play a useful part. On the other hand, for this very 
reason, as well as because of the large number of works involved, 
the author was physically unable to devote space to detailed 
analysis of any given author. Basically his plan is as follows: 
set forth clearly the principal premises of the theory; indicate the 
basic conclusions which follow from these premises; sketch in the 
contours of the theory on the basis of the foregoing. Occasionally 
the theory is linked to other theories, equivalences between 
theories, if any, are established. A more detailed exposition would 
have meant enlarging the book by a factor of four or more, a size 
hardly conducive to legibility and possibly defeating the author's 
principal objective — to give the fullest possible review of the 
present state of the problem of gravitational waves. As to the 
author's own studies, they have been described objectively, and 
inserted in their proper place, on equal terms with other 
approaches, and connections with the work of other authors 
(mainly Zel'manov's theory) have been indicated, 

A theory of gravitational waves within the framework of 
Einstein's theory of gravitation — speaking now of the fundamental 
difficulties encountered in constructing the theory — involves 
answering at least three questions (which entrain, incidentally, 
other questions as well): a) if one attempts to introduce integral 
invariants based on the ordinary theorem of Ostrogradskii 
[Gauss], the gravitational field energy is expressed not in terms 
of tensors (which would have been expected) but in terms of 
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geometric objects which can be made to vanish in a chosen coordi- 
nate system; b) when one defines the field energy in terms of 
tensor quantities and equates the field and the space-time continu- 
um according to Einstein, it becomes necessary to disentangle the 
concept of "energy of space and time’; if one discards the concept 
of energy and momentum of the gravitational field (a standpoint 
also to be found), any analogy with other fields must be deleted in 
advance; c) since the concept of a privileged coordinate system 
does not occur in the general theory of relativity, difficult problems 
of coordinate interpretation arise in concrete cases. These and 
other problems Zakharov presents acutely enough to his readers, 
affording them a purely constructive grasp of the principal diffi- 
culties in present-day theory of gravitation. Certain authors of 
popular publications either skirt the problem or simply fail to 
understand it, painting an optimistic picture in which the general 
theory of relativity is made to appear complete and capable of 
answering all questions put before it. 

Today one of the oldest branches of physics — the science of 
gravity — is undergoing a period of reevaluation based on contem- 
porary experimental research, and one may assume that the theory 
of gravitation will be presented very soon with more substantial 
data. In particular the latter would be of use in studying the very 
interesting phenomenon of gravitational waves — if these exist. It 
would then be possible to reevaluate the hypotheses discussed in 
Zakharov's monograph, picking out such elements as have been 
confirmed experimentally. 

Is it possible, from the material given by the author and, in 
general, from the present situation, to arrive at a preference for 
a particular theory of gravitational radiation? The author of this 
foreword believes, from data known to him, that this cannot be 
done at present. Experimental proof is lacking; a more funda- 
mental base is needed, and the theories under consideration are 
heuristic in nature. Strictly speaking — so the author of the fore- 
word believes — none of the structures considered in the book will 
merit the name of theory, in the strict sense of the word, until it 
has been verified experimentally and has led to experimentally 
detectable consequences, that is until it has ''begun to function" 
within this plan, ; 

Furthermore, new theoretical constructions introduced as the 
monograph was being sent to press (e. g., based on the principle of 
modelling of gravitational fields) have not beeh included. This 
unavoidable occurrence does not detract from the overall value of 
the book and the overall good marks which it deserves. 

Mention should be made of the truly tremendous bibliographical 
task accomplished by the author, which required considerable 
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knowledge of modern physics and of the mathematical apparatus. 
The bibliography (about 450 items) is of exceptional value thanks to 
the thematically well thought out selection and it is a great time- 
saver for the reader who would otherwise drown in the sea 

of information. With the possible exception of the last two years, 
it reflects practically all major work on gravitational waves. 

In evaluating Zakharov's monograph it should also be recalled 
that it is of immediate relevance today. Of this the author is well 
aware, and he unfolds a sweeping and varied view of the researches 
of many different theoreticians (Soviet scientists being well repre- 
sented) while deriving his book from the rigid frameworks of text- 
books now largely outdated as regards presentation of problems 
and especially as regards physical data. Today, with the flow of 
information threatening to engulf the reader, the need for such 
books, for review works on particular problems, is especially 
acute. Although not entirely free from flaws (Somewhat Schematic 
exposition, subjectivity in the evaluation of certain hypotheses, and 
so on), Zakharov's monograph has the unqueStionable merit of being 
the first attempt in the literature to sketch a sweeping picture of 
the various theoretical constructions related to the problem of 
gravitational waves. There is at present no work of this type deal- 
ing with gravitational waves, despite the obvious need for one. 

For this reason there is no doubt whatsoever that Zakharov's 
monograph will be swift to find readers. 


A. Z. Petrov 
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AUTHOR'S PREFACE 


In the present state of the problem of gravitational waves it is 
impossible to present a mathematical exposition of the theory (in 
practice, as yet unconstructed) of gravitational waves, or of 
laboratory methods for their detection. The material at the 
author's disposal consists of disconnected ideas and approaches, 
rigorous or approximate in varying degrees. Any review of these 
must of necessity be fragmentary to some extent, the more So as 
many have not been developed fully enough, despite the extensive 
literature in the journals. 

This being the situation, with subjective under- or over-evalua- 
tion of a given result not being ruled out, the author has decided to 
concentrate on concepts admitting of mathematically rigorous 
exposition, This has led to a definite selection of themes anda 
definiteemphasis. More specifically, the central theme of the book 
is the problem of the generally covariant formulation of a criterion 
for gravitational waves, resting upon mathematically rigorous 
methods for solving Cauchy's problem for Einstein's equations of 
gravitation. Cauchy's problem and the principal results of the 
theory of characteristics for Einstein's equations are given in 
Chapter 2. 

In Chapter 3 the statement of the problem of gravitational waves 
is formulated and two main methods of describing them — the one 
algebraic and the other differential — are presented. Necessary 
information from the basic mathematical apparatus used in the 
theory of gravitational waves — Petrov's algebraic classification of 
gravitational fields —is also given here. Both these chapters are 
introductory with reference to Chapters 4—8, which are devoted to 
algebraic methods of invariant description of gravitational waves. 
Chapter 12 and parts of Chapters 7 and 8 are devoted to differential 
methods. 

The second most important theme of the book consists of prob- 
lems of propagation of gravitational wave fronts, as well as of an 
investigation into the asymptotic properties of fields of gravitational 
emission whose sources are isolated material systems (Chapters 9 
and 11). 
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Finally, exact and approximate solutions of Einstein's 
equations describing gravitational waves are analyzed throughout 
the book (except for Chapters 2 and 3), solutions describing plane 
gravitational waves being specially examined in Chapter 10. 

These three themes constitute the main outline of the book and 
have determined its overall aim: to serve as a ''working review," 
singling out the most complete and, in the author's estimate, most 
promising approaches to the study of gravitational waves, In this 
context mention should be made of two other themes discussed in 
the book: approximate methods of investigation of wavelike gravita- 
tional fields (their review is the subject of the introductory chapter), 
and analysis of methods and prospects in the experimental detection 
of gravitational waves; the latter plays the part of physical conclu- 
sion to the book as a whole (final chapter — Ch, 13). 

This particular orientation of the book also determined the treat- 
ment of the problem of the energy transmitted by gravitational 
waves. Since at this point the problem of the theory of gravitational 
waves intersects another no less complex problem — that of the 
energy of the gravitational field —it is also the focus about which 
are grouped the fundamental difficulties encountered in the modern 
theory of gravitation. 

The author's personal opinions have obviously played some role 
in the selection or in the elucidation of material. The reader may 
discern a tendency on the author's part to reduce a rapidly evolving 
branch of science to a dead mathematical skeleton, the result of 
the vivisection of a living, evolving organism. Here the author may 
be exposing himself to sneering remarks, in the spirit of Goethe's 
Mephistopheles: 


Wer will was Lebendigs erkennen und beschreiben, 
Sucht erst den Geist herauszutreiben, 

Dann hat er die Teile in seiner Hand, 

Fehlt leider! nur das geistige Band. 


However, the author did not set himself the lofty aim of ade- 
quately depicting ''the splendid colors of living nature,'' mindful as 
he is of the fact that the ''tree of knowledge" is not the ''tree of life." 
Without a ''mathematical skeleton" there can be no theory; this is 
why the question of the interaction of gravitational waves and the 
atomic structure of matter, in particular, have been regorously 
excluded from the book. A satisfactory solution of this problem 
will only be achieved, it seems, providing there is substantial 
generalization or extension of Einstein's theory to the phenomena of 
the microuniverse. The creation of such a ''quantum theory of 
gravitation'' will be the concern of the future; in our book the 
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problem of gravitational waves is considered exclusively on the 
basis of the classical Einsteinian theory of gravitation. 

It is the author's hope that the publication of this book will to 
some extent make up for a severe deficiency in the systematic 
reviews on gravitational waves, a deficiency which is acutely felt 
today inspite of the fact that the reader has at his disposal many 
outstanding monographs dealing more or less with the same 
material as the present book. Thus in Weber's "General Relativity 
and Gravitational Waves" the problem of gravitational waves is 
examined only insofar as needed by the author for possible labora- 
tory methods of detection; moreover, the book's publication (1961) 
dates back to a period when most modern lines of investigation had 
not yet been defined. Braginskii's recently published thorough 
monograph, ''Fizicheskie eksperimenty s probnymi telami" 
(Physical Experiments with Test Bodies) deals with the problem at 
an up-to-date level but is limited to the experimental approach to 
gravitational fields. Other reviews (lectures by Pirani and Sachs, 
Petrov's review and so on) are either of necessity very brief, or 
devoted exclusively to the research of the authors themselves. 

The book assumes previous acquaintance with the foundations of 
the general theory of relativity and is intended mainly for students 
at higher levels and candidates for a degree, as well as for all 
specialists in gravitation to whom it may be of use as a reference 
manual. For this reason the author has sought to provide as com- 
plete a bibliography as possible, without claiming, of course, to 
list all works on gravitational waves. Well aware of the risks he 
assumed when he took upon himself to write a broad monographic 
review of so vigorously evolving a problem as that of gravitational 
waves, the author apologizes in advance to those colleagues whose 
works might have significantly enriched the book but escaped the 
author's attention. 

The author is grateful to reviewers A. Z. Petrov, 

I, D. Novikov and N.V. Mitskevich and editor V,N. Zakharov for 
reading the manuscript, as well as to Ya, B. Zel'dovich, 

K.P. Stanyukovich, K.S. Thorne, V.B. Braginskiiand A. L. Zel'manov 
all of whom contributed in some form to the book's publication. 
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INTRODUCTION 


The problem of describing gravitational waves theoretically, 
which is intimately related to problems arising in their experi- 
mental study, has become one of the most pressing and interesting 
problems not just of gravitation but of modern physics in general. 
Emerging almost simultaneously with Einstein's theory of gravita- 
tion (the first analysis of the problem was carried out by Einstein 
himself in 1916 — 1918), it has so far defied fully satisfactory 
solution. 

In the last ten to fifteen years (roughly since 1957), interest in 
this problem has increased considerably owing to the creation of a 
powerful new mathematical apparatus. This is the Petrov classi- 
fication of gravitational fields, which has initiated many new 
theoretical approaches. On the other hand, recent advances in the 
experimental field and in particular Weber's experiments hold out 
prospects of laboratory detection of gravitational waves. 

Heightened interest in the question of gravitational waves has led 
to numerous publications in the journals; these now number in the 
hundreds. Such works may be divided into several groups reflect- 
ing different lines of research. 

The first group of works seeks to give a rigorous definition of the 
the concept of gravitational waves, i.e., to formulate, in generally 
covariant form, the necessary and sufficient conditions which the 
Space-time metric must satisfy in order to describe a wavelike 
gravitational field. To this group belong the works of Pirani /69, 
262, 263/, Lichnerowicz /62, 87—90, 264/, Bel /56, 68, 75—80/, 
Debever /66, 81, 109, 265—268/, Trautman /55, 269, 270/, 
Petrov /271/, Ehlers and Sachs /185, 272 —274/, Hely /111—117/, 
Roy and Radhakrishna /96/, Zakharov /94, 100, 101, 275, 276/, 
Staruszkiewicz /277/, Parizet /278/, Zund and Levine /118 —120/, 
Misra and Singh /127, 128/, Maldybaeva /121, 123, 279/, Sokolik 
and Konopleva /92, 93, 124/, Aichelburg /280/, Lukaéevié /281/, 
Coburn /282/, Yadav /283/ and Nikolaenko /126/. 

The above workers treat the definition of gravitational waves in 
its purely geometric aspect, on the basis of the albegraic properties 
of gravitational fields as defined inPetrov's classification /57, 64, 65/. 
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In the second group of works the authors start from a chosen 
definition of the energy of the gravitational field and proceed to 
define gravitational waves. These works, in contrast with the 
geometrical works of the first group, are physical in approach, and 
include the researches of Infeld /8, 284—287/, Synge /83/, Peres 
and Rosen /288 —290/, Arnowitt, Deser and Misner /291 —293/, 
Geissler, Treder and Papapetrou /294/, Araki /295/, Brill /296/, 
Méller /297—299/, Gutman /300 —302/, Shirokov and Bud'ko 
/303 —305/, Petrov /306/, Wu Tthan Khiet /307/, Denisov /308/, 
Signore /309—311/, Isaakson /38/, Rodichev and Dozmorov /312 — 
314/, Zakharov /315 — 316/. It should be mentioned that certain 
authors (e.g., Synge, Wu T'han Khiet, Mgller (1961), Rodichev 
and Dozmoroyv, Denison, Gutman, Isaacson) employ a generally 
covariant (or tetrad) definition of the gravitational field energy. 

In this respect (covariance of the criterion) their approach may 
equally well be classed with the first. A traditional 'pseudoten- 
sorial'' approach was used by Peres, Rosen, Mgller (1958) et al. 
Lastly, Araki, Brill, Geissler, Treder and Papapetrou determine the 
energy of the gravitational field ina specially chosen coordinate system. 

The third group comprises works which study either waves ofa 
specific form (plane, spherical), or the gravitational emission of 
isolated systems of sources. These include: plane waves — Rosen 
/155/, Boardman and Bergmann /156/, Bondi, Pirani and Robinson 
/143/, Kundt and Ehlers /137, 145—147/, Weber and Zipoy /95, 
317/, Kerr and Goldberg /318/, Avez /157/, Newman /150/, 
Penrose /144/, Chevreton /149/, Johari /154/; spherical waves — 
Robinson and Trautman /129, 131/, Cahen and Leroy /319, 320/, 
Foster and Newman /138/, Marder /321/; emission of isolated 
systems — Bondi /20, 171, 322 —325/, Stachel /188, 203/, Janis, 
Newman, Torrence and Couch /182, 326/, Hawking /202/, Biéak 
/327, 328/, Van der Burg /329/, Isaacson, Winicour and Derry 
/330/, Le Denmat /331/, Madore /45, 46/, Persides /197/, 
Halliday and Janis /332/, as well as works of Sachs /110/, Newman 
and Tamburino /333/, Szekeres /134/, Unti and Torrence /334/, 
and Collinson and French /335 —336/, 

The fourth group is composed of works on exact or approximate 
solutions of Einstein's equations describing gravitational waves in 
the sense of a certain criterion. Exact wave solutions known today 
and analyses of these may be found in the works of Einstein and 
Rosen /187, 337, 338/, Takeno /102, 153, 163, 168, 339 —347/, 
Petrov /57, 97, 348/, Weber and Wheeler /193/, Marder /190 — 
192, 349/, Lichnerowicz /89/, Geissler and Treder /350 —352/, 
Kompaneets /189, 353, 354/, Peres /108, 160/, Robinson and 
Trautman /129, 130/, Pandya and Vaidya /164, 355—357/, Sciama 
/358/, Bonnor /359, 360/, Friedlander /361—363/, Nordtvedt and 
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Pagels /104/, Krishna Rao and Pandey /194—196, 364, 365/, 
Harrison /366/, Leroy /367/, Wyman and Trollope /72, 73/, 
Zakharov /91, 101, 103, 105, 107, 170, 368/, Johari /154, 369/, 
Misra /370, 371/, Bartrum /135/, Foster and Newman /138/, Lal 
and Prasad /372/, Dangvu /159, 373/, Hoffman /374/, Dozmorov 
/375—377/, Szekeres /378/ and Aichelburg /379/. 

Approximate wave solutions are investigated in the works of 
Rosen and Shamir /19/, Bonnor /11, 12, 17, 380—382/, Pirani 
/383/, Peres /384/, Lias /385/, Mehra, Vaidya and Kushwaha 
/386/, and Murenbeeldand Trollope /387/. A general method for 
constructing approximate wave solutions is given in the works of 
Choquet-Bruhat /388, 389/. 

In the fifth group of works gravitational waves are treated by 
approximation methods: either using linearized equations of gravi- 
tation (Einstein /1, 2/, Eddington /23/, Matte /82/, Dirac /390/, 
Vavilov /391/, Gertsenshtein and Pustovoit /392, 393/, Bonnor /6/, 
Carmeli /394/, Cooperstock /24, 243, 395/, Rotenberg /21, 25, 
29/, Campbell /22/), or by representation of the equations in an 
approximate form of specified order of smallness (Bonnor /11, 12, 
17, 396/, Fock /9, 397/, Infeld /8, 398/, Papapetrou /399—401/, 
Tonnelat /402, 403/, Isaacson and Winicour /37, 38, 44/, Treder 
/404/, Cooperstock /405/, Unt /406, 407/, Zerilli /408/, Vish- 
veshwara /179/, Couch, Kinnersley and Torrence /180, 181/), or, 
finally, by deriving the wave equation from modified (''Maxwellized") 
equations of gravitation (Rumer /409/, Kroki /410/, Mavrides /411, 
412/, Synge /413/, Berger/414/ ). 

Investigations in the sixth category deal with the gravitational 
emission of elementary particles. First and foremost is a mono- 
graph by Stanyukovich /415/ which reviews previous works in this 
field by the author (/416—419/), followed by works by De Witt 
/420/, Kundt and Thompson /186/, Halpern, Laurent and 
Desbrandes /254, 421/. Inview of the fact that Einstein's theory 
of gravitation is inapplicable to the description of microscopic 
systems, Stanyukovich replaces Einstein's equations with others 
involving a variable ''constant'' of gravitation. As mentioned 
earlier, works belonging to this trend are not touched upon in 
this book. 

The seventh group consists of works devoted to problems of ex- 
perimental investigation of gravitational waves, primarily detection. 
Workers in this field include Weber /31, 95, 212, 227—229, 422 — 
429/, Braginskii, Rudenko, Rukman /213, 237, 430 —434/, 
Gertsenshtein and Pustovoit /435, 436/, Kopvillem and Nagibarov 
/247—251, 255, 256/, Bashkov /437, 438/, Mironovskii /241, 439/, 
Petrov /306/, Forward and Berman /215/, Heintzman /242/, 
Winterberg /244/, Zipoy and Bertotti /245/, Dyson /258, 259/, 
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Slabkii /440/, Vodyanitskii and Dimanshtein /465/, Lavrent'ev 
/252, 253/, Dautcourt /441/, Melosh /225/, Wick /442/, Papini 
/443/, Boccaletti, de Sabbata, Gualdi and Fortini /444/. 

A number of works in this category deal with estimates of the 
power of gravitational radiation from cosmic sources, and 
prospects for laboratory study of the latter; these include work by 
Wheeler /36/, Fowler /445/, Zel'dovich and Novikov /34/, Thorne 
/30, 218—222, 446/, Cooperstock /32, 33/, Carmeli /35/, Weber 
/31/, Boccaletti, de Sabbata, Gualdi and Fortini /239, 240/, 
Weinberg /447/, Shklovskii /216/, Greenstein /230/, Kafka /232, 
233, 448/, Sciama, Field and Rees /231, 234, 235, 449/, Kaufman 
/450/, Peters /451/, Alladin and Sastry /452/, Chandrasekhar 
/453 —455/, Ezawa /456/, and Chau and Henriksen /457, 458/. 

The purely astrophysical aspect of the problem and in particular 
the role of gravitational emission in the energy balance of cosmic 
objects has been allotted limited space. A detailed exposition of 
this specialized astrophysical question would have stood out against 
the general outline of the book, which is devoted primarily to 
theoretical aspects — questions of principle — of the problem of 
gravitational waves, A rapid summary of experimental results and 
problematics in this category is given in Chapter 13. 

Works discussing the velocity of propagation of gravitational 
waves — by Finzi /59/, Arifov /459/, Krzywoblocki /460/, 
Mitskevich /461/ — constitute a group apart. Various aspects of 
the problem of gravitational waves are the subject of reviews by 
Pirani /70, 71, 262/, Sachs /148/, Trautman /131/, Petrov /462/, 
Bondi /463, 464/, et al. (e.g., the review section in /165/). 

In conclusion we note that the bibliography consists chiefly of 
works published prior to 1971. 
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Chapter 1 


APPROXIMATION METHODS FOR THE 
INVESTIGATION OF GRAVITATIONAL WAVES 


1. Linear approximation 


The fact that radiative and wave processes are a natural element 
of the modern theory of gravitation — that is, of Einstein's general 
theory of relativity — becomes apparent when one considers the 
case of weak gravitational fields. Here the overall picture becomes 
so strikingly simplified that many complex concepts and relations 
of this theory become analogous to known ones in classical field 
theory. Although the main subject of this treatise is the problem of 
the invariant formulation of the concept of gravitational waves, it 
would not be complete enough or motivated well enough unless we 
were to carry out a preliminary analysis of the methods of approx- 
imate description of wavelike gravitational fields. 

Einstein's equations 


Rag = — 4 (Top — 5 Tae) (1.1) 


were first treated in the weak-field approximation by their author 
/1, 2/. Here R,, is the Ricci tensor, Tf,, is the energy-momen- 
tum tensor of "matter," i.e., of matter and of all fields other than 
the gravitational one, T = g*T,, anddis Einstein's gravitational 
constant. The summation convention is employed throughout the 
book; the Greek indices run from 0 to 3, and small Latin indices 
from 1 to 3. 

If the metric tensor g,,g is only slightly different from the 


: 2 . _ (00) 
Minkowski metric gas, 


(00) 
Bap = Bap + hap, . (1.2) 
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then the quantities h,, are Small compared with unity. We assume 
that their partial derivatives are of the same order of smallness:* 


apy ~ Rapp ~ hag. (1 <3) 


Theninthe linear approximationin haggthe Riemann-Christoffel 
tensor is given by 


1 
Rapys =~ o (hay,p8 + hes,va = has+B —— Ney, a8). (1 4) 


We find the expression for the Ricci tensor from the above and 
write equation (1.1) in the linear approximation: 


4 (00) 4 (00) 
9 (g°®hay,ps i Ray — hove a hy as) =—h (Tax m7 ‘A Bax) . (1 . 5) 


Here we exploited the obvious fact that in the linear approximation 


(00) 00). 


a a a Osa (00) 
ge gh — n,n = he ge, h= gM hag. 


( 
Multiplying (1.5) by a, we arrive at the scalar equation 
Dh + hi = — AT, (1.6) 
where** 
(09) 
Cl = —‘g*d,d, = A — 0,0, (1,7) 


is the d'Alembertian of the special theory of relativity and A the 
Laplacian: 


* Strictly speaking, derivatives of Ay, of different orders have different dimensionalities; their 
orders of smallness should therefore be evaluated with reference to the characteristic linear 
dimensions (such as the radius of curvature) of the space-time background. But as, in this case, 
the latter is flat (infinite radius of curvature), the orders of smallness of the first and second 
derivatives can be regarded as the same in the approximations under consideration. Estimates 
for a nonflat background metric may be found, for instance, in the works of Isaacson, discussed 
later in the present chapter, 

In most sections of the book, partial derivatives with respect to the coordinates are denoted by 


a i , ae 
indices following a comma (for example, Qap 1 = 5yy Sap? while covariant derivatives are 
7 x 


7) 
denoted by indices following a semicolon (for example, Quy. = a Qu» Ty ius isa) 
H x 


however, wherever it is necessary to bring out the operator nature of these indices, the more 
explicit notation Qyg y= 9Qgg and Qjg., = V,Qa, is used. Other instances in which a 
special notation is used for differentiation are indicated in the text. 
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3 
A= >)0,0;. (1.8) 
i 


(00) 
Multiplying equation (1.6) by g,, and introducing the expression 


(00) 
obtained for ATg,, in (1.5), we arrive at a system of equations which 
is conveniently written as follows: 


(0) 
C ‘pap + ‘pi ay + "PE ay a Zap wv = 2hT apy (1 9) 
where we introduced the quantities 


(00) 


4 
Pap = hap — = has. (1.10) 


The system (1.9) can be simplified further if we recall that 
Hilbert's conditions can always be satisfied in the weak-field 
approximation /3, 4/: 


Pa,p = 0. (1.11) 
The field equations then assume the standard form 
Crpas = 207 ap, Yap = 0. (1.12) 


Thus in the linear approximation Einstein's equations are wave 
equations for the potentials ,,, their right-hand sides describing 
the sources of the gravitational field. Consequently, in the linear 
approximation the equations of gravitation describe the propagation 
of gravitational waves with the fundamental velocity* ec (z° = ct); 
conditions (1.11) are the analogues of the condition of gauge invari- 
ance in classical field theory. 

Let us choose the origin of the Cartesian system of coordinates 
within the volume V occupied by the sources. Let # be the radius 
vector of an arbitrary point P lying outside V and ~, the radius 
vector of an arbitrary point 0 within V. The general solution of the 
system of equations (1.12) for zero initial values (tug =0 and ‘thyp9 = 0 
for x°=0) is given by** 

T tt—lr 
es Te (1.13) 
v 
* The system of units with ¢ = 1 will be used throughout; the only exception will be the beginning 
of Section 2 of the present chapter. ; 
** The mathematical theory of the inhomogeneous d'Alembert equation is given (for example) in 
Sobolev’s book /5/. 
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where r=a—£&. Thus the solution of Einstein's equations in the 
linear approximation is given by the retarded potentials (1.13), 

By expanding the integrand in (1.13) in powers of |E|/|z], one can 
study the multipole emission of the material system to all orders. 
The corresponding components of the radiation are characterized 
by the multipole moment tensors. The rank s of the tensor of the 
25 -pole moment characterizing the mass distribution of the sources 
is determined by the number sg of the corresponding term in the 
multipole expansion. Thus, introducing the following three-dimen- 
sional tensors (Bonnor /6/ ): 


M=\ Ty.aV, Ay =\ To, aV, S.3=\ Tyav, 
Vv Vv Vv 


Mum =\ Too babs.» Em dV, 

y (1,14) 
Aapna..m = \ Tosbibs- «Em OV, 

v 
Sisntim =) Tikiks En dV, 

Vv 


and defining the corresponding traceless multipole moment tensors 
in terms of these, 


My = 3M — 81M pps 
Mim = OM kim == 81M ppm —— 8imM kpp mai SmxM pip, 
Ait = 3A ix a Six Apip» 


(1.15) 
a 4 
Sai = Sige + 94; (Skpip — 2S pp) + > Sri (Sppii—38 spp) an 
1 
+ 5 Six (Spot — 3Sipip)s 
one obtains the multipole expansion of the components of the 
retarded potentials »°: 
4S ij Ax, (3 if 
Y= — sta] — Bye Sat = ee 
4 ‘ 2 2 . 3A, 
Poi = — x27 4 (4 ilk i ae Mas) se 4 (4 ilkl + ie a ar +s, (1.16) 
oe 4M 2oR zy BM ee 3M a 
tom — Te ane (tn + Tt + Be) 
2r Krom 6M kim, 15M itm 15M itm 


aan [2 ie (11 sim ec 


[z| * [ep + qeh Jhon, 
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where the dots above a symbol designate differentiation with respect 
to the retarded time t—r. The first term of the expansion for y,, 
is obviously the Newtonian potential, while the succeeding terms 
correspond to the quadrupole and octupole moments of the mass 
distribution. The absence of the dipole term shows that, in particu- 
lar, a spherically symmetric system of sources cannot emit 
gravitational waves, Thus the solution (1.16) of the wave equation 
(1,12) describes gravitational waves in the linear approximation, 

the first ''radiative' term in the expansion of the potentials tog 
being the quadrupole term. 


2. Higher-order approximations 


To obtain nonlinear approximations of the equations of gravita- 
tion, we can use the classical Einstein-Infeld-Hoffman method of 
approximations /7/, which was first used to study the equations of 
motion and is based on expansion of the potentials of the gravita- 
tional field in the small parameter 1/ c¢ (e being the fundamental 
velocity)*. This method has been applied to the problem of gravi- 
tational radiation by Infeld (see /8/), Fock /9/ and Bonnor /6, 11/. 
Expansion in the parameter 1/c makes it possible to obtain the 
equations of Newton's theory in the zeroeth order, which in turn 
makes it possible to substantially simplify the treatment of higher- 
order approximations. Thus by choosing a suitable coordinate 
system one can ensure that in the expansion of the components hy 


Toys = Dy 07" heyy (24) 
n=0 (n) 


the first nonzero terms will be of second order for hy, and hj, and 
of third order for hg. 
Correspondingly, the first nonzero term of the expansion of %,, 


will be J; for tp; it will be };, and for »,,;, the terms py, of the 
(2) (3) (A) 
fourth order of smallness. The field equations in these approxima- 


tions become 


APoo = 2T oo, Aoi = 27 oi; Api; = 2nT 5; + Nay, (1 ah 7) 
(2) (0) (3) (2) (4) (4) (4) 


* To the formal expansion in the parameter 1/¢ corresponds an actual expansion in a dimension- 
less parameter of the type U/c?, where U is the Newtonian gravitational potential, or in v*/c?, 
where v is the velocity of one of the bodies forming the system of sources. 
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where N,; are the nonlinear terms, first appearing as corrections 
(4) 
to terms of the fourth order. 

As Gupta shows /10/ (see also Bonnor /11/), analogous 
equations can be obtained in the n-th order approximation by 
another method based on expansion of the metric ga, and the energy- 
momentum tensor [7,, in the parameter A (gravitational constant): 


(00) X n ~ 
Bap (x, 4) = Bap + Dy i hes Tap (22, 4) = DV N"P ag (2°), 


n=1 n=o () 


where hap and Jos are independent of 4. By inserting these ex- 
nr n 

pansions in the field equations one can express the latter as equa- 

tions (with zero on the right side) for analytic functions expandedin 

series in 4. Equating the coefficients of A" to zero, we obtain the 

equations of gravitation in the n-th order approximation, Making 

use once more of the definition (1.10) for 3, we can write (for 

n > 2) 


Atha = 2 [Tas — Nag], (1.18) 
() (@—) (nN) 


where N,, denotes nonlinear terms dependent only on combinations 
() 


of ~., of order lower than n, and independent of tg and A. 
(n) 


Equations (1.18) are taken, by definition, to be Einstein's equations 
of gravitation in the n-th order approximation (Havas and Goldberg 
/13, 14/). 

In approximations of arbitrarily high order, equations of the 
type of (1.17) —(1.18) may no longer be interpreted — generally 
speaking — as wave equations. However, in the case of isolated 
distribution of matter it is natural to assume that the classical 
concept of mass-energy of the system of sources used in the linear 
approximation is also applicable for higher-order approximations. 
Using this assumption and starting from the solution of the 
equations of gravitation of arbitrarily high order, one can calculate 
the variation of the energy of the system over a characteristic time 
interval (e. g., over a complete period of 2° -pole oscillations) with 
the aid of the energy-momentum pseudotensor which determines the 
energy transfer by gravitational waves in the linear approximation 
(Eddington /15/, Landau and Lifshits /16/). 

Another method which can be used to calculate the energy loss 
Am=m,—m, consists of comparison with the Schwarzschild field, 
provided the stationary states of the system at the instants ¢, and 4, 
corresponding to the mass values m, and m,, can be described by 
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Schwarzschild's solution. Obviously, this method is applicable 
only to those source distributions (and to such systems of 
reference) which go over to Schwarzschild's metric at sufficiently 
large distances from the system of sources. Then the change Am 
in the mass of the system over the oscillation period At = #4, — t, 
can be interpreted as a consequence of the transfer of energy by the 
gravitational waves. 

With this aim in view, the so-called method of double 
parameter approximations, which extends the foregoing 
method of expansion of the metric in powers of the gravitational 
constant, was put forward by Bonnor and Rotenberg /17, 18/ in 
application to sources of the insular type(in vacuum 7,, = 0). If 
the quantity m characterizing the total mass of the system is used 
as a parameter to replace the gravitational constant 4, and the 
characteristic parameter a of the system, with the dimension of 
length, is chosen as the other parameter which arises naturally on 
definition of the multipole oscillations of the system, then the 
solution of the equations of gravitation which can be expanded in 
converging Taylor series in m anda in the neighborhood of m=0 
and a=0 may be written as 


co oo ( 3) 
£e= >| bs mPas Zaps (1.19) 


p=0 8=0 


where the coefficients ek are independent of m and a. Ordinary 
expansion in powers of the gravitational constant is equivalent to 
expansion in the single parameter m(s= 0). 

Inserting this expansion in the equation for the field in vacuum 


Rug = 0 


and setting the expansion coefficients before m?a* equal to zero, 
we obtain a system of ten second-order differential equations known 
as the ''field equations in the ps -approximation": 


(ps) (ps) (ar) 
Dap (8uv) = Vag (8yy) (qxp—i, rss). (1,20) 


(ps) 
The left-hand sides of the above are linear in Luv (and their 
(qr) 
derivatives) while the right-hand sides are nonlinear in ae (and 
their derivatives), the latter being.already known from the approx- 
imation of preceding order. 
The 00 approximation obviously corresponds to the flat space- 


: (00) : : : 
time metric: g.5 = g.s.. All approximations of order 1s are linear 
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(1s) . 
and homogeneous in g,g and their derivatives, 


(18) 
ag = 0, 


and are therefore equivalent to the Einstein-Infeld-Hoffman approx- 
imations discussed earlier. Approximations of ps -orders with 
p> 2 are nonlinear; 2s-approximations correspond to the second 
order of smallness (s= 0, 1, 2, ...), 3s-approximations to the 
third order (s=0, 1, 2, ...), and so forth, 

The solution of the equations of gravitation in the linear approx- 
imation in the form of a multipole expansion is given by formulas 
(1.16). In our notation (1.19), it can be written as 


(00) = O) 
Bap = Baa -+ >) ma° Bap, 
=0 


(11) (10) 
where dipole terms are absent ( gap = 0), monopole terms gag (i. e., 


(00) (10) 
the static part, ge, + mgag ) correspond to the linear approximation 


. . ; (2) (18) (18) 
to Schwarzschild's metric, while the terms gag, Sags--+) Zap» are 


known respectively as quadrupole, octupole, ..., 2° -pole wave 
solutions in approximations of ls-orders, This corresponds to the 
definition of the 2° -pole moment of an axially symmetric system of 
sources distributed linearly along the axis of symmetry 


Q® (u) = mah (u), (1.21) 


where u=t—r, andthecoefficients h®(u) are independent of 
m and a, 

To determine the order of the approximation in which an isolated 
system of sources will exhibit secular variation of mass due to 
emission of gravitational waves, consider a nonstationary system 
which (for a>0) admits of passage to the limit of the stationary 
field of a point source (Schwarzschild field), The simplest system 
of this kind is an axially symmetric distribution of finite length, * 
described (in the spherical coordinates r, 8, 9, t) by the metric 


ds? = —A dr* — r* (B d0* + C sin? 0 dq) + D dt’, (1.22) 


where A,8B,C,D are functions of r,0,¢. Writing down the field 
equations (1,20) in the ps -approximation for this metric and 


* An example of such a system would be two equal point masses, connected by a spring and 
describing symmetric oscillations. 


1. APPROXIMATION METHODS 9 


integrating them (Rosen and Shamir /19/, Bonnor /17/), we obtain: 


da P = VM, + 4M) dt — 
— Vf(Ly + ZL) —§ (Naa + *N,) at} d — (my +49) + 
+ (Gu + 1716,) dB — (x, + rx), 


C=—A cosec? a { [24+ ril{24 + r(\u dt + n)t ar + (1.23) 
+ r4]sin 6 cos 8 d8-+ cosec? 6 \( Ndt+ o\sin*8 d8-+ucosec? 6, 


_— c+ri({24 a r (\Mae+ n) tar + rr, 
D=A +rffared + rfl (L—( Wat —oy)a0 + x}dr-+ry, 


where (J is the d'Alembertian in the coordinates r, 0, 9, t; 

P, M, L and N are the right-hand (nonlinear) sides of equations 
(1.20), assumed to be known from the gr -approximation (¢<p-— 1, 
r<s): 


P=¥,,, M=Y¥y, L=Vy», N = Pao; 


for the remaining values of the indices a and f equations (1.20) 
become identical; lastly, n(r, 9), o (r, 8), x(7, 2, T(8, #), v (8, A, p (r, 2) 
are six integration constants, by choosing which one can satisfy the 
requirement of flatness at infinity and ensure the absence of singu- 
larity of the metric along the axis of symmetry. The subscript 1 in 
equations (1,23) designates differentiation with respect to r. 

Thus the solution of the equations of the ps -approximation for 
the metric (1.22) reduces to integration of an inhomogeneous wave 


( pa) 
equation for the functions A , after which the expressions for 


(ps) (ps) (ps) 
C, 'B and D follow automatically from the remaining relations in 


the system (1.23). In particular, for a linear approximation of 
ls-order the functions P,M,N and Lare equal to zero in view of 
(1.20) and we obtain the homogeneous wave equation examined 
earlier for the functions YY ,,. 

The change in energy of the system of sources in the ps -approx- 
imation is evaluated by Bondi's method /20/, which is based on 


(ps) 
expansion of the ggg in inverse powers of r: 


l . 
1r-"80,u)  (w=t—n). (1,24) 


n=1 


10 GRAVITATIONAL WAVES IN EINSTEIN'S THEORY 


The gravitational mass of a system can be evaluated by establishing 
a correspondence between a given stationary solution and the equi- 
valent field which would be set up by a certain Schwarzschild mass; 


(ps) 
inthe terms W,, it is sufficient to retain coefficients of order 1/r, 


(ps) 
while in the terms Gy of equations (1.20) it is sufficient to retain 
coefficients of order no higher than 1/r°, One should take into 


account only such terms as describe the secular variation of the 


(ps) 
state of the system over the period At, ignoring the terms gas, 


which do not change form as a result of oscillation. Thus all 


terms bes describe only the constant (unvarying throughout the 
oscillation period) component of the field, i.e., the approximation 
of pO-order to the (strictly) Schwarzschild metric (to the 

field of the central mass m for which the linear dimension a =0). 
From this it iS easy to see, in particular, why expansions in the 
parameter m or i’ only do not describe gravitational waves. 


Nonlinear terms of the order W, may contain only combinations 
of the form 


(41) (20) 
Sap * Sap 


and thus vanish owing to the absence of the dipole terms 


Ce = 0). It is therefore sufficient to apply expansion (1.24) to 


terms of the order of oe with s>2, exclusively. But Rotenberg 
has shown that, for isolated axially symmetric systems of sources, 
energy and momentum, determined by the pseudotensor 7%, are 
conserved /21/* inthe 1s-order. The absence of secular mass 
variation in 22- and 23-orders is easy to show by expanding the 
corresponding les and ha in series such as (1.24) and evaluating 


(22) (23) 
W,gand Weg in the 1/r approximation. 


Thus the 24-order approximation is the lowest one which can 
contribute to the secular variation of mass of a system. A rigorous 
solution of the equations for the 24-approximation has been given 


* On the other hand already in the (1s) approximation it can be shown that an isolated system of 
sources loses angular momentum due to the emission of gravitational waves (Campbell /22/). 
For comparison we might note that, in the 1s -order, a rod rotating about an axis of symmetry 
orthogonal to itself — unlike the system under consideration — loses not only energy due to 
gravitational emission (Eddington /23/) but also momentum (Cooperstock /24/, Rotenberg 
/25/), in addition to losing angular momentum due to octupole gravitational waves (Cooper- 
stock and Booth /26/). : 


1. APPROXIMATION METHODS tl 


(24) 
by Hunter and Rotenberg /27/, The nonlinear terms ¥,, should 
contain combinations of the form 


(40) (44) (41) (43) (42) (12) 
Bap * Saps Sap ° Saps 8ap * Sap, 


of which only the last (quadrupole-quadrupole emission) can make a 
secular contribution to the mass variation. Indeed, the second 


(aL) . 
combination vanishes identically ( gag = 0), while for the first 


(24) 
(monopole-2*-pole emission) the right-hand sides ¥,,o0f equations 
(1.20) will not contain secular terms of order 1/r. For quadrupole- 


(24) 
quadrupole emission, onthe other hand, thefunctions ¥,, already 
contain secular terms of order 1/r, causing the change Am in the 
Schwarzschild mass over the oscillation period: 


Am= —Z(a— a)\ (o” (u)}? du; . (1.25) 


here 0 (u) is the quadrupole moment (1.21) of the system, the 
primes denote differentiation with respect to the parameter 
u=t—r, and a and £ are parameters characterizing the nonsta- 
tionary character of the system (a +f =1). From (1.25) it follows 
that the system fails to emit energy only in the strictly stationary 
case a=fp= 15, or for oscillations of a very special type charac- 


terized by the condition on (u) =0. Asit happens, the mass-energy 
losses calculated from this formula correspond exactly to those 
obtained with the help of the energy-momentum pseudotensor in the 
linear approximation (Rotenberg /28/)*, 


3. Critique of the approximation methods 


A significant drawback of the foregoing methods of approximate 
analysis of wavelike gravitational fields (Bonnor-Rotenberg, 


* Bonnor's method of double-parameter approximations can be extended as-well to the case of a 
honempty space-time. Thus Rotenberg has obtained a solution of the ps-approximation equa- 
tions for spaces filled with electromagnetic radiation /29/. The corresponding inhomogeneous 
wave equation in the ps-approximation for the metric (1.22) describes the gravitational and 
electromagnetic radiation. Bonnor's method has also been extended to isolated axially sym- 
metric fields of gravitation of a more general form than the distribution with the metric (1.22) 
(Rotenberg /21/). However, it is more convenient to deal with gravitational emission of 
arbitrarily isolated axially symmetric systems using the general method suggested by Bondi, 
Metzner and Van der Burg (see Ch.11). 
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Einstein-Infeld-Hoffman, Fock, Bondi) is that in the zeroeth 


approximation they deal with the flat space-time metric ae which 
corresponds toabsenceofa gravitational field. The correction to 
the metric plays the part of an infinitesimal of specified order and 
therefore can describe only weak gravitational fields. Thus these 
methods make it possible to determine the state of weak gravita- 
tional emission only, and only against the background of a flat 
Space-time. 

Recent astronomical observations, however, point to the 
possible existence in the universe ofsources ofvery strong gravita- 
tional waves. Thus according to computations by Thorne /30/, a 
neutron star undergoing nonspherical oscillatory perturbations 
(pulsations) maylose, by emission of gravitational waves, anenergy 
of the order of 10° ergs (0.1% of the rest mass of the star itself) 
within a single period of pulsation (10-4107 sec), Similar power 
estimates have been obtained by Weber /31/ for pulsars, Cooper- 
stock /32, 33/ for quasars and binary stars, Zel'dovich and 
Novikov /34/ for collapsed stars, Carmeli /35/ for gravitational 
bremsstrahlung of the sun and Wheeler /36/ for the metagalaxy. 
Moreover, inactual fields of gravitation one should expect that, in 
addition to the gravitational waves causing strong perturbation of 
the space-time metric, the background metric itself may respond 
to the strong gravitational field. 

An important step towards the resolution of this difficulty was 
taken by Isaacson /37, 38/ when he applied the Brill-Hartle method 
of approximations /39/ to the description of so-called "high- 
frequency’ gravitational waves propagating against the background 
ofa highly curved space-time. The idea of the method is based on 
expansion of the quantities h,,, i.e., the correction to the background 
metric y., in series in powers of the ratio ¢=//L, where I is the 
characteristic dimension of the wave perturbation (which can be 
interpreted as its ''wavelength'') and L the characteristic dimension 
of the gravitational field of the background, comparable condition- 
ally to the "radius of curvature’ of the space-time background, * 
This method is an extension of Bonnor's method. Since e is of the 
order of (l?m/r’)/;, the gravitational emission of an isolated system 
at large distances rfrom this system can always be regarded as 
"high-frequency emission,'' corresponding to the case where L—> o, 

According to the rigorous definition, a space-time V, with the 
metric gw describes high-frequency gravitational waves if it admits 
of a one-parameter family of coordinate systems (in the sense of a 


* Indeed, the smallness of the ratio e = U/L is interpreted as the relative shortwave character 
(high frequency) of the traveling perturbation 4A,,, i-e., waves on the background of the curved 
space-time. 
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one-parameter Lie group G with the parameter e) in which gy 
assumes the form 


Sv (©) = Tp () + ehy, (x, 2) = (e<1), (1.26) 

where 
tw =O(1), Oaty =O(1), Papp» = O(A), (1.27) 
hy = O(1), aly = O(8), Aaah == O(e~%). (1.28) 


The conditions (1.27) mean that the curvature of the background 
metric Y» is of the order of unity, i.e., the background space- 
time is curved and its Riemann tensor 


RQs — Rapys(Yys) = (0), 


However, it follows from conditions (1,28) that the order of 
magnitude of the total curvature of V; may even exceed the curva- 
ture of the background: 


Rapys (Yps + thy) = RY: + eR Sos + eRe s See (1 ¢ 29) 
here 
1 
Rvs = 9 (Ray 28 +f hes; ay = hey;a8 — has: ay + RQyshg —R soba )- 


In the above the covariant differentiation is with respect to the back- 
ground metric Yu, which is also used to raise and lower the indices, 
The second term in expansion (1.29), 


eRigrs = O(e"), 


dominates over the other terms in order of magnitude, since, for 
example, 


e® RE = O(4), e®RY . zw O(e), 


and so forth. Similarly, for the Ricci tensor Ry, of the total 
metric we obtain the expansion 


Rap (Yur + ehys) = Rog + eRog + e*Rep t+: 
where 


4 
RY = Rag (Yyv)s Ro = a ik (hoz:ap a hap;or . hra:pp = hrp:ap)s 
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the dominant term being again the second term 


eRQ = =O (e7), 
whereas 
e®R& ~ O(A). 


Hence the vacuum field equations for the first-order approximation 
will be given by 


R&B = 0, (1,30) 
those for the second-order approximation will be given by 
RQ = —2RR, (1.31) 


and so forth, Introducing the quantities 
“ 4 x 7 
Pus cs huy er Yuh, ~ = 1 paps 


where h = y“h,g, we can express equations (1.30) for the first 
approximation in the form 


. 4 z = bs 
buns — > Tw pap — Passe — Prue + 22P orp? + 
4 ROP + ROP. = 0, (1.32) 


where the covariant differentiation is again with respect to Yw. 
For the case of a weak gravitational field, given by the metric (1.2), 
the functions ,, transform into the p,, obtained earlier. 

In order to reduce equation (1,32) to the standard form of a wave 
equation for the potentials $,,, it is necessary to show that the 
second, third and fourth terms in this equation can be removed by 
gauge’ transformations, i.e., by choosing a suitable generator & 
of the Lie group G@ which will induce admissible transformations 


uy aa Ruy = Susy 2. aan 


and will not alter the first term in (1,32), 

To prove this we shall make use of formulas expressing the 
relation between the tensors Rag and Ragy and the barred tensors 
/40/ Rap and Ras 


H(i i 
RY } = RY ) — £: RY, Ros = Rs a LeRQ 85 
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where £¢. represents the Lie derivative in the direction of the 
generator & of the group G. Using the formulas for the Lie deriva- 
tives of the tensors Rag and Rasy (see, for instance, /41/), we can 
show that the quantities ¢RQ and ¢,R\, are of the second order 

in e, i.e., in the high-frequency approximation (e-+0) the tensors 


RY and RGs remain invariant under transformations of the group @ 
(are ''gauge-invariant | /37/). Using the law of transformation of 


the quantities Wy; a and p under infinitesimal transformations of the 
group) G (x* —> x +- £&4), 


Pura —* Pisa — 1 Eve i Roe, 


poh t 2b, 
we choose & so as to satisfy the system of equations 
a +. 1 - 
Y "basa — Rye" =e be =— z Pa. 


Then, obviously, the following conditions will be fulfilled in the 
new coordinate system: 


ee =0, p= 0, 
bringing equations (1.32) to the desired form 


= Ady = 1 Guna + 2RGaP? + RDP + Rp =O. (1.33) 


(00) 
In the case of a flat background metric, yy = g,; due to the 


equality Rina =0allterms apart from y*py;2s vanish and we 
recover the wave equation considered earlier, (1.12), For the non- 
flat metric y. the left-hand side of equation (1.33) is de Rahm's 
generalized topological d'Alembertian /42/, which lies at the base 
of the definition of wavelike gravitational fields* proposed by 
Lichnerowicz /43/. 

Thus gravitational emission on the curvéd space-time back- 
ground finds an elegant definition at the high-frequency limit. 

Isaacson's approach resolves two significant difficulties that 
arise when searching for a definition of gravitational waves valid 
in an approximation of specified order. The first difficulty, as we 
saw earlier, arises from the assumption that the gravitational field 
set up by the wave perturbation is weak, The second arises from 
the assumption of flatness of the background metric. This restricts 


* The rigorous approach to the definition of gravitational waves on the basis of de Rahm's operator 
is considered in Ch-8. 
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the validity of the approximation methods to asymptotically flat 
gravitational fields only (in particular, to fields of isolated sources). 

However, the Isaacson-Brill-Hartle method, together with the 
other methods described above, fail to resolve a third difficulty 
common to all attempts of defining gravitational waves on the basis 
of approximation methods: there exists no proof that the series of 
successive approximations converges. Thus Bonnor's method is 
based on the assumption of convergence of the series (1.19) in the 
neighborhood of m=0 and a=0, which, however, does not guaran- 
tee its convergence at sufficiently large distances from the source 
system. As Fock has shown /9/, upon expansion of the integrand 
in (1.13) for the wave function , sufficient convergence of the 
series is guaranteed only for ''moderately large" distances r, 
namely for distances large compared with the dimensions of the 
source system but small compared with the length of the emitted 
waves. Inside the ''wave zone" (i.e., the region situated at at a large 
distance from the sources compared with the length of the emitted 
waves) convergence of the series cannot, in general, be guaranteed 
unequivocally. Yet this very case constitutes the region of appli- 
cability of the Isaacson-Brill-Hartle method. This critical situa- 
tion is further aggravated by the absence of reliable experimental 
data concerning the properties of gravitational waves, whence the 
particular urgency of a rigorous proof that the series of successive 
approximations does indeed converge inside the wave zone, 

Lastly, every one of the methods of approximation considered 
above presupposes the selection of a definite system of coordinates 
(oraclassofadmissible coordinate systems) the existence of which 
cannot be guaranteed a priori in the specified gravitational field. 
For example, when the gravitational emission of an isolated 
axially symmetric system of bodies is treated by Bondi's method, 
the harmonic coordinate system becomes inapplicable, although it 
was applied successfully by V. A. Fock to the description of gravi- 
tational waves by the method of expansion of the potentials in the 
parameter 1/c. The reason for this is the appearance of the loga- 
rithmic term r log r in place of 1/r in the expansion of the gravi- 
tational potentials, precluding passage to the limit of the Schwarz- 
schild metric (Bonnor /11/), Isaacson and Winicour /44/). Conse- 
quently, for fairly large distances from the isolated source of 
gravitational emission the wave solution of the linearized theory of 
gravitation in harmonic coordinates cannot serve as a first approx- 
imation to the exact wave solution, although at distances signifi- 
cantly smaller than the length of the propagating wave this des- 
cription proves satisfactory /45, 46/. 

At the same time, the equations for the gravitational field are 
generally covariant. Therefore every physical consequence of the 
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theory must admit ofa generally covariant formulation. Hence 
our next task will be to consider rigorous generally covariant 
methods of describing gravitational waves. Our starting point will 
be the Cauchy problem for Hinstein's equations of gravitation. 

Henceforth it will be convenient to distinguish between the term 
gravitational waves andtheterm gravitational 
emission, With the former we will associate the field in the 
true wave zone, where it does not interact with the sources, The 
latter will mean the overall gravitational field of the source which 
is inducing the waves. In the case of an empty space-time, we 
will refer to the gravitational waves as free. 


Chapter 2 


THE CAUCHY PROBLEM FOR EINSTEIN'S EQUATIONS 


1. Einstein's equations as a system of the hyperbolic type 


Rigorous formulation of the problem of gravitational waves in 
the general theory of relativity became possible only after De 
Donder /47/ and Lanczos /48/ had proved that Einstein's system 
of equations was of the hyperbolic type,i.e., that its characteristics 
were identical with those of a wave equation of the form* 


Lh = as On (V = ge*00¥)=V. (2:1) 


Indeed, Einstein's equations in empty space 


can be reduced by identity transformation (see, for instance, /9/) 
to the form 


1 “ 
> 8 Buia + Ty — Ly = 0, (2.3) 
where we set 
Tuy =S (uy + Tu) — Til oy (2.4) 
le fel) FP Se" as (2.5) 


so that the ZL, are expressed only in terms of components of the 
metric tensor and their first derivatives: 


* Henceforth the d'Alembertian symbol will be Interpreted in the sense of (2.1). 


18 
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Luv = g%8 go BL ev,0 = g% 2POT oD yap. (2 .6) 


4 


We now pass in the given space-time region to a special coordi- 
nate system, choosing four solutions of equation (2.1) as the new 2s: 


S os 1 6; 
or=P= = (V—£8"),9 = 0. (2.7) 


It is customary to refer to such systems as the above as 
harmonic coordinate systems. Inthis system equations 
(2.2) become 


4 
8 Buy, — Lyy = 0. (2.8) 


It is well known that the characteristics of a system of partial 
differential equations are determined exclusively by the coefficients 
of the higher derivatives. In our case it is the components of the 
metric fensor g#® themselves that play this part. But the g 
matrix can be reduced to the canonical form characterized by the 
signature (+, —, —, — ) by a nonsingular transformation of the 
coordinates at any point in space-time. This proves that the 
system of quasilinear equations (2.2) is hyperbolic (see /49/, p. 61). 


2. Hadamard's discontinuity 


In the classical theory of partial differential equations (see, for 
instance, /50/), wave propagation in space is characterized by 
Hadamard's discontinuity inthe solution of the equations 
on the initial hypersurface S. We will show below that the hyper- 
surface of discontinuity S of the field functions (of their deriva- 
tives), known as the wave-front surface, is the character- 
istic hypersurface of the field equations. Our next task is there- 
fore to determine the characteristic manifolds ("characteristics") 
of Einstein's equations. However, before applying the concept of 
Hadamard's discontinuity directly to Einstein's equations, we will 
first illustrate it with the help of the scalar equation (2.1). 

Assume that the function is continuous in each of the regions 
1 and 2 into which the surface S divides the space-time region 
under consideration, and that it tends to the limits p and yp as 
attends to a point P, (z@) on S from the regions 1 and 2, 


20 GRAVITATIONAL WAVES IN EINSTEIN'S THEORY 


respectively. Then the following function of the point P, will be 
called the Hadamard discontinuity of the function on the surface 8S: 


[pl (Po) = vo} — vy. (2.9) 


Now let the function be everywhere continuous near § while 
certain of its first derivatives p, have finite discontinuities on S: 


[y] =0, [hal #0 (2.10) 


We construct the total differentials of ,,) and ,) on S: 


Api == Puy, ade", — APay = Pa, ada. 


Their existence and continuity was demonstrated by Hadamard /51/ 
by considering the limits to S from regions 1 and 2, Subtracting 
the one from the other, we obtain, due to (2.10) and the continuity 
of dp on S: 


(Ya,a — Ye, «) @2* = [pa] dat = 0. (2.11) 


Let the surface S be given by g(z*)=0. For 0.9 — the normal 
to the surface § —the relation 


P,adx* = 0 (2.12) 


holds if the increment dz* belongs to the surface §. Comparing 
(2.11) and (2.12), we conclude that [p,] and 9, are proportional to 
each other: 


ltpal = XP, (2.13) 


If the first derivatives ofthe function are continuous, it can be 
shown in a similar fashion that the discontinuities of the second 
derivatives are expressed by the formula 


Mapl = X09,» (2.14) 


and so on (cf. Hadamard /51/, pp. 81 — 89). 

Thus the problem of the investigation of gravitational waves as 
solutions of Einstein's equations must be linked to nonanalytic solu- 
tions of the Cauchy problem for a system of quasilinear’ hyperbolic 
partial differential equations; in other words, it is to be presumed 
that the coefficients of the equations, the initial data and the 
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solution itself may be functions of a finite order of smoothness C’ 
while derivatives of these functions of orders higher than r may 
have Hadamard discontinuities on certain hypersurfaces. (The 
function p is said to be of class C” if it has continuous partial 
derivatives up to the order r inclusive.) If » is of class C™1 in 
the neighborhood of S while its r-th derivatives have Hadamard 
discontinuities on S, is said to be a piecewise smooth function 
of class C"’, or piecewise* C”, 

The solution of Cauchy's problem will depend not only on the 
order of smoothness of the functions in which the initial data are 
given and the solution is sought, but also on the character of the 
initial hypersurface S. More specifically, the solutions to the 
problem will be substantially different depending on whether the 
hypersurface Sis free or characteristic. ** 

For Einstein's equations (2.2), the important case from our 
standpoint, Cauchy's problem may be stated as follows: 

Let the initial hypersurface § be given by an 
equation of the form 


S: 9 (z*) = 0, (2.15) 


and let the functions gyg(x*?) and their first deriva- 
tives gup(x?) be given on it; determine functions 
gap (x) and their first derivatives which reduce to 
the initial data on S and satisfy equations (2.2) 
for gap Outside S. ; 

In Lichnerowicz's analysis of the Cauchy problem for Einstein's 
equations, the choice of a convenient coordinate system played an 
important role. We will henceforth assume that the coordinate 
system is harmonic so that equations (2.2) assume the form of 
(2.8) while the initial data obey the condition 


r’ =0. 


The exact solution of Cauchy's problem for Einstein's equations in 
harmonic coordinates has been given by Lichnerowicz in /54/, 
which we will be following henceforth. 

Within the limits of the condition of harmonicity throughout the 
space-time region inside which the solution of Cauchy's problem is 


* "Par morceaux" in Lichnerowicz's original terminology (see /52/, pp.27 — 35). 
** A general formulation of Cauchy's problem and definitions of free and characteristic hyper- 
surfaces will be found in the monographs of Petrovskii /49/ and Bers, John and Schechter /53/. 
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being sought, it is possible to choose local coordinates z* on S 
such that in these coordinates the equation for S reduces to 


p(z*) = 29 = 0. (2.16) 


It is easily seen from the formulation of Cauchy's problem that 
the solutions of interest for the problem of gravitational waves are 
functions of class C! (piecewise C*) for sufficiently smooth (but 
nonanalytic) initial data. In other words, g,, and their first deri- 
vatives g,> will be continuous on the hypersurface S, while cer- 
tain second derivatives g,. will have Hadamard discontinuities on 
S. Let us determine which of the second derivatives of g,, can 
have a discontinuity on the hypersurface S given by equation (2.16). 
According to Hadamard's formula (2.14), discontinuities of second 
derivatives can be expressed as 


[Zuv,2g] = AvP,oP,p, (2.1 7) 


where a,, are the so-called "coefficients of discontinuity" (see also 
the works of Trautman /55/ and Bel /56/). From this it follows 
that of all the second derivatives of g,,on S, only yo can have a 
discontinuity: 


Cee (2.18) 


But ofthese, the derivatives of 824do not appear in Einstein's equa- 
tions. Therefore the important point for us is that ofall the derivatives 
of g,, entering into equation (2.2), only second derivatives of the form 
Zijoo canbe discontinuous inthe hypersurface §. Inthis case Cauchy's 
data for the problem (2.2) reduce to specificationon S only ofthe values 
ofthefunctions g,, and gw,»9, which we assume tobe thrice and twice con- 
tinuously differentiable with respecttothe z' coordinates, respective- 
ly. The remaining first derivatives of g,,, as wellas the second deriva - 
tives of go, and go, are uniquely determined by differentiation with re- 
spect to z'of the previously knownCauchy dataon S., 

The following question arises in this connection: when do 
Einstein's equations (2.2), together with the Cauchy data, uniquely 
determine the derivatives 8ij,0 as well? 

As Lichnerowicz shows (/52/, pp. 31 — 33; see also /57/, p.264), 
the system (2.2) in the harmonic coordinate system is equivalent to 
a system of the form 


: 
8 Bij 00 + 243 = 9, (2,19) 
S=0, r=, (2.20) 


2. THE CAUCHY PROBLEM 23 


where 
Sh = Ri — + ROG; (2,21) 


here Q,;, S, and [* do not contain derivatives of the form ij, and 
hence are completely determined by the Cauchy data, Fulfillment 
of conditions (2.20) on the hypersurface § guarantees their fulfill- 
ment in the neighborhood of § as well, provided equations (2.19) 
hold on § and outside it (Lichnerowicz's theorem regarding 
involution of the system of differential equations /52/). 

This means that conditions (2.20) serve only for determining the 
Cauchy data, which cannot, therefore, be arbitrary. At the same 
time, system (2.19) serves for integration of the system (2,2) over 
x, i.e., for determining the unknown functions g,,. Thus in empty 
space Cauchy's problem for Einstein's equations breaks down into 
two parts: 1) determination of the Cauchy data that satisfy 
equations (2.2); 2) integration of the system of equations (2,19) 
over 2°, . 

Let us assume that we have Cauchy data satisfying the conditions 
(2.20). Then under the condition g® +0 equations (2.2) admit ofa 
solution which does not have the Hadamard discontinuity on S; 
this corresponds to the case where equations (2.19) taken together 
with the Cauchy data uniquely determine the derivatives &ij,00 as 
well, On the other hand, if g°° =0 in the neighborhood of S, then 
the second derivatives ij 0 — and therefore also the curvature 
tensor components Roig3— Cannot be uniquely determined on S by 
the Cauchy data and the field equations, i.e., they have the 
Hadamard discontinuity on S, * 


3, Characteristic hypersurfaces of Einstein's equations 


The condition g® =0 governing the Hadamard discontinuity for 
the curvature tensor on the initial hypersurface S, specified in the 
form (2.16), can be formulated in generally covariant form. Indeed, 
passing to arbitrary coordinates in which 2° = 9 (zx), let us write 


* By definition, discontinuity of the derivatives of functions expressing the solution of Cauchy's 
problem for a certain system of equations subject to continuity of the solution itself is known 
as weak discontinuity of the solution of the given system of equations. Thus the 
condition g® =0 denotes a weak discontinuity of the first order for the solution of the system 
of equations (2.2). 
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down the condition 


Ox° 0x0 reac} _ 0 
az’* Oz'8 a 


go — 


in the form of a general equation for the hypersurface S admitting 
a weak discontinuity in the solution of the system (2.2) (primes 
have been omitted): 


g8Q,aP,p = 0, (2,22) 


i.e., the eikonal equation from geometrical optics (/16/, 
p. 173). 

But equation (2.22) is a necessary and sufficient condition for 
isotropy of the hypersurfaces S (see /58/, p.57). Consequently, 
first-order weak discontinuity in the solution of the system of 
equations (2.2) (Hadamard's discontinuity in the curvature tensor) 
is possible only if the initial hypersurface Sis isotropic. 

We see that the solution of Cauchy's problem for Einstein's 
equations in empty space depends radically on the nature of the 
initial hypersurface §. If the hypersurface S specified by the 
equation g (2*) = 0 does not satisfy the condition (2.22), 


g2Q oP 2 0, (2.23) 


it is called free, and Cauchy's problem for equations (2.2) in 
functions of class C! (piecewise C*) admits of a unique solution. But 
if S does satisfy condition (2.22) itis termed characteristic, 
and the uniqueness theorem of the solution of Cauchy's problem 
does not hold. 


4, Leray's theorem 


It is also possible to tell something of the solution of Cauchy's 
problem from the general existence and uniqueness theorems fora solu- 
tion ofa system of quasilinear hyperbolic partial differential equations. 
Thus if the coefficients entering into the equations are all analytic 
and the solution is sought in analytic functions, then given analytic 
initial data ona free hypersurface the Cauchy problem has a solu- 
tion and moreover a unique one (Cauchy~Kowalevski theorem; see 
/53/, p.56). However, in the case of interest to us — the solution 
of Cauchy's problem in C!(piecewise C*) functions — the Cauchy- 
Kowalevski theorem is inapplicable. 


2. THE CAUCHY PROBLEM 25 


As Lichnerowicz points out /54/, in our case it is possible to 
use Leray's theorem /59/ according to which a system of quasi- 
linear hyperbolic partial differential equations will have a nonana- 
lytic and moreover a unique solution provided that the initial 
hypersurface S is free (spacelike) and that the Cauchy data on this 
hypersurface are given by sufficiently smooth functions, With 
regard to equations (2.2) all the conditions of Leray's theorem are 
satisfied if one uses harmonic coordinates in which Einstein's 
equations assume the form (2.8) and the initial data are determined 
by the solution of the system of equations (2.20). 

The following important result follows from Leray's theorem: 

The Riemann tensor Rags in empty space-time 
can have an Hadamard discontinuity only on the 
characteristic hypersurface S of Einstein's equa- 
tions (2.2), which is determined by the eikonal 
equation (2.22), 

Indeed, the presence of an Hadamard discontinuity ina Riemann 
tensor on a certain hypersurface S§ is evidence that some at least = 
of the second derivatives 8 ,.. cannot be determined uniquely from 
the field equations in the neighborhood of S from the initial data 
specified upon it. This in turn means that the hypersurface 
cannot be free. 


5. Bicharacteristics of the equations of gravitation 


We saw that in the Riemannian space V4 of signature —2 the 
characteristic manifold of Einstein's equations in the absence of 
matter 


Ry =0 (2.24) 


is the isotropic three-dimensional hypersurface Vs; (2.15), the 
function » satisfying the eikonal equation (2.22), 

The envelope of the tangent planes to the characteristics passing 
through a given point is called the characteristic cone /50/ at this 
point. And as the characteristic hypersurface of Einstein's equa- 
tions is isotropic (i.e., the metric is degenerate), the character- 
istic cone of the system of equations (2.2) coincides with the light 
cone at the given point (/52/, pp.33—35), Equation (2.22) asa 
characteristic equation for the system (2.1) was first considered 
by Finzi /59/. 
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According to the definition of bicharacteristics (also 
known as rays) fora system of quasilinear second-order equa- 
tions (see /50/, p. 551), the bicharacteristics of Einstein's equations 
(2.2) are identical with the flow lines of the isotropic vector field I 
orthogonal to the characteristic hypersurface S, 


In = g8 (pp, (2.25) 


and may be expressed by the equations 


a = 8° Q,a5 (2,26) 


where t is a parameter on the curve. From (2.26) it obviously 
also follows that 


lg = Spa =~ . (2.27) 


The functions z*(t) and /,(t) can be specified by the canonical 
system of equations 


dx OH a, aH 
dy Oly,” a tC® 42528) 
where the Hamiltonian 
H (2*,lg) =~ gapvl? (2.29) 


is identical with the characteristic form of Einstein's equations in 
empty space /50/. Both the solutions 2*(t) of the canonical 
system (2,28) determine the extremals of the Lagrangian 


=o 808 Ge ae 


oh 
In actual fact, passing from the variables (z=, e) to the variables 
(x*, 18), we obtain the classical relation between L and H: 
dx* OL 
a= (Fe) 


dx* 
a( z Jax 
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Since the first integral of the system (2.28) is 2L= C = const, 
these solutions will be extremals for 


az dxB 
7 = a 
Vi 898 Ode 


as well, i.e., they will be geodesics of the Riemannian space V, 
with metric gap (/52/, pp. 33—35). 

As we know from the theory of partial differential equations /50/, 
bicharacteristics belong to the characteristic surface, i.e., the 
tangents to them are the generators of the characteristic cone, 

But, as we saw earlier, the characteristic cone for Einstein's 
equations is identical with the light cone. Consequently, the 
bicharacteristics of Einstein's equations are iso- 
tropic geodesics. 


6. The Cauchy problem for the Einstein-Maxwell equations 


In analyzing the definition of gravitational waves it may be of 
considerable benefit to make use of the experience acquired in the 
investigation of electromagnetic waves, which are described by 
the energy-momentum tensor 


Tap = 7 6apF iP — FEF pp, (2,30) 


where Fy, is the Maxwell stress tensor. It is therefore useful to 
examine Cauchy's problem for the Maxwell equations in the Rie- 
mannian space-time of the general theory of relativity. The 
general solution of this problem has been investigated by Lich- 
nerowicz (/52/, pp.43—52), 

The self-consistent system of Einstein-Maxwell equations has 
the form 


Da = 8 Fayip = 0, (2.31) 
4 
Et =~ nPOoF gg = 0, (2732) 
Qup = Rap — 5 Roop + tag = 0, (2.33) 


where 778° is the discriminant tensor* 


* More precisely axial tensor, since under coordinate transformation its sign will depend on 
the sign of the Jacobian of the transformation (see, for instance, /60/, p.11). 
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nrbrd = = erbyd | Nays = — V — 8 &apys, (2. 34) 


and euws is the wholly antisymmetric Levi-Civita symbol, equal to 
+1 for e.5, and all even permutations of the indices, to -—1 for all 
odd permutations and to 0 for the remaining cases. 

Let us assume that the field F,, is of class C° (piecewise C*). 
Then Cauchy's problem for the Einstein-Maxwell equations (2.31) — 
(2°33) may be formulated as follows: let the gravitational and 
electromagnetic fields be specified on the initialhypersurface (2.15); 
find the fields outside S ifthey satisfy equations (2.31) — (2.33). 

In the coordinate system in which the equation for S has the 
form (2.16), the Cauchy's data consist of gaa (continuously differ- 
entiable with respect to z' at least three times), go, and Fup 
(twice differentiable with respect to these coordinates). Only the deri- 
vatives 8apo0and F,,, maybediscontinuouson S. Separating the in- 
dex 0 for convenience's sake, we can write equations (2.31) inthe form 


D,= 8° Foi sg 2 Frio +...=0, (2,35) 
Dy = 8"Fing + ...=0 (2.36) 


(henceforth the three dots are employed so as to avoid writing out 
explicitly those terms which are uniquely determined by the Cauchy 
data). 

Lichnerowicz shows that the system (2,35) — (2.36) is equivalent 
to the system comprising equations (2.35) and 


D®=0, (2.37) 


where the quantity D®° (=g'D; + gD), as is evident from equations 
(2.35) — (2.36), is uniquely determined by the Cauchy data. Simi- 
larly, equations (2,32) are equivalent to the system 


{ ne?! F egg +... = 0; (2.38) 


fF 41,3 = 0, (2.39) 


where E° is determined by the Cauchy data. Lastly, equations 
(2.33) are equivalent to the system 


Qu = Ry— > Rey + dry =0, (2,40) 


Q2 = RY — +- R68 + ark =0, (2.41) 
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where Q; is given by the Cauchy data. (We recall that outside of 

the sources the trace of the energy-momentum tensor t; =0.) Thus 
Cauchy's problem may be broken down into two parts: 1) determina- 
tion of the Cauchy data satisfying conditions (2.37), (2.39) and 

(2.41) on §; 2) integration of the system of equations (2.35), (2.38) 
and (2,40) with respect to the coordinate #° for the obtained initial 
data. 

We assume from the first that g°=—0 throughout S. Then Fiz» 
is uniquely determined from equation (2.38), Fo; is uniquely deter- 
mined from (2.35) and gijoo is uniquely determined from (2.40), Thusall 
the second derivatives of gigand all the first derivatives of Fy, are 
uniquely determined by the initial data and initialsystem, and there- 
fore the latter has a unique solution, 

Now let g® =0 0n S. Then ij, and Foi. may be discontinuous 
on § and consequently S$ is a characteristic manifold of Maxwell's 
equations. The characteristic equation, as we know, has the form 
(2.22); in other words, the characteristic manifold of Maxwell's 
equations is identical with the characteristic manifold of Einstein's 
equations. 

Thus we have proved Lichnerowicz's theorem (/52/, pp. 50 —52): 

The characteristic manifolds (2.15) of Einstein's 
equations and of Maxwell's equations in VY, are 
identical and are determined by the solutions of 
the eikonal equation (2,22), 

An immediate consequence of this theorem is the fact that the 
bicharacteristics of Einstein's equations are identical with those 
of Maxwell's equations. 


7, Gravitational wave front and "rays" of gravitation 


On the basis of the foregoing it may be stated that the character- 
istic manifold of Einstein's equations is the hypersurface on which 
the Riemann tensor has a Hadamard discontinuity, i.e., this hyper- 
surface plays the role of a wave front understood as the sur- 
face of discontinuity of Ragy /55, 59/. The bicharacteristics of 
Einstein's equations, on the other hand, are the trajectories of an 
isotropic vector orthogonal to the characteristic hypersurface 
(wave front) and therefore play the part of wave vector. Since the 
characteristic manifold and the bicharacteristics are invariants of 
coordinate transformation (/50/, p.555), the three-dimen- 
sional characteristic hypersurface of Einstein's 
equations can be regarded as invariantly deter- 
mined gravitational wave front while the bicharac- 
teristics of Einstein's equations can be regarded 
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as invariantly determined gravitational rays, i.e., 
the trajectories of propagation of the wave front. 
Similarly, the electromagnetic wave front in the space Y, is 
determined by the characteristic hypersurface of Maxwell's equa- 
tions; by Lichnerowicz's theorem it is identical with the 
gravitational wave front. The trajectories of electromagnetic 
wave propagation — the electromagnetic rays — can be defined as 
bicharacteristics of Maxwell's equations; they are identical with 
gravitational rays. 


Chapter 3 


GRAVITATIONAL WAVES: SUBSTANCE OF THE PROBLEM 


1. Various aspects of the problem 


From the standpoint of the results given in Chapter 2, fields 
describing free gravitational waves are determined by the solutions 
of Einstein's equations (2.2) for initial data ona characteristic 
hypersurface. As arule, however, partial solutions of Einstein's 
equations are obtained without specifying definite boundary condi- 
tions, and thus the wave character of a concrete solution may not 
be apparent. Yet when studying gravitational fields we have only 
partial solutions of Einstein's equations at our disposal. Hence 
the problem arises of determining ina generally covariant manner 
the class of gravitational field corresponding to Hadamard's dis- 
continuity in solutions of Einstein's equations with initial data on 
the characteristic manifold. 

No definite solution of this problem exists, despite numerous 
variants which will be reviewed in the following sections, From 
the point of view of the theory of differential equations, the diffi- 
culty of the problem arises from the complexity of the nonlinear 
structure of Einstein's equations and the absence for the latter of 
universal boundary conditions, From the differential-geometric 
standpoint the difficulty resides in the absence of a generally 
covariant d'Alembert operator explicitly following from Einstein's 
equations. Lastly, from the physical standpoint the problem is 
that the absence ofa generally covariant expression for the gravi- 
tational field energy in the general theory of relativity makes it 
difficult to tackle the question of the possibility of energy transfer 
by gravitational waves, as well as to describe the waves themselves 
in terms of the free transfer of field energy. 

Like all physical problems the problem of gravitational waves 
calls for experimental as well as theoretical treatment. In other 
words, gravitational waves can and should be viewed as an aspect 
of physical reality susceptible of experimental measurement. But 
if the experimental data are to be given a correct theoretical inter- 
pretation, the physical theory must contain a few general premises 
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independent of concrete experiments. Only then will the theory be 
whole and closed, i.e., possess an intrinsic logical base (see, for 
example, /61/). Otherwise any comparison between the theory and 
the data of a given experiment would be tautological. 

In this sense the general theory of relativity is exceptional 
among physical theories. The rigorous geometric foundation on 
which Einstein's theory of gravitation rests allows us to hope that 
it will be possible to arrive at a general and rigorous substantia- 
tion of the concept of gravitational radiation. It is to elucidate 
this possibility that we undertook to discuss the Cauchy problem 
for the equations of the gravitational field in the preceding chapter. 
There we found that the solution of Cauchy's problem for Einstein's 
equations of gravitation exhibits profound similarity at a number of 
important points with the solution of the corresponding problem for 
the equations of the electromagnetic field. In particular, the 
characteristic manifolds and bicharacteristics of Einstein's equa- 
tions and those of Maxwell's equations in the space-time V, are 
identical, In classical electrodynamics the characteristic 
manifold of Maxwell's equations describes the electromagnetic 
wave front, and the bicharacteristics of Maxwell's equations 
describe the trajectories of propagation of electromagnetic radia- 
tion. In the light of this analogy one is justified in assuming that a 
gravitational wave front can be defined theoretically in terms of the 
characteristic manifold of Einstein's equation and its propagation 
trajectories in terms of the bicharacteristics. 

The basic problem, however, remains that of defining the field 
of gravitational radiation, or the gravitational wave zone. The 
foregoing considerations do not rule out the possibility of finding a 
satisfactory definition. The search for initial premises could 
clearly proceed along the lines of further pursuit of the analogy 
between the gravitational and the electromagnetic field. It should 
obviously be borne in mind that this analogy cannot be pursued 
indefinitely. Nevertheless all the methods for describing gravita- 
tional radiation discussed below make some use of the analogy 
between the electromagnetic and gravitational fields. 

This analogy is detectable by a variety of methods. The first of 
these follows directly from comparison of the differential structure of 
the equations of gravitation with that of the equations of electro- 
magnetism, i.e., from comparative analysis of the solution of 
Cauchy's problem for Einstein's and Maxwell's equations. This 
analysis is easiest to perform if one takes the curvature tensor of 
the space-time V, to be analogous to the electromagnetic tensor F,, 
in the gravitational equations, Einstein's equations 


4 
Rag = — Mag (Uap = Tap — > TEa9) (3.1) 
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and Bianchi's identities 
Rapysie + Raporss + Ragsox = 0 (3,2) 
yield the relations* /60, 62/ 
Rigess = — WO arp) (20 tas] = Vaya — Upia). (3.3) 


Relations (3.2) and (3,3) exhibit a striking analogy with Maxwell's 
equations 


Pays + Pop + Fron = 9; (3.4) 
Fay = — Iw (325) 

and the tensor 
J yap ae 2M 4fa;9) (3.6) 


can be interpreted as the gravitational analog of the electro- 
magnetic current. However, if (3.2) and (3.3) are regarded as 
equations for the components A,,,s5, a problem arises: under what 
conditions will Einstein's equations (3.1) follow from equations 
(3.2) and (3.3)? The problem can be formulated with greater rigor 
as follows: what conditions must the initial data of equations (3.2) 
and (3.3) satisfy in order for the complete class of their solutions 
to define the set of all gravitational fields of Einstein's theory ? 

This question has been resolved by Lichnerowicz, who shows 
/62/ that if, on an initial hypersurface S oriented in space, the 
initial data of equations (3.2) and (3.3) (i.e., the components of the 
Riemann tensor) satisfy relations (3.1), the latter will also be 
satisfied in the neighborhood of §.** In other words, equations 
(3.2) and (3.3) are completely equivalent to Einstein's equations of 
gravitation if the Cauchy data for equations (3.2) and (3.3) on the 
Spacelike initial hypersurface S are coupled by the conditions (3.1). 
In this case equations (3.2) — (3.3) are called quasi- Maxwell 
equations of gravitation /63/. 

However, the ''quasi-Maxwellian approach" leads to fundamental 
difficulties in the description of gravitational fields of the wave 


* Here and below square brackets denote antisymmtetrization with respect to the indices they 
enclose; correspondingly, parentheses denote symmetrization. 

** This theorem is the analog of the theorem concerning the involution of Einstein's system of 
equations, discussed in the preceding chapter. 
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type which are defined by initial data on a characteristic (isotropic) 
hypersurface. Thus it is possible to construct for the gravitational 
field a wave equation corresponding to the equation for the electro- 
magnetic field tensor F,, which follows from Maxwell's equations 
in Riemannian space-time (Tolman /63/),. However, this analog 
turns out to be an identity for the Riemann tensor Ry,,s and thus 
cannot serve to separate wave fields as a special class of gravita- 
tional field. 

Aside from differential methods, the analogy between the gravi- 
tational and electromagnetic fields can also be detected algebra - 
ically. This makes it possible to construct a definition of fields of 
gravitational waves based on the similarity in algebraic structure 
between the electromagentic field tensor Fy, and the Riemann 
tensor Rays. AS we will show in the next chapter, a characteristic 
trait of fields of electromagnetic radiation which sets them apart 
from the set of all electromagnetic fields can be expressed by 
purely algebraic conditions, namely the vanishing of the invariants 
of the electromagnetic field tensor F,,. From this point of view the 
definition of electromagnetic radiation is based on the algebraic 
subdivision of all electromagnetic fields into two types, interpret- 
able physically as wave and nonwave fields. This suggests the 
application of algebraic classification to the definition of gravita- 
tional wave fields as well. However, here it is immediately seen 
that the algebraic analogy between the electromagnetic and gravita- 
tional fields is not complete. 

Indeed, due to differences in algebraic structure between the 
electromagnetic field tensor F,, and the space-time curvature 
tensor Ray, gravitational fields subdivide not into two types (as 
do electromagnetic fields) but rather into five algebraically distinct 
types, defined by Petrov's classification (three basic types,of which 
two can be both degenerate and nondegenerate). Asa result the 
gravitational wave field possesses diverse algebraic properties and 
correspondingly there is a multiplicity of algebraic criteria for 
distinguishing wavelike gravitational fields from all gravitational 
fields. This is another difficulty accounting for the fact that the 
problem of gravitational waves in the general theory of relativity 
has yet to receive a universally accepted theoretical solution. 
Petrov's algebraic classification of gravitational fields /64, 65/, 
however, has a very important role to play in what follows. It is 
therefore useful to consider it in some detail. 
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2. Algebraic classification of gravitational fields. 
Einstein spaces 


Let a given Riemann space V, be an Einstein space, i.e., let it 
be described by the equations 


Rap = ¥Gap- (3.7) 


It is easily shown that one will then have 


n= +R = const. 
Proceeding according to Petrov's method (/65/, pp.113—117), we 
map the Einstein space at every point on the centered affine bi- 


vector space By of dimensionality 
N=C, = n(n—1)=6 (m = 4), 


placing a single collective index in the space Byin correspondence 
with each antisymmetric pair of indices of an arbitrary tensor in 
the Einstein space. Then to an arbitrary bi-tensor (i.e., tensor 
whose indices subdivide into antisymmetric pairs) belonging to the 
Einstein space there corresponds a tensor ofhalfthe valency in the 
Space Bg. 

We metrize the bi-vector space B,, introducing into it the 
tensor gy(a,b=1, 2, 3, 4, 5, 6) as the image of a fourth-rank 
tensor in the Einstein space: 


Sab > Saprs = Bares — Sad pr» 
(3.8) 
(a8) — a, (78) —> b. 


Assuming that V, has the signature -—2 (i.e., +, -, -, -—) and 
fixing the enumeration of the bi-vector space indices 


10> 1, 20> 2, 30> 3, 23+ 4, 31 +5, 12 >6, (3.9) 


we obtain the canonical form of the metric of the space R, in the 
chosen orthogonal frame 
2 0 
| ai . (3,10) 


35 


36 GRAVITATIONAL WAVES IN EINSTEIN'S THEORY 


where ¢ is a unit 3X3 matrix. From this follows, in particular, 
the nondegeneracy of the matrix | gap]. 

Writing down equations (3.7) in the orthogonal frame we arrive 
at the conclusion that in Einstein spaces the matrix ]R,| of the 
curvature tensor is block-symmetric: 


& E# 

(Ran) = (5¢ 9) (3.11) 
where the blocks ¢é and # are symmetric 3X3 matrices the 
elements of which satisfy the relations 

3 3 
D>) es = —%, Dy gg =O. (3,12) 


s=1 s=1 


Then, decomposing a A-matrix of form | Ra — Ag y|| into two three- 
dimensional complex conjugate matrices, we arrive at Petrov's 
fundamental theorem: there are three and only three types of gravi- 
tational field, defined in R, by the characteristics of the 4-matrix 
of the curvature tensor: 


Type 1 Type 2 Type 3 
(144, 111] [24, 24] [3, 3] (3.13) 
respectively. Here the bar indicates elementary divisors with 
complex conjugate bases. For type 3 the elementary divisors have 
real bases (no bar). All spaces of constant curvature defined by a 
characteristic oftheform [(111, 111)] belongtotype1 (/65/, p.119). 

Petrov shows (/65/, p.$19) that the matrix |R»] of the curva- 
ture tensor in a canonical nonholonomic orthogonal frame is reduc- 
ible to the form (3.11), where for type 1 fields 


a, 0 0 B, 0 0 
#-( ot, o), -(s 8, 0) (3,14) 

0 0 a 0 0 B, 
Die =—%, DB, = 0; (3,15) 


for type 2 fields 


a, 0 0 2, 0 0 
1-(1 dd). °0 ) (0 B, ') (3.16) 
0 0 o—1 0-4 p.7- 


t+ 2a,=—x, 8, + 28, = 0; (3,17) 
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and for type 3 fields 


1 
=| 4 —3y* O |, #=[0 0-4}. (3,18) 
0 OB tote 0-1 0 


Here a, and B, are respectively the real and imaginary parts of 
the bases of the elementary divisors 


Os = Gs -}- iB, (s = 4, 2, 3), (3,19) 


which are identical with the eigenvalues of the matrix || Ral. 

Since there are three types of Einstein spaces, let us agree to 
designate them by *7,, where i=1, 2, 3 denotes the type of 
gravitational field. An empty space-time R,, =0, i.e., the 
space *7; for x= 0, will be denoted by 7;. Petrov's classification 
was subsequently formulated within the frameworks of other 
formalisms used in gravitational wave research. Thus Debever 
gives a detailed description of types and subtypes of fields after 
Petrov which is based on the mutual orientation of isotropic vector 
fields in physical space-time/66/, while Penrose investigates the 
spinor properties of the Riemann tensor from the standpoint of the 
algebraic classification according to Petrov /67/. Penrose 
gives the following diagrammatric representation of Petrov's 
systematics:* 


— TI 


tod 
Roy: (3,20) 


lle V0 
i Ts *"T, 


Here I,D,O are subtypes of *7,, distinguishable by the following 
properties: for I the three eigenvalues in the blocks of the matrix 
|ARa»l] are all different; for D two eigenvalues out of the three are 
identical, e.g., 


Gz = 3, B, = Bs; (3,21) 


* Inasmuch as Penrose's notation has achieved wide currency in works on Petrov’s classification 
thanks to the clarity of diagram (3.20), we will henceforth treat the field typesI, D, 0, II, 
N and III as if they were simply the corresponding Petrov types, without calling attention to 
departures from Petrov’s original terminology /65/. 
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for O the three eigenvalues are identical and, due to (3.15) are real: 
0, =Ag=A=—-—x, 8, =8, = 8, —0 (3,22) 


(for x =0 type O contains only flat space-time). 

In the second column ofdiagram (3.20), II and N are subtypes of 
*7,; Il is the ''nondegenerate second type,'' for which the eigenvalues 
o, and o, are different, and WN is the ''degenerate second type," 
for which o, and o, are identical and, from (3.17), real:* 


O=O,=—5x, B= B,=0, (3,23) 


The only eigenvalue, always real, in both blocks of the matrix 
Rall for *7, is given by —fx. 
3, Classification of gravitational fields of general form 

In view of the fact that the classification of Einstein spaces *T7; 
rests entirely on the algebraic properties of the curvature tensor, 


it is expedient to use Weyl's conformal curvature tensor (see 
/58/, p.115) to classify general gravitational fields (Ry, gas)! 


Capue = Radue — —-(uaRov — Barltoy + Bevan — Soytter) + 
a ~— R (SuaBin — Brahe) (3.24) 
It possesses all the algebraic properties of the Riemann tensor: 
Cappy = — Cpayy = —- Capp = Cyvap,  Cafayyy = 0. (3,25) 
It is easily shown that 
Cay = Capwe* = 0, (3,26) 
i.e., algebraically, the Weyl tensor of an arbitrary space VY, 


* We note that Bel /68/ employs a somewhat different notation: types I, , II, IIl and NW are 
spoken of respectively as cas 1, cas 2a, cas 2b, cas 3a and cas 3b. 
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behaves like a Riemann tensor in empty Space, and it is 
identical with the latter in the case (2.2). Determining the charac- 
teristic of the A-matrix ||Cap — Agap|| of the Weyl tensor in the 
bi-vector space R, and repeating the arguments of the preceding 
section, we conclude that there exist three and only three types of 
V, of general form, corresponding to the characteristics (3.13). 
One can show (/65/, § 20) that the matrix ||Call in the canonical 
orthogonal frame assumes the same form (3,11), (3.14) — (3.18) 
as ||Ry||, it now being necessary to set x =0 throughout. 

Penrose's diagram (3.20) remains unaltered for the matrix || Call 
as well. Type O obviously represents conformally flat spaces V,j, 
for which Cyagw = 0 always, while type NW is described by the 
matrices (3.16) for 


4, Petrov's classification and isotropic vector fields 


Debever shows /66/ that the Riemann space V,of signature -2 
admits of the canonical form of the matrix ||C,)||in the bi-vector 
space R, with respect to at least one and not more than four iso- 
tropic field vectors /* +0 satisfying the equations 


lCayotslel@l" = 0. (3.27) 


The formulation of Debever's theorem in the form of equations 
(3.27) is due to Sachs /110/. The importance of the researches of 
Debever and Sachs in what follows impels us to dwell on their 
results in somewhat greater detail. 

Type I according to Petrov is characterized by the fact that all 
four vectors liy WN 1, 2, 3, 4) are different; for type D they are 
identical in pairs (two independent vectors), for type II there are 
three independent vectors (two out of four are identical), for type 
type III there are two independent vectors (three out of four are 
identical) and, lastly, type N is characterized by the fact that all 
four vectors are identical, i.e., they determine the same 
direction.* To see this it is enough to verify equations (3.27) in the 
canonical orthogonal frame, specifying suitable values of the 


* The vectors /* in these equations are determined only up to collinearity (multiplication by an 
arbitrary scalar), Therefore a difference in the vectors I/jy, is interpreted as a difference in the 
isotropic directions specified by these vectors. 


40 GRAVITATION AL WAVES IN EINSTEIN'S THEORY 


components I*. Thus for type N one should take lt =8+8{. We 
will refer to the vectors ly, which satisfy equations (3.27), as 
Debever vectors. 

This brings us to a new invariant characterization of algebraic 
types of gravitational fields of general form. Designating the 
number of coincident (collinear) vectors [jy by a number in square 
brackets, one can systematize the foregoing considerations as 
follows: 


Petrov type I D Ii WN Ill 


‘Debever-Sachs symbol [1441] (22) [244] (4) (34) (3.28) 


The above system can be used in the algebraic classification of 
gravitational fields according to Petrov as an alternative to 
diagram (3.20). In this formulation of Petrov's classification the 
types of gravitational field are distinguished by the mutual orienta- 
tion of the Debever vectors in actual physical space-time. 

The type of mutual orientation of the Debever vectors is deter- 
mined by the specific form of equations (3.27); the most general 
form of these equations characterizes the ''most general" case of 
vector orientation, i.e., typeI. For the other types equations 
(3.27) transform into more stringent relations. The complete list 
for all the types in system (3.28) is as follows: 


Petrov type Equations for Debever vectors 
N or [4] Comal = 0 (3.29) 
III or [31] Capy[slthyy = 0 (3.30) 
II, D or [211], [22] Capypstayl?l* = O (3.31) 
I or [1111] FG phen (3,32) 


It is easily shown that a vector Jy, satisfying any one of 
equations (3.29) — (3.32) will automatically satisfy all the following 
equations as well. The assignment of a gravitational field toa 
given type is therefore governed by two factors: 1) the Debever 
vector Uy satisfies the given equation of the series (3.29) — (3,32), 
and 2) this vector satisfies none of the preceding equations in 
the series, 

So far we have been considering general gravitational fields, the 
classification of which is characterized by the algebraic structure 
of the Weyl tensor C.sys. Inempty space-time the Weyl tensor is 
identicalwiththe Riemanntensor Rjg.3; accordingly, Rags willoccur 
inthe classification (3.28) —(3.32) in place of Cygys. 
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Henceforth we will refer to type I fields, corresponding to the 
maximally general form of mutual orientation of the Debever 
vectors, as algebraically general, and to fields belonging 
to the other types —D, II, N and III —as algebraically 
Special. The meaning of these terms is evident from the 
Debever-Sachs approach, 


Chapter 4 


PIRANI'S CRITERION 


1. Isotropic electromagnetic field — null field 


The first attempt at a generally covariant geometric definition of 
the concept of gravitational waves in empty space based on Petrov's 
classification is due to Pirani (1957 /262/; see also his works 
/69 —71/). Pirani's definition is based on two postulates in con- 
formity with the concept of a wave in terms of Hadamard's 
discontinuity: 

1) the state of free gravitational waves is completely charac- 
terized by the Riemann tensor; 

2) the front of a gravitational wave manifests itself as a dis- 
continuity in the Riemann tensor on an isotropic three-dimensional 
hypersurface. 

The first postulate merely means that in Pirani's approach the 
Riemann tensor Rg, acts as field function (intensity) of the gravi- 
tational field, whereas the metric tensor g.,, often associated 
with the gravitational field potential, does not play a primary role 
in the analysis of gravitational waves, From the physical stand- 
point the second postulate, according to which the gravitational 
wave front lies on the characteristic hypersurface of Einstein's 
equations, means that gravitational waves in empty Space propagate 
with the fundamental velocity. 

Pirani's definition presupposes an analogy between electro- 
magnetic and gravitational waves. As a third postulate — not expli- 
citly formulated — it is assumed that in gravitation as in electro- 
dvnamics wave fields can only be isotropic. In this connection 
Pirani gives an extension of the concept of isotropic electromag- 
netic fields to the case of gravitational fields. 

We know that the energy-momentum tensor t,, of the electro- 
magnetic field (2,30) has four pairwise equal eigenvalues: k, k, —k, 
—k. Moreover, . 
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Ale = D* + 9, (4,1) 
where 
O=F,FY, P=PiFe, (4.2) 
and *Fy, is a tensor dual to the Maxwell tensor: 
"Puy = > Nvpal™*. (4,3) 


The electromagnetic field is termed isotropic (null field) if 
the eigenvalues of the tensor ty, are zero: k= 0, i.e, D=W= 0, 
From the physical standpoint an isotropic electromagnetic field 
corresponds to wave propagation of electromagnetic energy withthe 
fundamental velocity, as here the observer is unable to track the 
field. By definition an observer will be able to track the electro- 
magnetic field if for him the flux of the Poynting vector 


P,= (8: _ Vav®) Tpev? (4 .4) 


equals zero across all possible two-dimensional surfaces (here yx 
is the timelike velocity 4-vector of the observer). But in order 
for the flux of P* to be zero at all times, it is necessary that the 
vector P* itself be zero, According to (4.4), this means that 


Taal? == (TooV"V) Va (4,5) 


From condition (4.5) it follows that if the observer tracks the 
field, his velocity 4-vector will be the eigenvector of the matrix 
It ||. However, as is known /62/, the eigenvectors of the tensor 
ty of an isotropic electromagnetic field can be isotropic (and 
moreover identical) or spacelike, but they cannot be timelike. 
Hence it is impossible to make P* vanish in the case of the iso- 
tropic field, and in order for an observer to track the field his 
velocity 4-vector must degenerate into an isotropic one, i.e., 
the observer must travel at the speed of light. 

Thus one can give the following definition of an isotropic electro- 
magnetic field: an electromagnetic field is termed isotropic if the 
matrix |tywl| has no timelike eigenvectors. 
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2, Principal Riemannian vectors. 
Tracking the gravitational field 


Pirani's definition is based on an extension of the notion of an 
isotropic field (in its latest formulation for the electromagnetic 
field) to the case of gravitational fields. Such a generalization is 
difficult to accomplish directly, however, as Einstein's theory of 
gravitation lacks a true energy-momentum tensor for the gravita- 
tional field. To remove this difficulty Pirani defines the concept 
of tracking of the gravitational field in a way which differs in prin- 
ciple from electromagnetism, introducing the so-called principal 
Riemannian vectors. 

The principal vectors of the Riemann tensor 
Ropys are the directions of intersection of pairs of two-dimensional 
surfaces defined by the eigenvectors of the tensor R,, (a, b =1, 2, 

., 6) in the bi-vector space R, (and therefore by the bi-vectors 
in the physical space V,). The number and orientation of the prin- 
cipal Riemannian vectors can be determined if one knows the 
eigendirections of the tensor R,, in Rs. Thus for type I of 
Penrose's diagram (3.20) there exists a single (timelike) principal 
vector; for type D there are two isotropic and one timelike prin- 
cipal vector; and for types II, III and N there exists a single 
(isotropic) principal vector. 

We now introduce the following definition, in accordance with 
Pirani: an observer will be tracking the gravitational field if his 
velocity 4-vector coincides with a timelike principal Riemannian 
vector. 

Obviously, anobserver canonly tracka gravitational field the 
principal vector of which is timelike, i.e., afield of type lor D. 
As to type II, II] and N fields, an observer tracking such gravita- 
tional fields would have to have an isotropic velocity 4-vector, i.e., 
he would have to travel at the speed of light. 

We can now take the following definition of an isotropic gravita- 
tional field: a gravitational field is isotropic if the curvature 
tensor A, ,,s has no timelike principal Riemannian vectors, 

It follows from the foregoing that isotropic gravitational fields 
comprise fields belonging to types II, III] and N, and only to these. 
We now come to the formulation of the criterion for the existence 
of gravitational waves according to Pirani. 

Pirani's criterion. In a given region of empty 
Space-time V, free gravitational waves are present 
if within this region the Riemann tensor belongs 
to type II, WN or III of diagram (3.20); otherwise 
gravitational waves are absent, 
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3. Example. Wyman-Trollope gravitational wave fields 


Thus the definition of the state of a free gravitational wave field 
based on Pirani's criterion (and on other criteria which will be 
discussed later) is closely linked to the determination of the type 
to which a given gravitational field belongs within Petrov's classi- 
fication. It is therefore worthwhile to examine those known solu- 
tions of Einstein's equations in empty space which belong to 
types II, MN and III. 

Of these solutions a considerable number (those of Peres, 
Takeno, Petrov, Robinson and Trautman, Kundt et al.) will be 
discussed in subsequent chapters when we analyze other suggested 
criteria for gravitational waves (in a vacuum, or ina medium 
filled with electromagnetic radiation). To illustrate Pirani's 
criterion, let us look at a class of solutions of Einstein's equations, 
recently obtained by Wyman and Trollope /72, 73/, which does not 
coincide with the above-mentioned solutions or generalizes them 
considerably, 

To the Wyman-Trollope class of solutions corresponds the metric 


2a 1 B 1 
141 0 0 0 

8a = Q et O | (4.6) 
y¥ 00 e 


where 


= E46, pone + no 


and §., 60,0 @ and +t are functions of the coordinates 2°, z', x; 
p is a function harmonic in stand 2°: 


Bar = Biss = 0, 


and y is a function harmonically conjugated to BB. The metric 
(4.6) was obtained by using the expansion of ggg in the orthogonal 
frame along four real vectors (see, for instance, /74/), one of 
which — the isotropic vector /*— is assumed to be harmonic, 

Wyman and Trollope succeeded in integrating the field equation 
in empty space (2.2) with the metric (4.6) for a number of special 
cases, They single out three special cases: 


A, & +7° = 0, 
BLO 4? +0, ao. = 0, 
C. 8 4n° 0, ta +433 = 0. 


46 GRAVITATIONAL WAVES IN EINSTEIN'S THEORY 


It can be shown* that cases A and C correspond to type III fields 
in Penrose's diagram, and case B to type II fields, If the functions 
£ and yn are suchthat B and y are independent of z’, the given 
metric (in a vacuum) belongs to type N (degenerate type 2, accord- 
ing to Petrov's classification); inthis case /* coincides with the 
Killing vector, which defines the group of motions of this space- 
time along the coordinate lines 2z?. Geometrically the trajectories 
of this vector are interpreted as bicharacteristics of Einstein's 
equations, 


* The corresponding calculations have been carried out by L.B.Grigor’eva. She demonstrates that 
the metric discussed in Trollope’s work /73/ and corresponding to the case of-a nonharmonic 
propagation vector 1% is identical with the well-known Robinson-Trautman metric (see 
Chapter 9). 


Chapter 5 


BEL'S CRITERIA 


1, Superenergy tensor 


The criterion for the existence of gravitational waves proposed . 
by Bel /56, 68, 76 — 80/ (see also Debever /81/), like Pirani's 
criterion, is based on an analogy with the theory of electromagnetic 
waves. But unlike Pirani's criterion, Bel's first criterion, to 
which this section is devoted, rests on the definition of the ''energy 
tensor'' (more precisely the ''superenergy tensor") of the gravita- 
tional field. By analogy with the energy-momentum tensor of the 
electromagnetic field (2.30), such a ''superenergy tensor" should 
obviously be quadratic in the Riemann tensor Ryg,s. 

Let the given V4, then, be an empty space-time — so that Ein- 
stein's equations have the form (2,2) — and let the indices af and yé 
be the antisymmetric pairs of the tensor R,,,3. We introduce two 
new tensors, duals of the tensor Rag: 


- 1 oe * 1 ao 
Rapys —) Zz Napea Rey) Rapys = > Thr 890 ReeaB (5. 1) 
It can be shown that in Einstein spaces (3.7) 
* def » 


*Raprs = Rapys = Rapes (5.2) 


(see Appendix I, Theorem 1). 
We now define Bel's superenergy tensor /76/ as the 
fourth-rank tensor 


pee _ + er Ree + Rere Ree ). (5.3) 
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It is easily seen that in itS properties the Superenergy tensor 
(5.3) in empty space exhibits a close analogy to the electromagnetic 
field tensor 1t,,. First ofall, it is completely symmetric, thanks 
to the fact that the tensor Baus in empty V,is symmetric with 
respect to the index pairs af and yd. Secondly, as in the case of 
the tensor t,, it gives zero upon contraction with the metric 
tensor: go,/**"=0. Thirdly, it satisfies the covariant equation 
of continuity, which is analogous to the equation of continuity for 
the tensor ty in a space-time without sources /80/: 


Tvauv:a = O. (5.4) 
Fourthly, the analogy between the tensors 7s and t, also 


emerges when one compares the eigenvalues of 1 and the invari- 
ants of 7”, Tensor tj, satisfies the relation /81/ 


tite = Meer, (5.5) 


where & is the eigenvalue of the tensor tag, having the form (4.1). 
One can show /81/ that the superenergy tensor satisfies the anal- 
ogous relation 


: 1 
Peale = sj Rear, (5.6) 
where 
Kk = (Rapys R°**)2 + (Raars Rw es (5, 7) 
The close algebraic analogy between Bel's tensor and the 
energy-momentum tensor of the electromagnetic field can be 
employed to define the ''energy and momentum density'' of the gra- 
vitational field. Let a timelike unit 4-vector y« be Specified at 


every point M of space-time. We construct the scalar 


Ww=7T" yy ql glyUby. (5.8) 


It is easily shown /80/ that this scalar can be expressed as 
W (ut) = + (Bap8" + Hap Hl™, (5.9) 


where the symmetric tensors 6,, and #,,, introduced by Matte 
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/82/, are defined as 


Bar = Raprylur, (5.10) 

Har = — RapryuBur. (5. 1 1) 

These tensors are obviously oriented in space, in the sense that 
the timelike vector yw is their eigenvector (corresponding to zero 


eigenvalue). Hence it is easily shown that the square of each of 
these tensors is nonnegative: 


Bup8** >0, FHagdt*® > 0, (5.12) 


the equality holding only when the tensor is zero, 

It is moreover possible to prove the following theorem /80/: if 
for a certain vector u* the quantities @,, and #,, are simultaneous- 
ly zero, it follows in a necessary and sufficient manner that 
Rapys = 9, i.e., the space-time is flat.* Indeed, let the following 
equalities hold: 


RaprayuPur — 0, Rapryubur. = 0. 


According to the lemma proved in Appendix I, the second of these 
is equivalent to the relations 


(Uy Raga + UrxrRapyy + UpRapyr) ue = 0. 
Multiplying them by wu’ and exploiting the condition u? =1, we obtain 
Rapryw = 0. (5.13) 
From this follows (for Einstein spaces) the equality 
Riprywe® = Rupr? = 0, 
which, in turn, is equivalent to the relation 


Uy Roary + UaR ory + UpRrary = 0. (5, 14) 


* We note that in /82/ Matte writes down the equations of the gravitational field in an empty 
space in the language of the quantities (5.10) and (5.11), so that in the first approximation they 
are analogous to Maxwell's equations, in which it is the quantities @ and # that play the 
role of electric and magnetic field strengths. From Matte's point of view this analogy is suffi- 
cient demonstration of the reality of gravitational waves. 
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Multiplying (5.14) by ut and using (5.13), we obtain Regy=0, which 
proves the theorem. 


2. Energy and momentum of the gravitational field 


While discussing the concept of gravitational waves from the 
standpoint of energy transport by them, Synge /83/ formulates 
two necessary conditions which must be satisfied by the function 
F (u*) expressing the energy density of the gravitational field: 

1) F(u%) >0, and 2) if F (ut) =0, then Resy =0, i.e., the field 
energy vanishes only in the absence of the field itself. It is easy 
to show that the scalar W (u*) satisfies both these conditions: the 
first is satisfied due to relations (5.9) and (5.12), the second by 
virtue of the theorem proved above. 

Thus the scalar W (u*) can be taken to be the definition of the 
"energy density’ of the gravitational field. * As the vector u* is 
timelike, one can choose a local coordinate system in which 
uz= 65. In this system, obviously, W= T°, just as the concept 
of nongravitational energy is linked to the component 7” of the 
energy-momentum tensor of ''matter'' in Einstein's equations. 

Consider now the vector /76, 84/ 


P* = (85 —u,ut) TP uguauy, (5.15) 


which can be termed, by analogy with Poynting's electromagnetic 
vector (4.4), the Poynting vector (or Ssuperenergy 
flux density) of the gravitational field. 

In the same coordinate system, obviously, 


P=0, p=, 


As Bel shows /76/, in the linear approximation the following rela- 
tion holds at every point: 


Wy =—Pi, (5.16) 


whence from Gauss's theorem 


* With allowance, of course, for the density properties required for correct application of such 
concepts (the energy density according to Bel, say, would be — in the physical sense — the 


quantity Y—g W). 
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a, | Wav =— | Pin, ax, (5.17) 
Vv b>} 


where > is the two-dimensional boundary of the given three- 
dimensional volume V and n' is the unit vector of the outward 
normal to 2. Formula (5.17) means* that the gravitational super- 
energy flux across an element dz is proportional to Ptn;. 
Consequently, in order for the superenergy flux across any 
surface 2 surrounding a given point to be zero, it is necessary 
and sufficient that P*(uz) =0. Thus we come to the formulation of 
the criterion for the existence of gravitational waves according to 
Bel. ; 

Bel's first criterion. The presence of a superenergy 
flux is necessary and sufficient for the exis- 
tence of free gravitational waves. Conse- 
quently, gravitational waves are present in the 
neighborhood of an arbitrary point of empty 
Space-time V, if, for any timelike unit vector uw 
at this point, Psuy+~0 If P*u%)=0, there are no 
gravitational waves in the neighborhood of this 
point. 


3. Equivalence of the Pirani and Bel criteria 


We will now prove the rigorous equivalence of the first Bel 
criterion and Pirani's criterion (this proof was obtained by Bel 
/80/). Thanks to their covariant character, it is sufficient to do 
so in the local coordinate system used above, where u* = §. In this 
system, obviously, 


Pi = Ri, R™ — —S Riso Rojor, P=, 
ik 


or, differently expressed, 


‘ 4 a 
P= = Cime™ F 


* Strictly speaking the partition (5.16) is noncovariant and therefore all further statements 
should be qualified by the phrase "in the given coordinate system.” 
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where e™i is the three-dimensional Levi-Civita symbol. The three- 
dimensional antisymmetric tensor Cj, has the form 


Cim = Dy (8 jx Hem — 8 mnFe3)- 


& 


As Bel shows, ina bi-vector space with basis bi-vectors 
constructed on the natural unit vectors of the coordinate 
system chosen by us (the enumeration of the indices in the bi- 
vector space corresponds to the choice (3.9)), 


Gay =e Ato Sa =a Alo 8) = &) A es 
A Xn n~ 
Ea) =&a)/\%ay 8) = 4a ACa, Fo =em/A las 


the matrix of the curvature tensor can be expressed as 


Pe = a 


where @ and # are the matrices ||&;,|| and ||#;,|| of the spatial 
components of the tensors (5.10) and (5.11). Obviously, if 

P+(u*) = 0, then P' =0 and therefore C;,, =0, i.e., the matrices 
|$in|] and ||#:,|| commute. But in order for two three-dimensional 
matrices to commute, it is necessary and sufficient for them to 
reduce simultaneously to the diagonal form in a certain basis. 
Thus it follows necessarily from the condition Pt =0 and from 
formulas (3.14) that the corresponding field Rggys is of Petrov type 1. 
Conversely, if P*(ut) #0, i.e., the corresponding gravitational 
field describes gravitational waves according to Bel's first 
criterion, it will belong to type 2 or 3, i.e., it also satisfies 
Pirani's criterion. This proves our claim. 


4, Invariants of the curvature tensor in an empty space 


Bel's second criterion, formulated by him in /68/, is based — 
like Pirani's — on an extension of the concept of isotropic fields, 
known from electromagnetism, to the case of gravitational fields. 
Unlike Pirani, however, Bel arrived at the definition of the iso- 
tropic gravitational field by generalizing the concept of an isotropic 
field as one whose invariants vanish, rather than by generalizing 
the concept of tracking the field. 
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Whereas the number of functionally independent scalars which 
can be formed from Maxwell's tensor is two, from the Riemann 
tensor it is possible to form 14 of which only four are nonzero in 
empty space /85, 86/: 


4 eee 4 <i Bie 
A = 3 RR, B — so RRM, 
; ; (5.18) 
oe oe ee ee eo # oo 
CHa RR Ri, = D= ag RRS Rp 


these Bel calls fundamental scalars. Then, defining the 
isotropic gravitational field by the condition 


A=B=C=D=0O, (5.19) 


we obtain a new criterion for the existence of gravitational waves. 

Bel's second criterion. The field of free gravitational 
waves is identical with the isotropic gravitational 
fields defined by condition (5.19) for Regys#O0. An 
empty space-time with Riemann tensor Ragy+ 0 
will describe free gravitational waves if all four 
fundamental scalars (5.18) vanish. Otherwise free 
gravitational waves are absent. 

Let us now determine which Petrov types of gravitational fields 
satisfy Bel's second criterion. Writing down the conditions (5.19) 
in bi-vector space and using the canonical form of the curvature 
tensor matrix (3.11), (3.14) — (3.18), we see that conditions (5.19) 
are satisfied by three out of the six types in Penrose's diagram: 
O,N and III. Dropping type 0, which is trivial ( Ags, =0), we see 
that Bel's second criterion is satisfied by all 
gravitational fields belonging to types N and III 
of Penrose's diagram and only by these. 


5. Debever vectors and Bel's second criterion 


Types O, N and IIIl in Penrose's diagram in empty space (for 
which the eigenvalues of the tensor R&,, equal zero) are known as 
degenerate types of gravitational fields. In Bel's classification 
/68/ they constitute the single type he calls cas 3. 

As Debever shows /81/, in order for anempty V, to belong to 
Bel's type 3, it is necessary and sufficient that it admit of the 
existence of a vector /* satisfying the equations 
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Rapysltl* = 0, 


Rapygl2l* — 0. (5.20) 


The vector field I* (Debever vector) is unitary and isotropic. 
Thus Bel's second criterion admits of the following equivalent 
formulations /75/: 

Bel's second criterion (new formulation). An empty space- 
time with Riemann tensor Rgegys~O describes free 
gravitational waves if and only if it admits of the 
existence of an isotropic vector field /* satisfy- 
ing equations (5.20). 


Chapter 6 


LICHNEROWICZ'S CRITERION 


1. Bilinear degenerate form of Maxwell's tensor 


The criterion for the existence of gravitational waves proposed 
by Lichnerowicz /87—90/ (for a comprehensive exposition see /62/) 
is also based on an analogy with the method for determining the 
state of electromagnetic radiation. The latter rests on the Solution 
of Cauchy's problem for the Einstein-Maxwell equations in the 
space-time YV,. Let us briefly consider the salient points of 
Lichnerowicz's approach to this problem. 

Let the electromagnetic field tensor Fz, belong to class C* 
(piecewise C?), According to Hadamard's formula (2.13), the dis- 
continuities in the first derivatives of F,, on the characteristic 
hypersurface @ (z*) are given by 


[Foa,yl = faply (, = ,y)» (6.1) 
where f,, are the discontinuity coefficients for the tensor Fy,. It 


then follows from the first group of Maxwell's equations that the 
discontinuity coefficients satisfy the equations 


Lafpy + Upfoa + lyfap = 0, (6.2) 
and from the second group of Maxwell's equations that 
i = 0. (6.3) 
Our main assumption will be that the discontinuities in the 
electromagnetic field tensor along the wave front are proportional 
to the field itself, i.e., that f,, ~Fa,g. Then the tensor Fap 
obviously satisfies equations of the type (6.2) and (6.3) 


ltaF py= 9, [Fag 0, (6.4) 
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whence it follows necessarily that the vector i* is isotropic: 

1,/*= 0, A second characteristic property of the vector l* which 
satisfies equations (6.4) is that the lines of the vector field 

form a congruence of isotropic geodesics /260, 261/. 

In conformity with Lichnerowicz, we will refer to the bilinear 
antisymmetric form satisfying equations (6.4) asa singular (or 
degenerate) form ofsecondorder.* One canthen formulate 
the following theorem (Lichnerowicz /62/, $7): the coefficients 
Fz, of a singular bilinear formina V, of signature -2 are given by 


Fop= labg — pba; (6.5) 


where b, is a certain vector orthogonal to J, (blt=0). An obvious 
consequence of this is that if the components of Maxwell's tensor 
Fa, are coefficients of a singular bilinear form, they define aniso- 
trepic electromagnetic field. Conversely, an isotropic electro- 
magnetic field F,, is the coefficient of a singular bilinear form. 
Earlier it was shown (Ch. 4) that the electromagnetic radiation 
field can be defined as the isotropic field corresponding to the 
vanishing of the invariants of the Maxwell tensor. Now we see that 
an isotropic electromagnetic field can be defined in turn as the field 
of the Maxwell tensor whose components F,, are the coefficients of 
asingularbilinearform. This establishes the following definition 
(equivalent to that given in Chapter 4): Maxwell's tensor 
Fy# 0 describes electromagnetic radiation if 
there exists a (necessarily isotropic) vector field 
IX satisfying equations (6.4). 


2. Degenerate double form of the Riemann tensor 
Let all functions in the metric g,,(x*) belong to class C'(piecewise 
c?), According to Hadamard's formula (2.14), the discontinuities 
in the second derivatives of g., on the characteristic hypersurface 
(2,15) are given by 
[8ap,c0] = Gaplols. (6.6) 
Inserting (6.6) in the expression for the Riemann tensor 


4 
Rapay ~ 3 (Za2,6n om Sppar — Sayer — 8ar,ay) fe Kasrys 


* 2-forme singuliére, according to Lichnerowicz's original terminology /62/. 
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where Kgg» does not contain second derivatives of the metric and 
hence has no discontinuity on S, it is easy to obtain the expression 
for the discontinuities of Regry: 


[Reprul =~ (darbely + apylals — daylela — aprlaly). (6.7) 
It obviously follows from here that 
La{ Rayay} + ls (Ryan) + ly [Roary] = 9. (6.8) 
If we assume further that there is no discontinuity at S in the 
energy-momentum tensor Tf,, appearing on the right in the Einstein 


equations and that therefore, byvirtue of (1.1), [Ag,g] =0, then, 
according to Lichnerowicz (/62/, $20), 


Ie [Rasa] = 0. (6.9) 
Let us make the following fundamental assumption: let the dis- 


continuities of the gravitational field described by the tensor Rag, 
on the wave front S be proportional to the field components: 


[ Rapys| ~ Raays- 


Then the curvature tensor should satisfy equations of the form 


L, Rays + la Rares + lpRravs = 0, (6.10) 
Pian 0; (6.11) 


from which it follows that the vector /* is isotropic for Rggy3+~ 0. 


Indeed, contracting equations (6.10) with J and taking (6.11) into 
account, we obtain 


(771) Repro = 0, 
which proves our claim. 


We now introduce the following definition /88/: we will say that 
every tensor 


Hapry (= — Hogan = — Hapy, = Arya) 


defines a degenerate double form* if there exists a 


* Double-forme singuliére in Lichnerowicz's original terminology. 
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vector /* such that Hg. 0 satisfies the equations 
ln apis = 0, td «OT = 0, (6, 1 2) 


Lichnerowicz draws the following three consequences from 
equations (6.12): 

1) i is isotropic (the proof is analogous to the preceding one); 

2) i= 0, i.e., the vector field i* defines a congruence of iso- 
tropic geodesics; 

3) the tensor Hag ys can be expressed in the form 


Haar = bartaly oe Doylalr =< baylalr aoe Dartaly, (6. 1 3) 
where 


bar = Dras banl® = 0. 


Then, contracting expression (6.13) with g8*, we find that the 
tensor Hg = Hagyg** has the form 


Hay =tlal, (v= 0%). (6.14) 


On the basis of the foregoing results, one can formulate the 
Lichnerowicz criterion for the existence of gravitational waves as 
follows. 

Lichnerowicz criterion, The space-time VY, describes 
the state of the total gravitational radiation if its 
Riemann tensor Repys(*0) forms the coefficients of 
a Singular double form, i.e., there exists an (iso- 
tropic) vector #+#0, satisfying equations (6.10) — 
(6.11). If no such vector exists, there is no gravi- 
tational radiation. 

The Ricci tensor for a gravitational field satisfying the 
Lichnerowicz criterion should, from (6.14), have the form 


Rasp = thal g. (6.15) 
Conversely, conditions (6.11) and the isotropy of the vector i 


follow from conditions (6.15) and (6.10). Correspondingly, (6.10) 
follows from (6.15), (6.11) and the isotropy condition for i /62/. 
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3. The Lichnerowicz criterion and Petrov's classification 


It is evident from formula (6.15) that for 1 =0 the Lichnerowicz 
criterion defines gravitational waves in empty space-time ("purely 
wavelike" gravitational field). Then (see Appendix I) equations 
(6.10) and (6.11) become equivalent, the one following uniquely 
from the other: 


ly Raps =0. I? Rass = 0. (6 wl 6) 


This means that a purely wavelike gravitational field is uniquely 
defined by any one of these systems of equations. 

As is known (/62/, §21; /68/, $6), an empty V,belongs to 
type N of Penrose's diagram if and only if there exists a vector 
field i* satisfying one of the systems (6.16). Thus all spaces 
V,of type N define purely wavelike gravitational 
fields; conversely, all purely wavelike fields 
are of type WN. 

For t #0 in (6.15) we obtain the case of so-called total 
gravitational radiation,* In this case classification of 
gravitational fields is based on the Weyl tensor (3.24), which was 
mentioned in Section 3 of Chapter 3. In order to determine what 
Situation in this classification is occupied by fields of total gravi- 
tational radiation according to Lichnerowicz, we will prove the 
following auxiliary statement: if ina V,of signature —2 there is a 
vector i* satisfying the Lichnerowicz conditions (6.10) — (6.11), 
then this vector also satisfies the equations 


lpCapivs = 9, 
l*Capap = 0, 


(6.17) 


i.e., the Weyl tensor in the given V,defines a singular double form, 
Indeed, inserting expression (6.15) in the Weyl tensor (3.24) and 
making use of (6.10) and (6.11), as well as of the isotropy of the 
vector /*, after a direct check we are convinced of the validity of 
equations (6.17). 

One can show (Lichnerowicz /62/, §21) that if a certain tensor 
Hap having all the algebraic properties of the curvature tensor in 
empty space, 


Hapys ant — F pays = — F aps, = Hysag, 


(6.18) 
F atavs) =0, A aaysg?® = 0, 


“Radiation totale in Lichnerowicz's terminology. 


XN 
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definesa singular double form, then in the bi-vector space A,the 
matrix || H,,| of this tensor reduces to the canonical form charac- 
teristic of type Nin Penrose's diagram (degenerate type 2 in the 
Petrov classification). 

Let a space V,satisfy the Lichnerowicz criterion. Then from 
the foregoing proof it admits of the existence of a vector /* satisfy- 
ing equations (6.17). Taking the algebraic properties of the 
Weyl tensor (3.25) — (3.26) into account and making use of 
Lichnerowicz's results, we arrive at the conclusion that the tensor 
Capys belongs to the degenerate type 2 according to Petrov's 
classification. Thus we obtain the following theorem: a space- 
time V, satisfying the Lichnerowicz criterion of 
total gravitational radiation belongs to the 
degenerate type 2 as regards the algebraic 
structure of the Bayley tensor /91/. 


4, Conformal representation of gravitational wave fields 


To illustrate the applicability of the Lichnerowicz criterion, we 
will consider a study undertaken by Konopleva /92, 93/ of the con- 
formal representation of purely wavelike gravitational fields, Let 
V,be the space-time defined by the conditions 


Rit = 0, 
Ri = —4VCinyR”, 


where A is the gravitational constant, and C’», isthe Weyl confor- 
mal curvature tensor.* As the Einstein spaces (3.7) satisfy these 
conditions trivially, we will assume that the given Y,is a non- 
Einstein space conformal to a certain Einstein space V4. The con- 
ditions for feasibility of such a representation can be reduced to 


iL footy 
Rien = 2eP Rye sy, 


* These relations were obtained by the method of compensating fields of the Yang-Mills type in 
the theory of gravitation with the Lagrangian 


1% » 
L=R+ ta Rw k™ . 


In the Einstein theory they are to be interpreted not as field equations but rather ‘as conditions 
for distinguishing a certain class of gravitational fields (including, in particular, all Einstein 
spaces). 
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where Ray is the Riemann tensor for the space V, and 0o,=40. 
Hence VY, satisfies the condition 


Ras, = 0, 


i.e., V,is an Einstein space of the Petrov degenerate type 2, 
satisfying the Lichnerowicz criterion for pure gravitational radia- 
tion; o, is an isotropic vector describing the propagation of the 
gravitational wave front. 

Just as in electrodynamics the equations 


Fryv;z] = 0 
enable one to express Maxwell's tensor F,, in terms of the vector 


potential A, and its first derivatives, so the initial equations 


ee 
Ruyteng = Reiss = 0 


make it possible to express the Ricci tensor in terms of the vector 
o, andits first derivatives: 


Ryr =>— 20a: + 20,5 + Bun. 


In view of the fact that the vector o, is a gradient and making 
use, by analogy with electrodynamics, of the Lorentz-type condition 
ol, = 2x/3, we obtain the wave equation for the vector oy: 


8*75u; ag = 0. 


Thus in the Einstein space V, the vector o, satisfies the 
Lichnerowicz condition for a purely wavelike gravitational field, 
while in the original V, conformal to it, it describes gravitational 
waves in the same sense in which the vector potential of the electro- 
magnetic field satisfies the covariant wave equation. 


Chapter 7 


THE ZEL'MANOV CRITERION 


1. Generalized wave operator 


The criterion for the existence of gravitational waves formulated 
in /94/* on the basis of a general idea of Zel'manov presupposes the 
use of the following covariant generalization of the 
wave operator: 


D=— gVV.. (7.1) 


The generally covariant wave equation for an arbitrary tensor field 
Qx:} will then be 


DOK = AX, (7.2) 


where AX} is a tensor containing no derivatives higher than the 
first of Qx!. An equation of the form of (7.2) was used by Tolman 
/63/ to describe electromagnetic waves in Riemannian space-time. 

It is evident that a homogeneous equation of the type of (7.2) is 
trivial for the metric tensor g,,, and that in the case of Einstein 
Spaces it is trivial for the Ricci tensor Rg, too. Hence the idea of 
writing such an equation for the Riemann tensor Ag,gys, which is 
what A.L. Zel'manov suggested. Note, however, that equation 
(7.2) for tensor Rggys, 


DRaays = Kapyts (7.3) 


* The earliest work on gravitational waves based on the wave operator (7.1) dates back to 1962 
(see Ivanenko's foreword to the Russian translation of the book by Weber /95/). Subsequently 
Roy and Radhakrishna independently suggested /96/ defining gravitational waves by means of 
the application of the operator (7.1) to the curvature tensor. 
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is also trivial (i.e., becomes an identity) in the case of symmetric 
Spaces, for which the Riemann tensor is covariantly constant, 
Rapys;c = O(in this case, of course, one must have Kggys =0). 
Equation (7.3), moreover, can also become an identity for certain 
special choices of the tensor K,,,3+-0. Thus if we take 


Kags = 4 Re Rages — Rigo! Rag. — 2 (Rete Ripe + Read. Rese’), 


then equation (7.3) is satisfied identically for all Einstein spaces 
(3.7) /97/ (see also Appendix I, Theorem 3), 

Thus if one approaches the concept of gravitational waves on the 
basis of equation (7.3), it is necessary to demand, first, that the 
Space-time V, not be symmetric, and second, that it be specified 
by a choice of the tensor Kg g,3. The Zel'manov criterion is based 
on the assumption* that Kaa, = 0. 

Zel'manov criterion, The space-time V, describes 
gravitational waves if and only if its Riemann 
tensor Aggy: 1) is not covariantly constant, i.e. 
Rapysic 0, 2) satisfies the generally covariant 
wave equation: 


DRapys = 0. (7.4) 


2. Characteristics of the generalized wave equation 


In order to establish the physical base of the Zel'manov 
criterion — which might seem at first glance to be somewhat formal 
and artificial — it is important to investigate the characteristics of 
the tensor equations (7.4) as a system of differential equations of 
second order in which the unknowns, strictly speaking, are the 
components of the curvature tensor Aya, . 

As Savel'eva shows /99/, the left-hand side of equations (7.4) 
can be reduced by identity transformations to the form 


a” Rapys; po — B°°O_ Raprs + Lapys (T°) oF Quays (Te) + Qapys (Re). (ts 5) 


where the following notation has been introduced: 


* In /98/, equation (7.3) is discussed under the assumption that K agus F 0. 
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Lagys (f°) as (rey <7 Tey) Rages + (T° a ie Ropys + 
5 (Is —— hs) Rasys + (Cs nae ITs) Rapys I Ragys, ols 


Quays (oa) = (Tyet'y — 8 Tay Dev) Rages + 
+ (Ty era— BVT ov) Roars + (Mies — a Thal ev) Roos + 
+ (Civ g's —— eT halo) Ragys a ge” (TivRapas, ¥ + 
+ TiaRaprs v-+ MisRapys, ») + Pfs apres, ») — 
~~ g (Cy Rapes, ¥ ar Pa Reprs, B te PvaRaoys, Us + IsRagyo, p)s 


Oasys (Ri) are [Ry Rapes + Ri Rapys + RaRaors + RsRapyol : 


Here, asearlier, [T° = — g*®Ij,. We see that in the identity 
(7.5) the term Lag, is expressed algebraically in terms of the 
quantities [f° and their derivatives and that it therefore vanishes for 
f= 0, i.e., in the harmonic coordinate system. The term Q, ys 
is expressed algebraically in terms of the Christoffel symbols I, 
and therefore canbe made to vanish at any specified point in space-time, 
Theterm Qagys is identically zero in empty space-time (Rap = 9). 

Let a harmonic coordinate system be introduced in the region of 
empty space-time under consideration. Then the terms L,,,s and 
Qasys in expressions (7.5) vanish identically. Let us further assume 
that the chosen coordinate system is locally geodesic at a certain 
specified point M, Thenat M, firstly, the term Q,,,5 vanishes; 
secondly, the matrix | g°*] assumes the canonical form for the 
signature (+, —, —, —). Thus equation (7.4), considered as a 
system of second-order partial differential equations for the com- 
ponents of the curvature tensor A,.,,s, iS hyperbolic. 

Examining the Cauchy problem with given initial data Rag, and 
Rapyso ON a hypersurface S of the form (2.16) (in a locally chosen 
coordinate system), we reduce the system of equations (7.4) to the 
form 


8° Ra pyd,00 ata 0, 


where the dots denote terms not containing derivatives Raggy. and 
therefore completely determined by the initial data, Hence the charac- 
teristic equation for the system (7.4) inthe given coordinate system is 
the equation g° =0. Writingit in general coordinates (as was done for 
Einstein's equations in empty space in Chapter 2), we arrive at the 
conclusion that the characteristic hypersurface of the 
system of wave equations (7.4) is identical with the 
characteristic hypersurface of Hinstein's equation 
(2.15), where the function @ satisfies the eikonal 
equation (2,22). 
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Consequently, the generally covariant wave equations (7.4) in 
empty space describe the propagation of the discontinuities of the 
second derivatives of the Riemann tensor R,,,3 along the bicharac- 
teristics (isotropic geodesics) of Einstein's equations. 

Furthermore, from the notation (7.5) it can be concluded that 
in empty space, in a certain coordinate system (namely the 
harmonic system) the generally covariant system (7.4) becomes 
anordinary system of wave equations for each component of the 
Riemann tensor in the neighborhood of any point M: 


go” Rapys, pe — 0. (7. 6) 


This links the Zel'manov criterion directly to the ordinary meaning 
of "local waves of curvature” in the neighborhood of a point M. 


3. Zel'manov criterion and Petrov's classification 


In order to determine the Petrov type of the gravitational fields 
satisfying the Zel'manov criterion, we will confine ourselves to the 
case of the Einstein spaces *7;: 


Rag = 8 o9° (Ta) 
In spaces *7,, aS is known, the following identity holds /97/: 
8°’ Rapys; oa + Rape Ryo. +2 (Resa-Row — Reog-Racy + *%Rapy) = 9 (7.8) 


(see Appendix I, Theorem 3). From (7.8) follows an obvious result: 
in order for the Riemann tensor of the space *7; to satisfy equa- 
tions (7.4), it is necessary and sufficient that it satisfy the 
conditions 


Rigs. Riso. + 2 (RisoeRapr — Resa’ Racy. -+ %Rapys) = 0. (7.9) 


Writing conditions (7,9) in the bi-vector space R, in a canonical 
nonholonomic orthogonal frame and making use of the canonical 
form of the matrix of the curvature tensor (3.11), (3.14) — (3.18), 
we arrive at the following conclusions /100, 101/. For the spaces 
*7, conditions (7.9) lead to the system of equations 


a, (a, — as) — By (6, — Bs) = 0, 
By (a, — a3) + a (B, — Bs) = 0 
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and four more equations, obtained from the above by cyclic 
permutations of the indices 1, 2, 3. These equations are the 
integrability conditions, written down in R,, of the equations 


Rapys,c = 0, (7.10) 


which define symmetric spaces (/65/, p.399). The new system of 
equations obtained by covariant differentiation of the integrability 
conditions is identically satisfied by virtue of the initial equations 
(7.10). Consequently, any space *7, defined by conditions (7.9) is 
a symmetric Space. 

For the spaces *f7, relations (7.9) can be written in the canonical 
orthogonal frame as a system of equations /101/ which are 
consistent only under conditions 


a, = a, = 0, 8, = §, = 0, x= 0. (7.11) 


As we saw, these conditions define an empty space 7, of the 
degenerate second type (type N on Penrose's diagram). Conse- 
quently, the spaces *7, defined by relations (7.9) can only be 7, 
Spaces (x =0) of type N. 

Conversely, it we write down relations (7.9) in the canonical 
orthogonal frame in the bi-vector space R, and recall conditions 
(7.11), we see that relations (7.9) are satisfied identically, i.e., 
any Einstein space Tf, of type N will satisfy conditions (7.9). 

Finally, writing down conditions (7.9) in the canonical orthogonal 
frame for the spaces *7,;, we see that for any x they lead toa 
contradiction. This means that the spaces *7, cannot satisfy 
conditions (7.9). 

As Petrov shows /57/, there exist only two symmetric spaces 
*T,; they belong to the degenerate type 2 of spaces Tf,(x =0) and in 
the special coordinate system they are described by the metrics 


ds* = 2dz°dz! — sh? 2°dx?* — sin? (x° + k) dx?” 
ds? = 2dzx°dz! + ch® 2%dx?* + cos* (2° + k) dz* (7.12) 
(& = const). 


Thus, Einstein spaces satisfying the Zel'manovyv 
criterion can only be empty spaces Y, (x= 0) of 
Petrov type WN. Conversely, any empty space-time of type N 
— except for the two symmetric spaces (7.12) — will satisfy the 
Zel'manov criterion. 
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4, Relationship between the Zel'manov and Lichnerowicz 
criteria. Examples 


The relationship between the Zel'manov and Lichnerowicz 
criteria in the case of empty V, is easily determined on the basis 
of Petrov's algebraic classification, We already know that free 
gravitational waves in Lichnerowicz's sense of the term corres- 
pond to type WN fields on Penrose's diagram. Consequently, any 
empty V, satisfying the Zel'manovyv criterion will 
also satisfy the Lichnerowicz criterion, Con- 
versely, any empty Wsatisfying the Lichnerowicz 
criterion — except for the two spaces (7.12) — also 
satisfies the Zel'manovyv criterion, 

For general V, spaces (Rag #%ga,), the general relationship 
between the Zel'manov criterion and Petrov's classification — and 
also between it and the Lichnerowicz criterion — still remains 
unclear, However, one can already state that there exist gravita- 
tional fields of the form Ras #88 that satisfy the Zel'manov 
criterion. A number of such solutions of Einstein's equations with 
a nonzero right-hand side 


Rag — > RE gg = — Maps (7,13) 


interpretable as the energy-momentum tensor of an isotropic 
electromagnetic field, are given in /102, 103/. 

The isotropy condition for the electromagnetic field Fy, 
formulated in Chapter 4 imposes a rigorous restriction on the 
energy-momentum tensor tg,, and therefore also — due to the field 
equations (7.13) — on the Ricci tensor, In order for the metric gy,» 
to describe a gravitational field with the energy-momentum tensor 
of an isotropic electromagnetic field, it is necessary and sufficient 
that the corresponding Ricci tensor R,, satisfy the following set of 
conditions: the Rainich-Wheeler algebraic conditions /36/ 


R=0, Raph = + 88 (RR) = 0 (7.14) 


and the Nordtvedt-Pagels differential conditions /104/ 


nore (Roy. a Ruy = Roy;o Ry) = 0. (7. 1 5) 
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In turn, ina number of cases /103, 105/ the existence of the 
electromagnetic wave front leads to an additional symmetry pro- 
verty of the space-time V,, interpretable in terms of groups of 
motions of the given VY, (continuous groups of coordinate transfor- 
mations preserving the functional form of the metric; see, for 
example, /58, 106/). Thatis, the metrics of these V, admit of 
5-parametric groups of motions G, (4 <r< 6) which leave unaltered 
the isotropic three-dimensional hypersurfaces *V, that act as the 
wave front surfaces. If any point of a certain surface in V, can be 
shifted via transformations of the group G, to any other point of 
this surface, the latter is called a surface of transitivity of the 
group G,. Thus the isotropic surfaces of transitivity of the groups 
of motions acting in such gravitational fields either are identical 
with electromagnetic wave fronts of (for r<3) belong to to them. 

Let us look at a few metrics satisfying conditions (7,14) — (7.15), 
i.e., describing gravitational fields generated by electromagnetic 
radiation and accordingly admitting the group of motions &G,(r>3), 
which acts transitively on the isotropic hypersurface of the wave 
front /103/: 

1) the space-time described by the interval 


ds? = 2dx°dz!+ a (x°) dx?” + 28 (x°) dx’dx3+ y (2°) dx?” (7.16) 
admits of the group of motions G3, G; and G, acting transitively on 


the isotropic three-dimensional hypersurfaces *V, (here ay — B? > 0). 
2) The space-time described by the interval 


ds*= dx°* —dx?* — 225dz? (dx® + dz?) + a (x9+ x) (dx®* + dx3*), (7.17) 
admits of the group of motions G, acting transitively on the iso- 
tropic hypersurface *V;. The requirement that the density-pressure 
of electromagnetic radiation be positive — Rog < 0 — imposes the 


following condition on the function a for a Lorentz signature at the 
point: 


2, 
a<0, (2) —205%>-41 w=2?+2%) 
3) The space-time described by the interval 


ds? = du dv — 22°du dz®4 a(u) dx®*® + 28 (u) dz*dx® +. y(u)dx*, (7.18) 


in other words, the space time with the metric tensor 
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1 0 —zx* 0 
0 —1 —2z' 0 

—x2? —2 o« Bf’ (7,19) 
0 o 3 


admits of the group of motions G, acting transitively on the iso- 
tropic hypersurface *V;. The requirement that Ry< 0 fora 
Lorentz signature reduces to the following conditions upon the func- 
tions entering in (7.19): 


a<0, m>O, 
am|[ Fa + (ie) — gene —1]— (gr) <9 


where m= ay — B?. Obviously, these conditions are easily satisfied 
by choosing suitable functions a, § and y. 

A direct check easily shows that the metrics (7.16) — (7.19) 
satisfy the Zel'manov criterion. Moreover, as shown in /107/, the 
spaces VY, with these metrics admit of an isotropic field J* = 8 
satisfying the equations (6.10) — (6.11), i.e., they describe gravita- 
tional waves in Lichnerowicz's sense as well. 

However, there are solutions of the gravitational equation 
(7.13) with +t .g 0 which satisfy the Lichnerowicz criterion but not 
the Zel'manov criterion. Examples of this kind — obtained in /91, 
107, 108/ — will be discussed in Chapter 10. 


Chapter 8 


OTHER CRITERIA FOR GRAVITATIONAL FIELDS 


1. The Debever criterion 


Debever's approach to the question of the wave properties of 
gravitational fields is based on their relationship in empty space- 
time with the isotropic vector fields. 

As we know from the works of Debever /66, 81, 109/, Bel /80/, 
Sachs /110/ et al., in an empty V, one can associate with the 
Riemann tensor two isotropic two-dimensionalhypersurfaces, which 
are jointly spanned at every point by four isotropic vectors: [ty (N = 1, 
2,3,4). Making use of the canonical form of the matrix ofthe curvature 
tensor in bi-vector space | R,,], one can show that in every empty 
V, there exists at least one and at most four isotropic vector 
fields [iy + 0 satisfying the equations 


LpRay py [slot = 0 (8.1) 


(Debever's theorem /66/ in the formulation of Sachs /110/). 
Let us define, in accordance with Debever, the superenergy 
tensor 


4 
Vaprp == Bs P [ (eo Pru + By Pap + 8arPou + BauPar + 
96 
+ Bar Pay + Bay Pox) ja L Loar] ’ (8.2) 


where p is an arbitrary nonzero scalar and we have introduced the 
completely symmetric tensor 


1 


Linn = 5 gy eb Pn Ui (8.3) 


(L, M,N, 


in which P denotes summation over all permutations of the indices. 
Moreover, in (8.2) we find the quantities 


10 
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x oQee 12 
L= Ly, Dnyy. = Leis Py = L Dy a Baw 


and it is assumed that L +0. 

The Debever tensor (8.2) in empty space-time has the following 
properties /81/: 1) it is completely symmetric; 2) all its contrac- 
tions with the metric tensor vanish: 


EV apru = 0; (8.4) 
3) the tensor V,,,, is conserved: 
Vary: a = 0; (8.5) 


4) the tensor V,,, satisfies the algebraic identity 
9 y 
VapayV™ = + p?(1 — 30) &:, (8.6) 


where 


__ 9 (bas)? (Lia)? + (Las)? (Fas)? + (he) (laa)? 


=a (lasha + laslea 4 lala? , 


a 
Lonny = Bag lanliny. 


We can easily see that the properties of the Debever tensor (8.2) 
reveal its profound similarity to the Bel superenergy tensor (5.3), 
Going over to the canonical orthogonal frame and taking the 
concrete form of the isotropic vectors Jz, we can show that, upto 
a scalar factor, the two tensors are identical: 


Vapay S¢ Lapa. 


Hence the Debever tensor, like the Bel tensor, can serve as a 
basis for defining the ''superenergy flux" of the gravitational 
field. Thus the Debever criterion is completely 
equivalent to the first Bel criterion. 


2. Gravitational waves of the integrable type 
(Hély and Zund-Levine criteria) 


In the works of Hély /111 —117/, and also of Zund and Levine 
/118 — 120/, gravitational waves are defined according to . 
Lichnerowicz but with certain additional restrictions upon the 
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space-time metric. Carrying into effect the idea of Lichnerowicz 
/62/, Hély examines gravitational fields satisfying the criterion of 
"oravitational waves of the integrable type.''* 

By definition, the space-time V, describes gravita- 
tional waves of the integrable type if it admits of 
a vector field & which satisfies, in addition to 
equations (6.10) —6.11), and the isotropy condition 
for & which follows from them, the following 
gpradient condition as well: 


le = Gal, (8.7) 


where 1 is a scalar function of the coordinates. 
An example of free gravitational waves of the integrable type 
is the gravitational field described by the metric /113/ 


(00) 
Eup = Bag + 29Xlals, (8.8) 


(00) 
where gas is the flat space-time metric, 


= __ ap of, 
ly aa Al == a -+ Yar ’ (8 9) 
Ml,=0, I%.=0, ° 


~ and f, are five arbitrary functions of the argument 7 and q is 
a scalar. 

An example of a gravitational field corresponding to the case of 
total gravitational radiation of the integrable type is considered 
in /114, 115/, The corresponding metric is written as follows: 


(00) 
Bap = Bap oP LU + lgua, (8. 10) 


where the vectors u, and 1, (8.7) satisfy a certain system of 
differential equations which follow from equations (6.10) — (6.11). 
Zund and Levine /119/ introduced a definition of "total gravita- 
tional radiation of specialtype’: the space-time VY, corre- 
sponds to the case of total gravitational radiation 
of special type if it 1) is conformally flat and 
2) admits of a covariantly constant vector field iz 
satisfying the Lichnerowicz conditions (6,10) — (6.11), 


* “De type intégrable" in the Hély-Lichnerowicz terminology. 
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as well as the equations 
I, = (In t),¢5 (8.11) 


where «+ is a scalar defined by the system of equations (6.15) which 
follows from equations (6.10) — (6.11). 

The class of gravitational fields satisfying this definition is con- 
siderably narrower than the class of fields describing total gravita- 
tional radiation of the integrable type according to Lichnerowicz. 
Indeed it is easy to show that the metrics which describe total 
gravitational radiation of the special type should be representable 
as /119/ 

(00) 
Bag = 79? (h) Bags (8.12) 
V=2°+x!, a = const +0, 


where the function @g (A) satisfies the condition 


@ (A) (aA? + 2BA + y)-*, 


and a,B and y are numerical coefficients. The components of the 
vector Jz are of the form (8.11), the scalar t has the form 


t = 2 g-§ [2 (g’)® — gq”] + const, (8.13) 


and the prime designates differentiation with respect to A. 

As Zund and Levine point out, not all gravitational fields which are 
fields of total gravitational radiationin Lichnerowicz's sense and 
the sense of equations (6.10) — (6.11) satisfy the Zel'manov 
criterion and equation (7,4). However, every total gravitational 
radiation of the special type in the sense of Zund and Levine, as 
one can demonstrate /119/, automatically satisfies the Zel'manov 
criterion. This point has a bearing on the question of the relation- 
ship between the Zel'manoy and Lichnerowicz criteria in the case 
of nonempty space-time. 


3. The Maldybaeva criterion 


The criterion for gravitational waves proposed by Maldybaeva 
/121/, like the Zel'manov criterion, is based on a generally co- 
variant generalization of the d'Alembertian. Instead of the oper- 
ator (7.1), however, the gravitational waves are described with the 
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help of the de Rahm operator /42/ 
A = dd + 44, (8.14) 


where d is the operator of exterior differentiation (see /122/) and 
8 the operator of co-differentiation (see /42/), i.e., for taking the 
divergence of skew differential forms. 

Due to its structure, the operator (8.14) is applicable only to 
polylinear skew differential forms. Consider, for instance, its 
application to the bilinear skew form 


F = 4 Figdat \ do? (8.15) 


where Fy, is the electromagnetic field tensor. Maxwell's equations 
for the electromagnetic field in a region free of sources can be 
written in the form of two groups of equations /123/: 


aF=0, $F =O, (8.16) 


from which it follows that the 2-form F satisfies the generalized 
wave equation 


AF = 0. (8.17) 


Equation (8.17) was first obtained by Tolman /63/. We note that, 
unlike the operator (7.1), the operator (8.14) has a number of funda- 
mental topological /42/ and group /124/ properties. The idea of 
applying this operator to the study of gravitational wave fields was 
put forward by Lichnerowicz /43/. 

To describe gravitationalwaves in terms of the Riemann tensor 
Rapys» one constructs a 4X4 square matrix || Q..,||, the elements of 
which are the bilinear skew forms 


Qaa = 4 Raps de \ det (8.18) 


(curvature 2-form" in Lichnerowicz's terminology /125/). 
According to Maldybaeva's definition, an empty space- 
time with the Riemann tensor Raggy defines gravita- 

tional waves if the corresponding curvature 
2-form Q,, satisfies the generalized wave equation 


AQaa = 9, (8.19) 
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where the operator Ais given by formula (8.14) and the 2-form Qu 
by formula (8.18). 

Expressing the operator A interms of the covariant derivatives 
with the help of the de Rahm formula /42/, we write equations (8.19) 
in the tensor form 


D Rays + 2Rovse Rig = 0, (8.20) 


where the operator D is given by formula (7,1). 
Extending the arguments to the case of the Einstein spaces (3.7), 
we obtain the equations 


DRaays + 2Rovse Ris. — 2xR, 5 = 0. (8.21) 
8 Bry 


Making use of the identity (7.8) and writing down equations (8.21) 
in the canonical orthogonal frame, we can show that they are satis- 
fied if and only if the Riemann tensor belongs to Petrov's 
degenerate second type (type N of diagram (3.20) for x =0, see 
/126/). 

Thus the Maldybaeva criterion is satisfied by 
Einstein spaces of the degenerate second typef, 
and only by these. In Einstein spaces, therefore, the 
Zel'manov and Maldybaeva criteria demarcate the same class of 
gravitational fields; symmetrical spaces (7.12) of the type N, which 
satisfy the Maldybaeva criterion but not the Zel’manov criterion, 
constitute an exception. 


4, Misra-Singh criterion 


The criterion for gravitational waves given by Misra and Singh 
/127, 128/ is based on the concept of an isotropic gravitational 
field definable with the help of the Matte symmetrical tensors (5,10) 
and (5.11) discussed above. It is easily seen, by passing to the 
coordinate system in which J, = 62 that the rank of each of the 
matrices || &,,|| and | #.,| does not exceed 3, Further, the tensors 
84, and #2. are oriented in space: 


Bar ur = 0, 
Har u> = 0. (8,22) 
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Thus in a certain sense the tensors 6,, and #., are analogous to 
the three-dimensional vectors E‘ and Hi? of the electric and 
magnetic fields /82/. It is therefore possible to construct a defini- 
tion of isotropic gravitational fields based on the analogy with the 
isotropic electromagnetic field (Chapter 4). 

In this approach the isotropic gravitational field is character- 
ized, firstly, by identity of the eigenvalues of the matrices | @,, | 
and | #.,|| (analogy with the identity of the moduli of the electric and 
magnetic field vectors) and, secondly, by the fact that one of the 
eigendirections is common to both matrices (analogy with the 
mutual orthogonality of the three-dimensional vectors E‘' and H*). 
Introducing a third requirement — that the eigenvalues of the 
matrices |@,, and | #a.l{ equal zero — we arrive at the Misra- 
Singh criterion: an empty space-time VY, having the 
Riemann tensor Rup, describes gravitational 
waves if and only if the tensors (5.10) and (5.11) 
satisfy the relations 


Barb” = Hark”, (8.23) 
Carpe” = Hakpt™ = 0. (8.24) 


In their study of the properties of gravitational fields in empty 
space satisfying the above criterion, Misra and Singh prove two 
important theorems. The first states that a space-time satisfying 
this criterion admits of the existence of the isotropic vector field 
le for which the following conditions are fulfilled: 


Redhat 0: (8.25) 


But according to the known result of Debever /66/ and Bel /80/ 
(see Chapter 3, Section 4), the empty space-time V, defined by 
condition (8.25) can only belong to types N or III of Penrose's 
diagram (3.20). 

The second theorem states the following: in order for the space- 
time V, to describe gravitational waves in the Misra-Singh sense, 
it is necessary and sufficient that its Riemann tensor satisfy the 
equations 


Rapys ee Rives = 0; (8 26) 


By writing down these equations in the canonical orthogonal frame 
for the canonical form of the matrix of the Riemann tensor, we can 
again show that they are always satisfied for fields of types N and 
III of diagram (3.20) and only for these. 
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Thus in order for an empty space-time VY, to 
satisfy the Misra-Singh criterion, it is necessary 
and sufficient that it belong to the Petrov type N 
or III. The Misra-Singh criterion is therefore 
equivalent to the second Bel criterion. 

The Misra-Singh criterion was extended by its authors /127/ to 
the case of nonempty space-times YV, (fields of total gravitational 
radiation), To arrive at this extension it is sufficient to formulate 
the definition of the isotropic gravitational field in terms of the 
Weyl tensor Cypy (see Chapter 5, Section 3). Indeed, by analogy 
with the Matte tensors @,, and #,., we introduce the symmetric 
tensors 


~ 


~ * 
Bar = Carry ®u®, Har = Capa eiu®. 


Defining the isotropic gravitational field in the general case with 
the help of the tensors @,, and #,, in precisely the same way as 
we defined it for the empty gravitational field with the tensors &,, 
and #,,, i.e., by means of the equation 


Capys CO’ Cures = 0, 


we can formulate the following theorem: in order for the space- 
time V,to satisfy the definition of total gravitational radiation of 
Misra and Singh, it is necessary and sufficient that it belong to 
type N or Ill as regards the algebraic structure of the Weyl 
tensor Cup,s. 


Chapter 9 


PROPAGATION OF GRAVITATIONAL WAVES 


1. Gravitational geometric optics 


So far we have been chiefly concerned with generally covariant 
formulations of concepts and criteria fundamental to the study of 
gravitational wave fields. The question of the propagation of 
gravitational waves, i.e., the analysis of definitions of wave front, 
oropagation trajectories (rays) and so on, has essentially been 
disregarded, 

By an analysis of the Cauchy problem for the system of equations 
of the gravitational and electromagnetic fields in the Riemannian 
snace-time V, it was shown (Chapter 2) that the fundamental ideas 
of geometrical optics are common to the electromagnetic and the 
gravitational field. Indeed, as in the theory of electromagnetic 
waves in Maxwell's electrodynamics, the law of propagation of 
gravitational waves is determined by the eikonal equation (2.22). 

The solution of this fundamental equation — 
the scalar function @g9(z) known as the eikonal — 
defines the hypersurface of the gravitational 
wave front (2,15), as wellasthetrajectories of wave propagation, 
which form a congruence of lines of flow of the isotropic vector J+ 
(2,25). The vector J« will be termed the wave vector ofa 
gravitational wave. 

To different solutions of the eikonal equation there obviously 
correspond different wave front types, which in turn physically 
define distinct types of gravitational waves. Hence the problem 
arises of establishing a generally covariant classification of gravi- 
tational wave types based on the properties of the wave front. 

Although there exists as yet no universally accepted definitive 
solution of this problem, ina number of instances it has proved 
vossible to give a generally covariant description of gravitational 
wave types based on certain specific wave front properties definable 
in terms of the wave vector /*, Indeed, by specifying in a 
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space-time V,an isotropic vector field which guarantees the 
existence of a solution of the eikonal equation — i.e., which is a 
gradient field as in (2.25) — we can classify gravitational waves 
according to the properties of this vector field. Our task is made 
easier by the fact that the solution of the problem of the existence 
of isotropic fields (Debever vector fields in particular) in the 
space-time V, is in many instances helpful in establishing the 
algebraic type of the space-time V, according to the Petrov 
classification. 

Two types of gravitational waves, plane andspherical, 
have been determined in this way. The names themselves are 
determined by the geometric analogy with plane and spherical 
electromagnetic waves, based on the similarity between the laws of 
geometrical optics for gravitational and for the electromagnetic 
fields. We will be making extensive use of this analogy in the 
physical interpretation of corresponding types of gravitational waves. 


2. Spherical gravitational waves. Examples. 


According to the Robinson-Trautman definition /129 — 131/, the 
metric of a space-time V, describes a field of spherical gravita- 
tional waves if the given V4 admits of an isotropic vector field / 
satisfying the equations 


Lra:p) = 0, (9.1) 
Darl ® — 5 (Ifa)* = 0, (9.2) 
20: . (9.3) 


The first of these conditions means /57/ that the vector field /* 
is a gradient field, i.e., that it can be represented in the form 
(2.25), Then, due to the isotropy of [, the equation 9 (z*)=a, 
where a is a parametric constant, describes the family of 
characteristic manifolds of Einstein equations satisfying the eikonal 
equation (2.22), 

As we saw earlier (Chanter 2), the trajectories of the vector 
field I* satisfying equations (2.22) and (2.25) constitute a family of 
bicharacteristics of the Einstein equations and are therefore iso- 
tropic geodesics in the space-time Vy with the metric gag. Then 
from the physical point of view condition (9,2) means the absence of 
distortion of the form of the shadow thrown upon’a screen by an 
opaque object illuminated by light rays moving along the 
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trajectories of the vector /* /131 —134/. Condition (9.3), in turn, 
implies magnification or reduction of the shadow compared with the 
object itself (unlike the case of plane electromagnetic waves, for 
which I, = 0 always /135/). 

Indeed, let a given congruence of isotropic geodesics with the 
tangent vector /* define the trajectories of propagation of light. 
Consider a small opaque object illuminated by the light rays and 
casting a shadow ona screen placed orthogonally to the rays. It 
can be shown /110/ that all parts of the shadow reach the screen 
simultaneously. We now shift the object by parallel displacement 
along the rays to the position occupied by the screen and compare 
its dimensions with those of the shadow. If the screen is placed at 
a distance dr from the object, its shadow will be subjected to the 
rotation wdr, the shear or distortion |s {dr and dilation or 
contraction edr; the scalars e, wo and o are given by /71, 136/ 


4 


1 1 
e= ya OF = zlpa; pl *, =F la;ml?—e*. 


The intervretation of the "optical scalars" e,w and o is based 
on the following hydrodynamic analogy. Let wu* be the field of the 
timelike unit vector by means of which we specify the 4-velocity of 
an ideal fluid. Then the kinematics of an infinitesimal volume 
element of the fluid will be described by the coefficients of the 
following expansion of the covariant derivative of the 4-velocity: 


1 ~ 
Ug: 8 = Wap +- Sap + Fz Nas — Ags. 
Here 
Nab = — gap i Ugls, 


a, = Ugpu® is the so-called ''4-acceleration vector" (equal to zero 
for geodesic congruence), @= u’,,, and the tensors 


Oap = NaN; + Sas = NANA Uo; « 


are respectively the tensors of the angular rotation velocity and shear 
rate of the volume element in the geodesic coordinate system (Fermi 
system) fixed along the trajectories of vector u*. The scalars of 


rotation 5 and shear (distortion) ¢& of the trajectories .u* are 
defined in terms of the tensors ,3 and oy: 
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~ ~, 1 
o? = es Wapo**, 6 = > Gago™#. 


As Kundt shows /137/, the condition o=0 is not only necessary 
but also sufficient for the fulfillment of condition (9.1), i.e., 
specification of a congruence of isotropic geodesics without 
rotation is equivalent to specification of the wave front surface. 

Robinson and Trautman /130/ have obtained the class of exact 
solutions of the Einstein equations in empty space describing 
spherical gravitational waves and comprising — as Foster and 
Newman demonstrate /138/ —all algebraically special 
gravitational fields, i.e., fields of Petrov types D, II, NM and III. 
This class of solutions is represented by the metric 


ds? = 2dz° da! + (K —2Hx!— =) dz — 
— «YP {{dr* + Q3dx|* + {dx? + Q,3ax°]%}, 
M = M(x), P= P (2°, z*, 2), @= @ (2°, x, 2°), (9.4) 


where the functions H and K are given by the formulas 


H = P™P,) + P (P™'Q)23 — PQ (P~),23 
K =P (P 92 + P33) =, (P\2)? pai (P,3)?. 


In this coordinate system, the gravitational wave front is 
expressed by the equation x®8=const while the isotropic vector I 
orthogonal to it is given by l#= gf, so that the trajectories of the 
gravitational wave (gravitational rays) are identical with the 
coordinate lines z?, 

The field equations (2.2) with respect to the metric (9.4) reduce 
to two relations: 


Mb page . (9, 5) 
wa — Kp = 4P7? (My — 3HM). 

The former allows one to make @Q vanish by a coordinate trans- 
formation (conserving the form of the metric), so that in the new 
coordinate system the field equations reduce to the single (second) 
relation in (9.5). 

Studying the Robinson-Trautman metric (9.4) in a coordinate 
system in which Q@=0, Bartrum /135/ demonstrates that it may be 
extended to the case of nonempty space-time, and more specifically 
to the case where an electromagnetic field is present. It can be 
shown, further, that the electromagnetic field is isotropic, i.e., 
both invariants (4.2) of the Maxwell tensor F,, vanish. 
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For the Robinson-Trautman metric, only one of the components 
Ty, is nonzero: 
Too = E*/0?, 


where 


ME? = +P? (Ka — K,23) — 2M 9 + (QM/P) Po. 


Making use of the theorem which he himself had proved earlier 
/139/ regarding the differential properties of the wave vector I+, 
Bartrum obtained a particular solution describing a self-consistent 
system of spherical gravitational and electromagnetic waves pro- 
pagating along general trajectories — the lines of the vector field [+: 


ds* = (2 — A) dr? 4+ Br’ (dg? 4+ sin’g dd?) — 
—2(1 — A) drdt— A d#, (9.6) 


where 


= )"(2tet)" se Mo, (9.7) 


n is an integer and the spherical coordinates r,,6,¢ are related to 
the ''Cartesian'' coordinates, in which @ has dropped out, by the 
usual transformation. A remarkable point is that in empty space 
for m =m, the metric (9.6) — (9.7) transforms into the Schwarz- 
schild metric. 

The Kerr-Schild solution /140, 141/ in empty space-time may 
also serve as an example ofa gravitational field admitting of inter- 
pretation in terms of spherical gravitational waves: 


(00) 
Bw = Bu + 2Hlyl,, (9.8) 


Here J, is anisotropic vector field: g, ly =0, whence follows 


(00) 
also the condition g/l? =0, and conversely; the function A(z*) 
together with the vector /* satisfies the field equations 


Qyabpy + (H + 2H) bea: 9) — He" la; ola;c +s blala=0, (9.9) 
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where we have used the notation 


dq = (H + 2He)q — 2H oli? — H1;.°; , 
b=2HH—e" Hap, e=—- ls, H=lH a. 


From the field equation (9.9) it follows /141/ that the vector [+ 
is tangent to the congruence of geodesics of zero shear: 6 =0. 
Assuming furtherthat this congruence has zero rotation, i.e., 
demanding that the vector /* be normal (ifa;g; =0), we obtain an 
empty gravitational field with gravitational waves propagating 
along the congruence /+, 

In order for this gravitational field to correspond to spherical 
gravitational waves (i.e., in order for the condition ¢e+0 to be 
satisfied), it is necessary and sufficient that it belong to Petrov 
type II or D(Mas /142/). If this condition holds, the metric (9.8) 
becomes 


(00) 
Bw = By + 2elyl,. 


On the other hand, in the case e =0 the corresponding gravita- 
tional field belongs to type N or III. This follows from the result 
of Kerr and Schild /141/, who demonstrated that the following 
relation is fulfilled for the metric (9.8): 


Rapys (°° = Hply. 


Moreover, from condition e =0 it follows that H =0 and thus 
the metric (9.8) satisfies the Bel criterion (5.20). In this last case, 
as we will see below, the Kerr-Schild metric describes plane 
gravitational waves. 


3. Plane gravitational waves. The Kundt definition 


A generally covariant definition of plane gravitational waves 
was given in two different ways by Bondi-Pirani-Robinson /143/, 
Penrose /144/ and Kundt /137, 145—147/. 

According to the Kundt definition, the metric of a space- 
time V, will describe a field of plane (or 'plane- 
front’) gravitational waves if the given Y, 
admits of an isotropic vector field /- satisfying 
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equations (9.1) —(9.2) and the condition 
lig = 0. (9.10) 


As we already know, these requirements reduce to the 
following: that the trajectories of the vector /* are isotropic 
geodesics. Here the first requirement, expressed by equation (9.1) 
(condition of existence of a wave front), means that the shape of the 
shadow cast by an object on a screen orthogonal to the rays (to the 
trajectories of I*) is not deformed as it would have been had the 
object undergone a rotation with respect to its true position. The 
second requirement, expressed by equation (9.2), means that it is 
not deformed as it would have been had the object undergone a 
displacement. Lastly, the third requirement, i.e., condition (9.10), 
means that the shadow is neither magnified nor reduced by com- 
parison with the object itself (see /71, 131, 148/). 

Thus in his definition Kundt accounts for the main properties 
familiar to us from the study of plane electromagnetic waves in the 
Minkowski space-time, where the front of a plane wave travels 
parallel to itself along trajectories of propagation of light orthogonal 
to it, without any distortion, 

Chevreton /149/ has shown that a field of plane gravi- 
tational waves satisfying the Kundt definition in 
an empty space V, can only belong to Petrov 
type O (trivial case, Rasy =0), type N or type III. 
If it belongs to type N, the vector /* will be covariantly constant: 
It = 0. 

All empty spaces V, of type N admitting of a covariantly constant 
vector field J* are known /105, 150/, and they constitute a class of 
solutions of the Einstein equations which have been determined up to 
the integration of a certain system of differential equations. Plane- 
wave solutions of this type will be discussed in detail in the next 
chapter. 

Kundt also obtains another class of solutions of the field equa- 
tions in empty space /137/, corresponding to plane waves in the 
sense of the definition which he introduced. Included among solu- 
tions of this class are fields of Petrov type N and type III. 

The Kundt metric was obtained in an isotropic tetrad, i.e., an 
isotropic quasiorthogonal frame {t*, 7, l¢, m*}, of which the vectors 
Ix and m* are real and the vectors ¢# and ?* complex conjugates; 
further, lm, = fi, = 1: 


ds? ~ |dz + Bdu P+2dvdu+H du’, (9.11) 


where z=z-+iy, and B and H are real functions of the coordinates, 


9. PROPAGATION OF GRAVITATIONAL WAVES 85 
such that 


B,,=9, Bax + By =9, H=—vB, +A, Ay = 9, 


A xx + Ayyy + 2 BB t3 BiB ix — 2 Bout By By = 2. i212) 

Here the lines u are timelike and the scalar +t is given by 
equations (6.15). Ina vacuum t =0, whence it follows, in view of 
the isotropy of J,, that for a nonempty space-time the metric (9.11) 
corresponds to the case of electromagnetic radiation. The general 
trajectories of propagation of electromagnetic and gravitational 
waves (lines of flow of the isotropic vector field J«) coincide with 
the coordinate lines v; in other words, if the v-coordinate is 
chosen as affine parameter on the congruence of the isotropic 
geodesics, the latter are given by the equations 


== [*, (9.13) 


and the wave front is given by the hypersurface u=const. 

The Kundt metric (9.11) is a special case of the class of metrics 
obtained by him in /137/. This class, according to Kundt, 
describes plane-front gravitational waves corresponding to more 
general conditions in which the vector /« is not necessarily normal, 
i.e., is not bound by the requirement that J,,:,; =0. 

In the case t =0 (empty space-time), the Kundt metric belongs to 
Petrov type N if B =0. According to the Chevreton theorem considered 
earlier, for B =0 the Kundt metric describes a space-time which 
admits of a covariantly constant vector field i. Such spaces VY, 
will be discussed in detail in the next chapter. 


4, Plane gravitational waves. 
Bondi-Pirani-Robinson definition 


As Goldberg and Kerr have shown /151, 152/, the asymptotic 
behavior of gravitational fields set up by isolated sources is very 
similar to the behavior of plane electromagnetic waves in the 
Minkowski space-time. 

This result may be regarded as the motivation which led Bondi, 
Pirani and Robinson /143/ to attempt to give a rigorously group 
definition of the concept of plane gravitational waves in empty 
space as a metric field satisfying two postulates: 1) the field is the 
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same at any point of the wave front; 2) the metric tensor, like the 
vector potential of plane electromagnetic waves, has a definite 
symmetry group. The group of motions G,, which leaves the iso- 
tropic hypersurface in V, unaltered, is sucha group. The equation 
of this hypersurface (which describes the wave front) is given, in 

a certain coordinate system, by 


x! — x° = const. (9.14) 


Petrov shows /57/ that there exists only one space V, of 
signature 2 admitting of a group of motions G, which acts transi- 
tively on the isotropic three-dimensional hypersurface. In the 
coordinate system of (9.14) its interval is written as follows: 


ds? = — Adz!” — 2D dxtdr? — B dz?® — Cdx®*4+Cdzx™, (9.15) 


where A, B, C and Dare functions of the retarded argument 2! — x° 
that satisfy the differential equation 


M* — 4M’ (In MY — M’ (nO) — AB’ —(D'? = 0, 


in which M = AB — D’? >0 and the prime denotes differentiation with 
respect to z!— 2°, 

Starting from results obtained by Takeno /153/, Johari /154/ 
extended the Bondi-Pirani-Robinson definition to the case of non- 
empty spaces V, satisfying equations (7.13). Hély /112/ and 
Johari /154/ examined the concept of plane-polarized gravitational 
waves and defined the parameters of polarization. (Problems of 
this kind had been discussed earlier by Rosen /155/, as well as 
Boardman and Bergmann /156/.) In this connection an interesting 
work which stands alone is that of Avez /157/, in which a definition 
is given of monochromatic gravitational waves. We will now con- 
sider this work in some detail. 


39. Monochromatic gravitational waves. 
The Avez definition 


Pursuing the general electromagnetic analogy, Avez /157/ calls 
a gravitational wave monochromatic if the correspond- 
ing wave vector (2.25) is harmonic, i.e., satisfies 
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the equation 
(V =e"). =(V—8 8" Pp),a = 0. " (9,16) 


Inthis casethe complex function U = exp (ig) is also harmonic: 


(V— 8 8 U p).a = 0; (9.17) 


this fact is invariant under the substitution g—/f(9), where fis any 
function of the class c?. Asa result Avez was able to study mono- 
chromatic gravitational waves as the field of a simple periodic 
function U (z*). 

Let the wave vector [+ satisfy the condition (9.1), and let, 
further, the isotropic geodesic congruence corresponding to it have 
zero dilation: 


It = 0. (9.18) 


Then, obviously, it will satisfy equation (9.16) and be harmonic. 
The corresponding plane gravitational wave (in Kundt's sense) will 
be monochromatic. 

Avez studied in greatest depth the spaces V, which admit of 
interpretation in terms of monochromatic gravitational waves. He 
started from the Einstein equations (1.1), taking the energy- 
momentum tensor of the ideal fluid plus electromagnetic field on 
the right (for Ty,): 


4 : 4 
Top = X |(0 + P) Uatlp — Pap +z Fiyl op — FapFi?]——z 8s. (9.19) 


Here p, p and uz are the density, pressure and 4-velocity, 
respectively, and X and @® are constants representing the arbitrary 
nature of the choice of the energy-momentum tensor. 
The solution of the problem is given by the following theorem of 
Avez: gravitational fields definable by the energy-momentum 
tensor (9.19) admit of monochromatic gravitational waves only if 
1) p =p = 0(i.e., the only source is the electromagnetic field); 
2) the gravitational wave front 9 (2*) is defined by the condition 


lesa = Lg&p f Lp8as (9,20) 


where e,is a certain vector orthogonal to 1; 
3) the wave vector I+ is an eigenvector of the tensor Fag: 


Fagl® = ely. , (9.21) 
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Noting that under these conditions 


PRapys = (Lap — Lt pa) lL, + Lpbax — labpy 
Hap = &6:a — babs 


and computing the quantities 
Ragys al oe Real l 


in accordance with Debever, we see that in the case of an empty 
space-time (7a, = 0) they vanish, i.e., they satisfy equations 
(5.20). This means, according to the second Bel criterion, that 
the given gravitational field belongs to type N or Ill. Thus, an 
empty space-time describes monochromatic 
gravitational waves only if it belongs to Petrov 
type N or Ill. 

To obtain exact solutions of the Einstein equations describing 
monochromatic gravitational waves, it is convenient to employ a 
coordinate system in which the components of the metric tensor 
are expressed explicitly in terms of the harmonic function 


U (x*) = exp (ig). (9,22) 
Then the g,,, which can be treated as the potentials of the gravita- 
tional field, are completely analogous to the vector potential of 
a plane monochromatic electromagnetic wave /151/: 


A* = a’ expi (kr — ot). 


Let the solution be sought in the form, say, 


(00) 
Buy = Bw + 2a (x*), (9.23) 


where a,, is a constant symmetric tensor (analog of the amplitude 
of periodic oscillations). As ''amplitude tensor" a,, we choose the 
tensor 


Qyy = Ley + lye, , (9.24) 


where J, is a constant isotropic vector and e, is a constant 
vector orthogonal to it.* Further, let the eikonal g in the chosen 


* In a coordinate system in which the metric g,, is given by (9.23), the relations expressing the 
covariant constancy and orthogonality of the vectors assume the simple form of the ordinary 
relations of constancy and orthogonality in a Galilean metric. 
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coordinate system be given by 


(00) 
Q = Lian = Zaal x", (9.25) 


so that 
U (xt) = A, cos (1.2%) + Boy sin (lax*) (Ag, By = const). (9.26) 


By substitution in equations (2.2) it can now be shown /159/ that, 
under the given assumptions, the metric (9.23) — (9.26) describes 
an empty V, corresponding to plane monochromatic gravitational 
waves. The analogy between such a gravitational field and the field 
of plane monochromatic electromagnetic waves will become yet 
more striking if we express the metric (9.23) —(9.26) in the equi- 
valent form /159/: 


(00) 
Buy = Bw + 2dyy, (9.27) 
Guy = yy EXP [i1z"] + dyy exp [— lax’), (9,28) 
Aus, yy = 5 (uly + Iyey) (Ay =F 1By)- (9,29) 


The well-known Peres metric /160/ 


ds?= dx®* —dx!* —dz?? —dz8* — 2f (dx® + dz*) (9,30) 


also belongs to the class of metrics (9.23) if one sets 


B= z = f(z, x, u), Po=uw=—2+ 2’, 


433 = oq = 43 = — 1 
(the remaining a,, components are equal to zero). The field equa- 


tions (2.2) then reduce to a single condition of harmonicity of the 
function fin the arguments z! and 2?: 


far + fee = 0. (9,31) 


Chapter 10 


PLANE GRAVITATIONAL WAVES DEFINED BY 
AN ABSOLUTELY PARALLEL VECTOR FIELD 


1, Plane waves in empty space-time 


As we already know from Chapter 9, the trajectories of propaga- 
tion (rays) of a plane gravitational wave in the Kundt sense for 
type WN fields in empty space are defined by the absolutely parallel — 
or, in other words, covariantly constant — isotropic vector field [+. 

The converse can also be proved /105/: every empty space- 
time admitting of an absolutely parallel vector field i= belongs to 
the Petrov type N, where /* is isotropic and unique. This 
follows from the conditions of integrability of the equations J, = 0, 
given by (6.11), and from the Debever result (3.29). 

Thus an absolutely parallel vector field in empty space defines 
a congruence of isotropic geodesics which represent the trajectories 
of propagation of plane gravitational waves according to Kundt, 

As Eisenhart shows /161/ (see also Kruchkovich-Solodovnikov 
/162/), the metric tensor of a space V, admitting of the unique iso- 
tropic absolutely parallel vector field Jl can be written as 


(10.1) 


where A,B,D, FE, F, Q@ are functions of the x, z*, 25 coordinates; the 
vector it in this same coordinate system is given by /* = 6. 

For the metric (10,1) only six of the twenty significant compon- 
ents of the Riemann tensor are nonzero: Roos. Rseo2» Rozoe Ase32+ Roses 
Roos - Out of the ten field equations (2,2), six are identically satis- 
fied, while of the remaining four three become 


Rosse = Rgeso = Raeog = 0, (10.2) 
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where 


4 
Rees oT D »5 7 (A,33 >. B 9) —_ + {A\s (A.3B _ BD) + 
+ By(B A — AsD) + Bs [A,D — A(2D.— A) — 
— (2D; — B.) [AB — D(2D,,— A,3)]}, 
4 4 
Rys20 = > (Dyoz — Ayos + Q,23 — Fe) + Te (Do+ F2—Q,s) X 
x (AD — AB») + Ayo (AB — BeD) + (Bs + Do — Fie) X 
x [D(2D,— A.) — AB) -+ By |A,D + A(Ag—2D,)), 
4 
Rges0 = a (Boz — Dos — Q,33 — F128) + 


+ Fe (Do + Os — Fa) (AB — BaD) + Bo(ABa— AgD)+ 


+ (Dy + Fx — Q3) (D(B,g — 2D,3) — ABs] + 
+ Apo [B(2D,3 — B,2) — BDI}, 


g = det} gaa] = AB— D?< 0. 


The last field equation links together the remaining components of 
the Riemann tensor: 


A Rosso = BRoeao = 2D Roeso = 0, (10.3) 


where 


Rore0 nae Q 02 a" > (Ajoo + E29) + z {(Do “~_ Fy a Qs) x 
x [24D + A (Q,s as Fy, —s Do)} a (Ay)* B + (Ey oe, 2F 1) x 
X [A.D + A(Ay — 2D,,)] + (EZ, — 20,3) X 
X [B3D + B(B.,— 2D,)}}, 


4 4 
Rosso = F 3 — > (Broo + E33) + Te {((Dyo + Os — Fs) X 


X [2B oD + B(Q3 + Fie — Do) — 
— A(Byo)* + (2F 9 — Es) [ABs + D(By —2D,3)\ + 
+ (Ey, — 205) [BsD + B(B, — 2D,)]}, 


Rorso = + (Fro a Qos = Ens _ D,o0)— 
— Fg (Dp + Bs — Fa) (BAy — D(Da + Fp — Qs) + 
+ Bo lA (Do + Fy2— Qs) — DAo] + (20,0 — Es) X 
x (DB, ad BA) + (2F is Es) (DA,s =~ AB,s)}. 


Thus the general metric (10.1) which satisfies the four field 
equations (10.2) — (10.3) defines the class of exact solutions of the 
Einstein equations in empty space satisfying the Lichnerowicz and 
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Zel'manov wave criteria as well as the Kundt definition of plane 
gravitational waves. Wecan show that this class comprises a 
number of known solutions of the Einstein equations belonging to 
the degenerate Petrov type II, in particular the solutions of Takeno 
/153, 163/, Peres /160/ — (9.30), Petrov /57/ — (9.15), et al. 

Thus the Peres metric (9.30) can be reduced by a simple 
coordinate transformation (rotation of the axes in the (r°, z’) plane) 
/164/ to the form 


1+ f(x®, 22, 23) 


0 
0 (10.4) 
at 

0 1 


The metric (10.4) is obviously a special case of the metric (10.1) 
for 


A=B=—1, D=Q=F=4), 
E=1+4+/ (2, z*, 2°). 
Here the field equations (10.2) are satisfied identically and equation 
(10.3) reduces to the condition of harmonicity of the function f in 


the arguments z* and z?: 


hee i hss = 0. 


The Takeno metric /153/ is 


ds? = — A dxt? — 2D dz'dz* — B da*” — (P — S) dx’* — 
— 28 dz'dx® + (P + S) dz”, (10.5) 


where A,8B,D,P,S are functions of the argument zx? — ¢° and the 
quadratic form 


dl? = A da'* + 2D da'dx® + B dz*” 
is positive definite and is a solution of equations (2.2) if 
Mss — (M,5)*/2M = A,3B,3/(D,3)" = 0, M —— AB _ D?. 


It can be shown /165/ that, like the Petrov metric (9/15), the 
above metric is reducible by a coordinate transformation and trivial 
relabelling to the form 
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ds? = 2 dx°dz' + A (x°) da®* + 2D (x) dx*dz? + B (2°) dz3?, 


whence it follows that it is a special case of the metric (10.1) if we 
choose E=Q=F-= 0 andif the functions A, B, D are regarded as 
dependent only on 2°, 

Lastly, another Petrov solution /65/, represented by the metric 


0 1-2 0 
10 0 0 
Be =| 23 0 A(x) D(x |? (10.6) 


0 0 D(z) B(2°) 


is again a special case of (10.1). For this metric AB — D*?<0, 
A, B-#=0, and among the field equations (2.2) the only nontrivial 
equation is (10.3). 

It is natural to suppose that the arbitrariness inherent in the 
definition of the metric tensor (10.1) of the Einstein spaces admit- 
ting of an absolutely parallel vector field is somewhat fictitious, 
i.e., stems from the choice of coordinate system; there should 
exist, then, admissible coordinate transformations * 


ae a 

zi = ft y (2, x. ?), (10.7) 
x aay p (2, x #°), 

x* = tp (2, 29, 2), 


which conserve the form of the metric gp» (10.1) and of the vector 
field J (= 83), and which may enable us to simplify the metric 
(10.1) and in particular to remove certain of the functions that 
enter into it. Moreover, the solution of the field equation (2.2) is 
itself defined by the form (10.1) only up to the integration of four 
second-order partial differential equations. This also reduces the 
number of independent arbitrary functions in the metric (10.1) 
quite considerably. In this connection it would be appropriate to 
attempt to integrate at least a few of the equations (10.2) — (10.3). 

Both problems were recently solved by Kaigorodov and Pestov 
/166, 167/, who show that in the new coordinate system (#) the 
components #3, Zy3» B33 Can be assigned the values 


* The tern admissible will be applied to transformations whose Jacobian does not vanish. 
Correspondingly, the arbitrariness of the functions % @ and pin (10.7) is restricted only by the 
following differential condition: J = Q,,5 — P,2P,3 F 0. 
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for a suitable choice of the functions g,p and X in the transforma - 
tion (10.7) (these functions are defined as integrals of a certain 
system of equations of the Cauchy-Kowalevski type). 

Thus in the general case the metric tensor of the Einstein space 
admitting of an absolutely parallel vector field can be represented 
as 


E 
1 

83 = Q ’ (10.8) 
0 


where A, Z,@ are functions of the coordinates 2°, x’, 23 satisfying 
equations (10.2) — (10.3), and the equation (10.3) becomes 


Rosso <a Roe = 0. (10.9) 


Making use of the necessary and sufficient conditions obtained 
earlier for fields of the degenerate type 2 — which involve only 
first derivatives of the vector fields, with the result that there is 
further simplification of the field equations /74/ — Kaigorodov 
reduces the arbitrariness in the definition of the functions £, Qg and 
A toa single second-order equation (10.9), which remains unin- 
tegrated. 

The final result of these investigations can be formulated as 
follows: in order for an HKinstein space to admit of 
an isotropic absolutely parallel field i, it is 
necessary and sufficient that in a certain coordi- 
nate system the vector & be expressed as f= 686, 
and that the metric Vs assume the form (10.8), 
where the functions £,Q,A are given either by the 
relations 


A = —41, Q = 2ex°, Ein + Eis3 = 4€ (e = 0, 1) 
(first class of solutions) or by the relations 
A = — [x5 + a]*, B = 2b (23° + 2a2%), 


Eo + (x° +4)? Eg — rer Eo + (2? + a) Bg + 


a os) [8A (r8 + a)? — 4x3 (2? + a) (Ade),9 — (23 + 2)? a,o9 — 
— 4734 95 (x3 + a) + 425a,,a,94} = 0, 
dX = const, a = a (2°, 2*), b = b (2°) 


(second class of solutions). 
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It is easy to see that the above-mentioned wave solutions of 
Peres, Takeno and Petrov are special cases of the two classes of 
solutions singled out by Kaigorodov and Pestov, or that they reduce 
to the latter by admissible coordinate transformations. 


2. Absolutely parallel vector field in a nonempty space-time 


The class of solutions of form (10.1) admits of extension to the 
case of nonempty space-time /91/. 

Let the space-time V, have a metric of the form (10.1) satisfying 
(aside from the condition dg,,/dxz' = 0) only the three equations (10.2). 
If one computes the components of the Ricci tensor Rg one soon 
sees that nine of its ten components vanish identically and that the 
only nonzero component has the form 


Roo ane — (AReoso + BRo220 = 2DRopgo)- 


A direct check will show that the Rainich-Wheeler conditions 
(7.14) and the Nordtvedt-Pagels conditions (7.15), which define the 
energy-momentum tensor of the isotropic electromagnetic field, 
are satisfied for this particular metric. 

However, the solution of the gravitational equations thus obtained, ” 
since it satisfies the Lichnerowicz criterion for total gravitational 
radiation (Ch. 6), describes the propagation not just of electro- 
magnetic waves but also of gravitational waves. Indeed, as the 
space V, with the metric tensor (10.1) admits of congruence of the 
lines of flow of the covariantly constant vector I, it follows that 
this vector satisfies the second Lichnerowicz condition (6.11). 

Next introducing the vector field 1, = 82 and the expressions for the 
components of the Riemann tensor into the first system of 
Lichnerowicz conditions, (6.10), one can see that it is satisfied 
identically. 

Thus this class of solutions of the Einstein-Maxwell equations 
describes a self-consistent system of free electromagnetic and 
gravitational waves in space-time without matter /91/. It com- 
prises a number of known wave solutions of the Einstein-Maxwell 
equations; these include, to name a few, extensions of the solutions 
of Peres and Bondi-Pirani-Robinson to the case of nonempty space- 
time, known as "fields of type P and H" (according to Takeno's 
terminology /102, 168/). 
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Special cases of these solutions satisfying not only the 
Lichnerowicz but also the Zel'manov criterion were discussed in 
Chapter 7. However, it should be noted that in the general case 
these solutions of the Einstein- Maxwell equations describe a self- 
consistent system of electromagnetic and gravitational fields 
(nonempty space-time) in which the gravitational field satisfies the 
Lichnerowicz criterion but not, as a rule, the Zel'manov criterion. 
This last point is relevant to the question of the general relation- 
ship between these two criteria in the case of nonempty space-time. 

As the class of solutions we have described satisfies the 
Lichnerowicz conditions (6.10) — (6.11), it follows, from the 
general theorem proved in Section 3 of Chapter 6, that they belong 
to the degenerate type 2 gravitational fields of general form. 

As we saw in Chapter 2, the characteristic manifolds of the 
Einstein equations and Maxwell equations — like their bicharacter- 
istics — are identical. Consequently, the trajectories of propaga- 
tion of plane gravitational waves defined by an absolutely parallel 
vector field /#, are at the same time the trajectories of light rays. 
At the same time, it follows from the coincidence of the gravita- 
tional and electromagnetic wave fronts that electromagnetic waves 
in a space-time admitting of an absolutely parallel vector field are 
also plane according to Kundt. 

It can be shown /169, 170/ that the class of solutions described 
above always admits of a group of motions with a singular operator 
/106/, the trajectories of which are identical with those of wave 
propagation (with the rays). Moreover, if the given group of 
motions is intransitive /106/, then its (always isotropic) surface of 
transitivity will either belong to the wave front or coincide with it. 


Chapter 11 


THE ASYMPTOTIC PROPERTIES OF FIELDS 
OF GRAVITATIONAL RADIATION 


1. Gravitational radiation of axially symmetric 
isolated systems. Bondi information function 


In studying the asymptotic behavior of fields of gravitational 
radiation, it is convenient to use the Bondi-Sachs method based on 
expansion of the quantities characterizing the field in powers of 1/r, 
where ris a parameter playing the role of distance from an 
isolated system of sources. Assuming that such an expansion is 
possible, we choose a coordinate system such that, firstly, the 
form of the terms which dominate for large values of r is 
simplified to the maximum extent, and secondly, their wave 
character is brought out. 

Bondi /20, 171/ analyzed the problem of the emission of gravi- 
tational waves by an axially symmetric system of bodies using the 
following space-time metric as initial metric: 


ds? = go dx” + 2godr%dz'+ Qgy,dx°dx? — gy_dz* — B33 4x5 °. (11.1) 


The form of the metric (11.1) was determined from the following 
physical analogy with the model of radiation in a flat space-time. 
Let a ray of light emerge from a certain point 0, surrounded 
by a small sphere on which the angular polar coordinates @ and 6 
have been specified. We introduce the timelike coordinate u, 
analogous to the ''retarded time" of the ordinary theory, and define 
the coordinates u,98,@ of an event £ to be the values of the corre- 
sponding coordinates of the point at which the ray OF intersects 
the sphere. As the coordinates u, 6, @ are constant along the ray, 
the trajectory of a ray of light in space-time is a coordinate line of 
the fourth coordinate r. The affine parameter r defined in this 
manner along the trajectory of the ray will be henceforth inter- 
preted as the distance from the source Q. In flat space-time, 
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choosing u=t-—r to be the timelike coordinate, we can express 
the metric as 


ds? == du® + 2 du dr — r? (d6? + sin? 0 dq). (11.2) 


Foran asymptotically flat Riemann space-time, the 
immediate generalization of this metric is the metric (11.1). 
Moreover, by choosing a coordinate system of the type of (11.2), 
Bondi was able to remove the so-called "logarithmic term" of form 
(In r)/r from the expansion of the Riemann tensor, the presence of 
which he regarded as a flaw in this approach. For the problem of 
radiation of an axially symmetric system, Bondi specifies the 
following components for the metric (11.1): 


Zoo = 1B exp 26 — r*A? exp 2y, 
Bo1 = XP 2B, Bop = Ar? exp dy. (11.3) 
Bo = — Tr? exp 2y, g5, = —r® sin? 8 exp (— 2y). 


Here f, y, A and B are functions of r, 6, u and are independent of gq. 
In the coordinate system so chosen, the axially symmetric metric 
(11.1), (11.3), as we shall see, contains no logarithmic terms in 
the expansion in 1/r. 

Of the ten field equations (2.2) three are satisfied identically for 
the metric (11.3): 


Rog = Rig = Rog = O. 
Of the remaining seven, four equations 
Ry = Ry = Re = Rag = 0 (11.4) 
are independent, one equation, 
Ry = 0 


is satisfied owing to (11.4), andthe components Ry, and Ry, due 
to the Bianchi identities, satisfy the relations 


(r? Ro), = 0, (7? Roo), = 0, (11.5) 


i.e., do not contain terms of order higher than the second in 1/r. 

The system of four equations (11.4) in the four unknowns 6, y, 4, 
B (‘fundamental equations,'' in Bondi's terminology) belongs to the 
group of equations R,; =0 which, as we saw in Chapter 2, 
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determine the behavior of the solution in the neighborhood of the 
initial hypersurface § depending on the initial data upon it. The 
system of fundamental equations has the following remarkable 
property: it contains no second derivatives of the type 8,0, and 
contains only one derivative of the type 8jrxo:1, namely Yo. This 
means that specification of the function y(r, 6, u) on the initial 
hypersurface u =const, 


y (r, 8) = y (r, 8, 4) |u~const » 


defines not only the functions y,f, A and B themselves, it also 
defines the derivative Y,0, i.e., the behavior of the solution in the 
neighborhood of the hypersurface uw=const, apart from the arbitra- 
riness in the choice of the integration functions. As the equations 
only contain derivatives in r of the sought functions, it follows 

that the five integration functions will depend only on the coordinates 
uand 6. One of these should be set equal to zero so as to satisfy 
the requirement of regularity of the metric for large r; a second 
function is eliminated by a coordinate transformation which con- 
serves the form of the metric. Asa result the expansion of the 
metric in powers of 1/r, 


y=Cu, O/r+C, (u, 9)/r?4+..., 
B = —C*/4r5 + ..., 


4 
U = —— (Cp + 2Cotg 6) + (11.6) 
4-4 [F (u, 8) + 300, + 4C*ctg 9] +... 
V=r+G6(u, 6)+4+..., 
is determined by three integration functions — C,Fand G. Since 


these functions are related by the two supplementary relations (11.5) 


Gyo = > (Cas + 3C2 Cbg 8 — 2C)o — (Co)*, (11.7) 


+ 


i 3F' 9 = Go + 3CC oe + 4CC 9 ctg 8 + C C2, (1 1.8) 


only one of them can be regarded as independent, say C (u, 8). Thus 
the value of C (u, 8) andofits derivative C, w.r. tothe ''retardedtime" 
u on the hypersurface wu = a-=const completely determines the 
solution in the neighborhood of the hypersurface, i.e., in a certain 
interval (a<u<b). Thus every irreversible change in the system 
of sources is determined by the quantity 
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Co= 2 Clu, 9), (11.9) 


which contains all the information concerning the behavior 

of the system with variation of u. The function C, is 

called the information function* of the system /171/. The 
most important property of the information function is that it 
determined the fraction of mass of the system borne away by 
gravitational radiation. 

Let a system of sources be stationary at all times until the 
instant u =a, and let it return to the stationary state after a period 
of nonstationarity — at the instant u=b. The definition of the 
information function makes it possible to obtain the mass lost by 
the system during the time interval a<u<b. The gravitational 
field described by the following Weyl metric is the stationary field 
corresponding to the axially symmetric metric (11.3) (see Synge 


FAt2}): 
ds* = exp (2 )dt? — exp (— 2 p)[exp (20) (dp? + dz*) + p* dg]. (11.10) 


Expanding the function »p in the multipole moments of the system 
of sources and writing down the metric (11.10) in the Bondi 
coordinates u.r,6,@, we obtain the relation between the functions 
G and F and the multipole moments of an axially symmetric 
system: 


=2M, F=2Dsino. (11.11) 


Here M is the mass of the system of sources (''monopole moment") 
and D is its dipole moment. Returning to the general metric (11.3), 
we can define the mass of the source system as the mean of G (u, @) 
over the angle: 


M(u)=+\6u, 8) sin 0 d0. (11.12) 


0 


Next, integrating the equation (11.7) we obtain the expression for 
the decrease of mass with time in terms of the information function: 


dM 1 
My = Fe = —+\ (Co) sin 9 d0 <0. (11.13) 
0 


* A term in frequent use is “Bondi news function." 
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Formula (11.13) means that the mass of an isolated axially 
symmetric system of sources will be constant if and only if the 
information function of the system equals zero. Otherwise the 
mass decreases monotonically. 

In view of this we can formulate the following criterion of 
gravitational radiation for axially symmetric isolated systems: 

Bondi's criterion. The gravitational field of an 
isolated axially symmetric system of sources, 
defined by the metric (11.1), (11.3), is a field of 
gravitational radiationif the information function 
(11.9) of the system is nonzero, Otherwise, there 
is no gravitational radiation.* 


2. Newman-Penrose formalism 


By computing the components of the Riemann tensor Rgyg,,s for 
the Bondi metric and expanding them in powers of 1/r, it can be 
shown /171/ that they break down into three groups of equal com- 
ponents (up toa sign). In the first group, the expansion begins 
with a term of the third order, 2 (G +CC)r-* , in the second group 
with a term of the second order, (Co. + 2C,ctg6)r-*, and in the third 
group with a term of the first order, Cyor-'. This means that for 
Cov # O the gravitational field decreases in inverse proportion to r 
at large distances. . 

In an attempt to remove the element of arbitrariness in this 
result which stems from the choice of a dominant" coordinate 
system, Newmann and Penrose undertook to study the tetrad com- 
ponents of the Riemann tensor in place of its coordinate components. 

In order to do this we will use the Newman-Penrose formalism 
/174/ which, in addition, yields an invariant formulation of the 
"laws of conservation" of multipole moments, 


* We note that the vanishing of the information function does not imply stationarity of the system. 
Indeed, from relation (11.8) it follows that for €,o=0 and G20 the quantity F,9 is nonzero, 
which, owing to the second formula of (11.11), means that the dipole moment of the system 
increases linearly with time. How much the monopole and dipole terms contribute to the 
energy and momentum lost by an isolated system is a nontrivial question which has attracted 
special attention. It was recently reviewed by Papapetrou within the framework of the linear 
approximation of Einstein's theory of gravitation /173/. It was found that the effective 
monopole and dipole terms which appear formally reduce in actual practice to the quadrupole 
moment of the source system (quadrupole moments of the so-called "electric" and "magnetic" 
types). At the same time, the quadrupole moment of a system cannot be regarded as an 
effective moment reducible to higher moments. 
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The Newman-Penrose method is based on the assumption that 
there exist in the chosen region of the manifold V, four isotropic 
differentiable vector fields* (as smooth as the components of the 
metric g.), two of which are conveniently chosen to be conjugate 
complex. In Bondi coordinates (2° =u, zt} =r, 2? =6,23 = @) these 
vector fields can be specified as follows: 


I] = dy, nh = bo > dy, 


mi = (85 + siggy 8S)» = (8b — arog 88) 


(11.14) 


It has been shown by Newman and Penrose /174/ and — inde- 
pendently of them and in a different form — by Kaigorodov /74/ 
that the metric Sw can be represented as the combination 


Suv = Lytey +- Nyly = myn, — Thum, 
with 


L,l¢ = nn! = fiym = 0. 


The ten independent real components of the Riemann tensor in 
vacuum (fy, =0) can be characterized uniquely by five complex 
scalars ("tetrad components of the Riemann tensor" **): 


Dy = Ragysl*milm®, 
D, = Rapysl2n8l*m, 
@, = Rasys menbl*ms, (11.15) 
Ds = Reays M*nhlrn8, 


O, = Rays menbmn’, 


Gravitational radiation is said to be emergent from a system or 
incident upon the system depending on whether the initial data (the 
information function in our case) are given on an isotropic hyper- 
surface in the absolute future, or in the absolute past (in the 
linearized theory these solutions are called retarded or advanced 
potentials), Assuming that only emergent waves are present inthe radi- 
ation (i.e., disregarding the scattering of the incident external radiation 


* In the general case the existence of such vectors is not at all guaranteed. However, for gravita- 
tional fields induced by isolated distributions of sources such vectors always exist, given 
certain asymptotic (boundary) conditions. 
** Such components can be defined for any tensor having all the algebraic properties of the 
Riemann tensor in a vacuum, e.g., for the Weyl tensor. 
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by the system*), and also specifying a definite class L of functions 
@,(A=0, 1, 2, 3, 4) in formulas (11.15), we can show /175, 176/ 
that in an asymptotically flat space-time one of the scalars @,, 
say ®,, admits of expansion of the form 


L 


Oy = > r+) (D4, + O (r-G+L)), 


n=0 
where the expansion coefficients ®% are independent of r. 
Let this scalar be @®,. According to Newman and Penrose /177/, 
it is sufficient to take it in C’: 
OD, = Dor + Mors +. Mor? +. O(r*). (11.16) 
Then, making use of the Bianchi identities and of the identity 


oe 
Rony. 6 0, 


which follows from them in the case (2.2), we obtain the asymptotic 
expansions of the remaining @,: 

®, = Ofr+*+ O(r3), OD, = OF 3+ 0(r>), (11.17) 
0,= OF +0(r4), @,=O%740(r), 


where @°, 2, 08, @{ are independent of r. For the Bondi metric 
the coefficients @% depend only on 6 and u, and 


oO = C (u, 9),00- 
In the case of axial symmetry the coefficients % of the 
expansions (11.16) — (11.17) can be expressed in terms of the asso- 


ciated Legendre functions /182, 183/: 


oo = s a, (u) P2 (cos 8), D$ = >) b, (u) Ph (cos 9), 


anne n=] 

5 S .18 
©. = Sic, (u) P2 (cos 6), O° = >} e, (u) Ph(cos 0), (11.18) 
sie n=] 


co 


Of = >) h™P? (cos 8). 


n=2 


* The scattering of gravitational radiation incident upon a system of bodies has been discussed by 
Zerilli /178, 408/, Vishveshwara /179/, and Couch, Kinnersley and Torrence /180, 181/. 
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In this form the expansion coefficients are analogous to the 
retarded potentials of the wave equation of the special theory of 
relativity. 

Indeed, in flat space-time the axially symmetric solution of the 
wave equation in polar coordinates is given by 


pm = pr(r, u) P, (cos 8), (11.19) 


(1) 


where the functions ’ are defined as the solution of the equations 


Wt — 2p + Ark (HP — VO) —rAn(n +A) y= 0 (11.20) 


and describe 2"-pole radiation. Inturn, the solution of equations 
(11.20) has the form of a polynomial in r-?: 


ype = >»; L (u) r-+1) | 
k (k) 
the coefficients of the polynomial satisfying the recurrence relation 
qd ™ (n) 
2(k +4) L =(n—h)(n—k+1)L, 
4 C+) (k) 


In view of this analogy, we can regard the quantities @, charac- 
terizing the Riemann tensor as being a solution of a generally 
covariant wave equation of the Zel'manov type (7.4). 

This analogy allows us to define the multipole moments of an 
axially symmetric system of sources in the following way: 


| (@: + B) P® (cos 6) sin 0d, (11.21) 


0 


1 


monopole (mass) M =— - 


dipole moment D = ——-\ @®P! (cos 6) sin @ d8, (11.22) 
if) 


nm 


\oeP2 (cos 6) sin 640, 


0 


4 


Quadrupole moment Q = — —- 


and so forth. 
For the Bondi metric (11.3) we obtain from (11.21) . 


M = My + 4-\(CO)osin 6 a0, (11.23) 


0 
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where Mg denotes the Bondi mass of the system, given by formula 
(11.12). Analogously, the law of conservation of the Bondi mass 
(11.13) becomes 


dM, 
du 


= —) C0 sin 648 + <- \(CO)oosin 6 a8. (11.24) 


Definition (11.21) should be regarded as more general than 
Bondi's definition (11.12), as it is not linked to the choice of a 
particular system of coordinates. According to the interpretation 
given by Newman, the quantity M can be treated as the total mass 
including the energy of radiation, whereas the quantity Mg, is the 
mass of the system ina stationary state (''Weyl'' mass). 


3. Gravitational radiation of arbitrary isolated systems. 
The Sachs metric 


We saw that the tetrad components of the Riemann tensor for the 
Bondi metric break down into three groups, the dominant terms at 
large distances being of order r-! in the first group, of order rin 
the second group and of order r-° in the third group. Thus it has 
proved possible to express the Riemann tensor as the sum of three 
terms, proportional to the first, second and third inverse powers, 
respectively, of the affine parameter r: 


Raps =r. Napys + r?. [Tlagys =f Lays -O (r-*). (1 1,2 5) 


It was found that Nagy, apy and Iggy are tensors having the 

Same algebraic properties as Rieman tensors belonging to types JN, 
III and I, respectively. An analogous decomposition was established 
by Robinson and Trautman for their solution of the equations of 
gravitation, which describes spherical gravitational waves: 


Rapys — Tr Nays + r [Lopes + r Days ; (1 1.2 6) 


here Dg gy is a tensor having the same algebraic structure as a 
type D Riemann tensor. It should be mentioned that the physical 
meaningfulness of these results is further substantiated by the fact 
that for electromagnetic radiation (as Goldberg and Sachs have 
shown /184/) we have a decomposition of the field tensor Fy, which 
is clearly interpretable in terms of near zone (induction zone) 

and far zone (wave zone), In other words, the general physical 
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analogy between the theory of electromagnetic radiation and the 
Bondi-Sachs approach described above indicates that it is legitimate 
to interpret the Bondi-Sachs results in terms of a wave zone of 
radiation (terms of types WN and III) and an induction zone (terms 
of types I and D). 

Indeed, the dominant terms at large distances from the emitting 
system belong to type N (wave zone), whereas at small distances 
from the sources they are terms of type lor D, which describe 
the properties of a stationary field (of the Weyl type for the Bondi 
metric and of the Schwarzschild type for the Robinson-Trautman 
metric). Here the term Ngg,ys for the Bondi metric is proportional 
to Coo, characterizing the relationship between the information 
function Cy, and the asymptotic behavior of the radiation field. 

Sachs /185/ put forward an extension of (11.25) to the case of the 
gravitational radiation of arbitrary isolated systems, 

It is easily shown that the coordinate lines g in the Bondi 
coordinate system are trajectories of the Killing vector field /58/ 
normal to a certain three-dimensional hypersurface. Assuming 
that the metric tensor g, can depend not only on u,r,6@ but also on 
the coordinate g, let us construct a more general coordinate 
system suitable for an arbitrary system of sources with no spatial 
symmetry whatsoever. 

Let u(z*) be a scalar field defining the isotropic three- 
dimensional hypersurface (2.15), i.e., satisfying the eikonal 
equation (2,22): 


gu .Ug = 0. 

Further, let there be a congruence of isotropic geodesic lines with 
direction vector J, =u,, orthogonal to the hypersurface u= const, 
We also introduce 6 and g —two scalar functions satisfying the 
equations 

Gol = 9,17 = 0, (11.27) 
Let us further define the scalar r 

rt — (K siney}, (11.28) 
where 


K = (g768 8 9) (g¥ wP,v) — (g290,.8.5)? > 0. 


The condition that K +0 is a consequence of equations (11.27); the 
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requirement that K >0 is sufficient in order for the scalar r to be 
real and positive. 

It can be shown /185/ that an asymptotically flat space-time V, 
will satisfy all of the requirements mentioned if and only if its 
metric is representable in Bondi coordinates as 


ds? = Br-* exp (2B) du® + 2 exp (2B) du dr — 
— Hy (dx* — Atdu) (dx” — A” du) (11.29) 
(a, b = 2, 3), 


where the two-dimensional quadratic form h,,dz*dz® is given by 


2haydx* dx” = [exp (2y) + exp (25)] do? + 
+ 4 sh(y — 8) sin@dédg + 
+ [exp (— 2y)+ exp (—2 4)] sin’ 6 dq’, 


and the six functions V, U*(a =2, 3), 6, 7, 6 depend on all four 
coordinates u,r,6@ and gq, 

Having required that the coordinate lines @(u, r, 6 =const) be 
trajectories of a certain normal Killing vector field, we obtain as 
a Special case the axially symmetric Bondi metric defined by the 
conditions 


y=6, ui =0, 0g,,/d9 = 0. (11.30) 


4. Geodesic rays. The peeling theorem 


One of the assumptions under which the Sachs metric (11.29) 
was obtained is the existence in V, of a congruence of isotropic 
geodesics (known as geodesic rays). We will therefore con- 
sider the types of gravitational fields that admit of geodesic rays. 

As we noted in Chapter 8, every empty space-time V, admits of 
at least one and at most four isotropic vector fields 2 
(Debever vectors) satisfying the algebraic relation (3.29) — (3.32), 
in which the Weyl tensor is replaced — in our case — by the 
Riemann tensor Rgays- 

Debever vectors are not aS a rule the direction vectors of con- 
gruences of isotropic geodesics. However, as Sachs shows /110/, 
all algebraically special empty spaces V,, i.e., Petrov type D, HU, 
N and III spaces, admit of a geodesic isotropic vector field I, 
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1218, = 0, (11.31) 


satisfying the condition (3.31) and hence necessarily the condition 
(3.32). The corresponding congruence of geodesic rays has zero 
shear (Goldberg-Sachs theorem, see /184, 195/): 


2? = lay li — 5 (If)? = 0. (11.32) 


If there exists a solution of the eikonal equation (2.22), then 
condition (11.32) means that there is no distortion of the shape of 
the shadow cast upon a screen by an opaque object oriented ortho- 
gonally to the rays. ' Depending on whether the expansion of the 
congruence 


+h, (11:39) 


e= 
equals zero or does not equal zero, the wave front will correspond 
to plane or spherical gravitational waves. * 

This suggests that algebraically special fields correspond to 
gravitational waves far from the source system while 
algebraically general fields represent a gravitational field near the 
sources, ''perturbing" the gravitational wave front. To substantiate 
this conclusion mathematically, let us examine the asymptotic 
behavior of the Riemann tensor of an asymptotically flat space- 
time described by a Sachs metric. 

Assuming that the functions that enter into the Sachs metric 
(11.29) are infinitely differentiable with respect to r (which, as we 
will see, is not generally obligatory) and expanding them in Taylor 
series in powers of r!, we can show by direct calculation with the 
field equations (2.2) that all the @, in (11.15) are nonzero and that 
their expansions begin with terms of varying orders, namely from 
r+ to r-5>, Taking the general Newman formulas (11.16) — (11.17) 


* If the space-time admits of an isotropic geodesic congruence with nonzero shear, the 
corresponding field of gravitational radiationisinterpreted as a field of cylindrical waves, 
Examples of axially symmetric gravitational fields of this type with an infinite linear distribu- 
tion of sources are the Einstein-Rosen solution /187/, and also the Jordan-Ehlers /188/ and 
Kompaneets /189/ solutions that develop it. Unlike the plane and spherical wave fields con- 
sidered earlier, they belong to the algebraically general type | and, like gravitational fields of 
isolated sources (Bondi and Sachs), admit of peeling into terms of type N, II and | (see the 
works of Stachel /188/ and Marder /190-192/). Researches by Weber and Wheeler /193/ and 
Krishna Rao / 194—196/ have been specially devoted to the wave properties of nonstationary axially 
symmetric gravitational fields. 
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for the asymptotic form of @M, into account, we arrive at the follow- 
ing peeling of the Riemann tensor (Sachs' "peeling 
theorem /110/: 


R = Nr-? +4 ollIr-* 4 olIr-? 4 gIr-t 4 Ir 4+... (11.34) 


Here, ,N, ,Ill,... denote tensors with the algebraic structure of 
the Riemann tensors of types JN, III, etc., respectively (indices 
have been dropped for conciseness), and the index 0 on the left 
means that these tensors are covariantly constant along the corre- 
Sponding geodesic rays. In the fourth and fifth terms, I' and ol are 
distinguished by the fact that the tensor I' has no geodesic rays 
while pl has geodesic rays. Further, the sum of the terms up toa 
certain given order inclusive is a tensor of the same algebraic 
type as its last term. The coefficient ,V is proportional to Cy, 
where C, is a function which goes over into the Bondi information 
function under the conditions (11.30). In particular, for C,»=0 for 
the Bondi metric ,N = ,III = ,ll =0, and we arrive at the stationary 
axially symmetric Weyl metric (11,10). 

On the basis of the peeling theorem, the gravitational fields of 
isolated sources can be given the following algebraic interpretation. 
A field set up by a material system in the empty space surrounding 
it belongs to the first Petrov type (Ior D). Ifa system emits, 
then at distances significantly greater than the dimensions of the 
system itself and the wavelength of its emission (i.e., in the wave 
zone) the gravitational field will be approximately of type N. In 
other words, from the standpoint of an observer lying at a large 
distance (in a fixed direction) from the emitting system, terms of 
type N will predominate in the Riemann tensor; for an observer 
whose distance from the system is small compared with the dimen- 
sions of the system itself and with the wavelength of its emission 
(near zone), terms of type lor D will predominate in the Riemann 
tensor depending on the nature of the wave front. Lastly, at 
distances great by comparison with the dimensions of the source 
system but small compared with the wavelength of its emission 
(transitional zone), the gravitational field is described by a Rie- 
mann tensor of type II or III depending on the nature of the source 
distribution. 

As shown by Sachs /148, 185/ and Persides /197/, given certain 
special assumptions (possibility of expansion of the metric in 
infinite series in r-!, Euclidean asymptotic form at infinity), the 
Riemann tensor peeling (11.34) is feasible for arbitrary gravita- 
tional fields in empty space possessing geodesic rays. This 
possibility is guaranteed by conditions (11.16) ~— (11.17) obtained 
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above which give the asymptotic form of the tetrad components of 
the Riemann tensor. For cases where the space-time does not 
admit of a congruence of geodesic rays, Lehman /198/ and 
Goldberg /199/ formulated and proved the following weaker state- 
ment: in these cases one of the isotropic Debever vector fields is 
asymptotically geodesic, i.e., 


laiel®=O(r™) (n> 2). (11.35) 


If we assume that the peeling (11.34) still holds for this case, it 
follows that the vector /* no longer remains the direction vector 
for the geodesic ray in Spaces corresponding to tensors JN, III, 
etc., of the expansion. However, there should then exist a con- 
gruence of isotropic geodesics with tangent vector l’/* which are 
asymptotes of the trajectories of the vector [*. One might then 
expect that the congruence /’* represents geodesic rays only ina 
Space-time corresponding to the first four terms of the peeling 
(11.34): 


ltpRatarps lo) 2° 1% = O (1°), (11.36) 


i.e., that it is asymptotically a Debever vector. 

Lastly, if the vector field l’* satisfies condition (11.36) and 
moreover is a gradient field (i,’=9,), then, as Newman and Pen- 
rose demonstrate /174/, the peeling (11.34) always holds. 


5. General algebraic structure of the Riemann tensor 


It is obvious that the Sachs peeling formula (11.34) is not gener- 
ally covariant since the parameter r serves as distance coordinate 
along the geodesic. However, the algebraic classification of 
canonical types of Riemann tensors makes it possible to formulate 
the peeling of a Riemann tensor of algebraically general structure 
ina generally covariant manner, 

To do this we use the Newman-Penrose formalism, In an iso- 
tropic complex tetrad (I¥, n', m*, m“) whose vectors satisfy 
the relations 


Ll? = ngn* = m_zm* = 1,m* = n,m* = 0, 
(11.37) 


[,n* = m,m* = 1, 


we define three simple bi-vectors 
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Vag = 21M}, Vag = 272M}, (1 1 .38) 
Mas = 2l tale] aa 2ffi[aMg}- . 


Then, as shown by Kammerer /200/ and Szekeres /134/, to all 
possible algebraic types of the Riemann tensor in a vacuum corre- 
sponds the following general combination: 


Raprs +t *Ragys = CiVagV ys + Ce (VagMys + 
+ MapVys) + Cs (MagM ys + UagVvs + VapU ys) + 
4 Cc (UspMs +My Ua) +C,Uagt ys, (11.39) 


where C,, ..., C; are arbitrary scalars. 

Using the canonical form of the matrix || Ral| of the Riemann 
tensor in the bi-vector space R,s, one can show that there is, to 
each type of gravitational field (Ch.3), an expansion of the Riemann 
tensor in the bi-vectors (11.38) corresponding to some particular 
case of expression (11.39), where some of the scalars (Cj, ..., Cs 
(not more than four) vanish. The problem of expressing the 
Riemann tensor in terms of bi-vectors for the case of empty 
spaces (and the Weyl tensor in the general case) was solved for 
all types of gravitational fields by Debever /81/. Comparing 
Debever's results with expression (11.39), we arrive at the follow- 
ing conclusions, Let a given gravitational field be algebraically 
special. Then C,=C,=0, i.e., the expansion (11.39), generally 
speaking, comprises only the first three terms; in this case the 
vector field (* defines geodesic rays. Further, let C,;=0; the 
Riemann tensor then belongs to type III. For the case C,= Cs= (y= 
= C,= 0 the gravitational field is of type N with the vector [+ as 
Debever vector. Lastly, the condition C;#0 characterizes fields 
of type Il and D. 

Comparing expression (11.39) with the Sachs peeling formula 
(11.34), we see that each of the first five terms of the asymptotic 
expansion of the Riemann tensor in r-1 has the algebraic structure 
of the corresponding term in expression (11.39). It is therefore 
of interest to study the asymptotic behavior of the coefficients 
Cy, ...,C, in (11.39). To do this we choose a canonical parameter r 
on the geodesic ray with tangent vector J#. If we assume that the 
condition r— oo always corresponds to the asymptotic value of the 
Riemann tensor, then in the general case there is no need to 
assume that the space-time is asymptotically flat. Moreover, 
instead of the rigorous Sachs condition regarding the infinite differ- 
entiability of the metric with respect to the parameter r, we con- 
fine ourselves to the assumption that the components of the Riemann 
tensor, as also those of the vectors [*, n*, m*, m*, are functions of 
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the coordinates of class C* (which corresponds to our aim of 
investigating the asymptotic behavior of the first five terms of the 
expansion in r-*), 

We insert expression (11.39) in the Bianchi identities and pro- 
ceed to regard the latter as equations for the functions (,..., 
C,. Studying the main (by order of smallness) terms of the solu- 
tion of these equations, we find that for r—-o the coeffici- 
ents Cy in (11.39), N=1, 2, 3, 4,5, tend to zero as 
l/r’, respectively. 

Thus the Sachs peeling theorem can be proved under far more 
general assumptions. A rigorous proof of this kind-was first 
given by Kammerer /200/; as he assumed C*-smoothness, his 
proof covers only the first four terms in (11.39). This corresponds 
to the fact that geodesic rays exist only for the first four terms in 
(11.39). Kammerer's result can be extended to the fifth term using, 
Say, considerations relating to asymptotically geodesic congruences, 


6. Asymptotic symmetries. Bondi-Metzner group 


From the peeling theorem we have seen that at large distances 
from the source system even the properties of asymptotically plane 
gravitational wave fields are algebraically complex. The nontrivial 
character of the asymptotic forms of Bondi and Sachs fields of 
gravitational radiation is particularly striking when one studies 
their asymptotic symmetries, 

Indeed one might have expected from physical considerations 
that on the asymptotic hypersurface r-oo there wouldacta groupof 
motions of the metric of the same order as the order of maximum 
mobility of the space V, of the corresponding type, i.e., of Petrov 
type N. However, due to the fact that at large distances the 
radiation field of an insular system is described only approxi- 
mately by the structure of the type N Riemann tensor, it turns 
out that there is no group of motions that conserves the asymptotic 
form of the field and the boundary conditions as well. However, 
there exists a group of coordinate transformations satisfying these 
requirements; for axially Symmetric isolated source distributions 
it is known as the Bondi-Metzner group /171/. This group contains 
the Lorentz group as subgroup (but not a normal one) and includes, 
moreover, the infinite-dimensional group of 'supertranslations." 

The fact that the information function C, has very Simple trans- 
formation properties with respect to the Bondi-Metzner group 
plays a very important part in the physical interpretation of this 
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group. It makes it possible to express the conserved integral 
quantities in terms of the group invariants. An extension of the 
Bondi-Metzner group to the case of arbitrary isolated systems was 
put forward by Sachs /148/. 


7. Asymptotic properties of the Einstein-Maxwell fields 


The Bondi-Sachs method has also been extended to the case of 
gravitational radiation by isolated systems in nonempty space- 
time (Kozarzewski /201/, Hawking /202/, Stachel /203/). In 
particular, Kozarzewski /201/ shows that the asymptotic behavior 
of a gravitational field induced by arbitrary isolated systems of 
electrically charged bodies is also determined by the Sachs peeling 
formula, This result seems natural since geodesic rays are the 
trajectories of propagation of gravitational and electromagnetic 
radiation alike. 

The analogy between the gravitational and electromagnetic fields 
is also clearly discernible in their asymptotic behavior. Thus 
starting from the integral form of the Maxwell equations Goldberg 
and Kerr /152/ established that the electromagnetic field of a 
bounded distribution of charges and currents admits of the following 
asymptotic expansion: 


Fyy =r No + Pt, + ry (11.40) 


Here r is an affine parameter which varies along the isotropic 
gradient directions I", i.e., along the electromagnetic rays; all the 
components Jy, = Jt. are functions bounded above and the anti- 


symmetric tensors N*, and III%, satisfy the relations 


Nib =0, lhl =a (a is a scalar), (11.41) 


in complete analogy with the algebraic relations characteristic for 
the Riemann tensors of types WNW or III, respectively. 

It follows from formula (11.41) that in a coordinate system in 
which the parameter r characterizes the distance from the system 
of sources, the field F,, becomes isotropic at a distance from the 
system, i.e., satisfies the relations 


. 


InP ay; = 0, I¢F, = 0. (11.42) 


However, as Chevreton /149/ shows, conditions (11.42) are 
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necessary but not sufficient for plane electromagnetic waves in 
Minkowski Space-time; in order for the electromagnetic field F,, 
to fit the case of plane waves, it is necessary and sufficient that 
the following be satisfied in addition to (11.42): 


lp Fey, « = 0, l¢ Fog g = 0. (11.43) 


Consider conditions (11.42) — (11.43) in the general case of a 
curved space-time. As Mariot /260, 261/ shows, the trajectories 
of the vector field I« satisfying relations (11.42) form a congru- 
ence of isotropic geodesics (11.31). Differentiating the relations 
(11.42) covariantly and recalling (11.43), we obtain 


Freydazo=9, Faglis=0, (11.44) 


whence it follows /149/ that the tensor ],, can be expressed as the 
product of two vectors: 


h;e = Aglh,, (11.45) 


where A, is a vector orthogonal to the vector Jl. It follows 
automatically from formulas (11.45) that for the congruence of 
geodesics defined by the vector /*, the rotation w, the dilation 
e andthe shear so tendto zero. This means that at infinite 
distances from a system of arbitrary charged sources the electo- 
magnetic and gravitational fields are fields of plane waves with 
common congruence of geodesic rays. 


Chapter 12 


GRAVITATIONAL WAVES AND CHRONOMETRIC INVARIANTS 


1, Chronometric invariants 


The preceding chapters were devoted mainly to the study of 
gravitational fields from the standpoint of generally covariant 
criteria for the existence of gravitational waves. A gravitational 
field which satisfies the generally covariant wave criterion will 
have a wave character irrespective of the choice of coordinate 
system. But it would be interesting to drop the requirement of 
general covariance andformulateacriterionfor gravitational - 
inertial waves covariant only with respect totransformations 
of the three-dimensional coordinate frames in which a chosen 
reference body is at rest. Further, sucha criterion would have 
to be invariant under transformations that conserve the coordinate 
lines of time z° since these are world lines of the reference body. 
In other words, this criterion reflects the observer's choice of a 
frame of reference and its fulfillment may therefore be regarded 
as a Sign of the reality of the waves for the given choice of refer- 
ence body; indeed such waves can generally be eliminated by 
passage to a different reference body, due to the dependence of the 
inertial properties on the state of motion of the observers. 

Let us then consider the transformations given by 


7 7 (2? gt 2), : (12.1) 


z' = # (2), 2*, 23), (12.2) 


zt 
which have the property that ss =0. The transformations (12.1) — 


(12.2), where z® is the timelike coordinate, are obviously the most 
general transformations between coordinate systems at rest 
relative to a given reference body. 
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In formulating the criterion we will confine ourselves to 
invariance under the transformations (12.1), i.e., chronometric 
invariance, and to covariance with respect to the transformations 
(12.2), i.e., spatial covariance. Thus for a specified refer- 
ence body, freedom of coordinate transformations is restricted 
by the chronometric invariance (12.1) and the spatial covariance 
(12.2). Transformations (12.1) — (12.2) constituted the base on 
which Zel'manov constructed his formalism of chronometric 
invariants /204 — 206/. 

Chronometric invariants, i.e., three-dimensional 
physical quantities invariant under the transformations (12.1), can 
be regarded as observables in the general theory of relativity, 

i.e., quantities directly linked to physical measurements. The 
chronometrically invariant approach to gravitational-inertial waves 
is thus all the more interesting as the waves so defined can be 
treated as a subject of immediate physical measurement. 

Following Zel'manov /204, 205/, we introduce chronometrically 
invariant operators of differentiation, designating them by means 
of stars to distinguish them from the usual ones: 


4 


i 
V goo 


‘ Boi 
Oo, 0, = 4; — > Oo. (12,3) 


We also introduce the chronometrically invariant spatial metric 
tensor 


8oi & 
b= — Bat, bik — — gik, b = det |] dj» |. (12.4) 


For the chronometrically invariant gravitational-inertial 
force vector F; and tensor A;, of absolute rotation 
of the frame of reference & relative to the locally 
co-moving geodesic system %X,, we have the expressions 
/206/ 


P= % -aWw —"AV)), (1229) 
Ain = "OG Vay + Fi Vi; (12.6) 


where W and JV; are the scalar and vector potentials of 
the gravitational-inertial field, respectively: 


W=(1— V £00): Vi = — Boil V Bo0- 
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The chronometrically invariant tensor D;, of deformation 
velocities of the three-dimensional space of reference of the 
system 2% relative to the locally co-moving system 2%, is given 
by the expressions /205/ 


Dix =>", D™ = — 4-adit, D = "any. (12.7) 


Here D = Dj is the rate of relative volumetric expansion of an ele- 
ment of space, 

We also define the chronometrically invariant analogs of the 
Christoffel symbols andthe operation of chronometrically invariant 
three-dimensional covariant differentiation /205/ 


1 a ° e 
Ais =D" ("O:by, + “Odie — *Oxdi,), (12.8) 


*V,03" = "0,03" — Abou +... + ARO (12.9) 
with 


*Vibix = 0, *v, by = 0, v0 = 0; 


Within the formalism of chronometric invariants, in addition to 
the dynamic quantity F;, kinematic quantity A;, and static quantity 
D;,, itis also possible to introduce a fourth geometric character - 
istic of the co-moving three-dimensional space, namely the spatial 
curvature tensor Kyj,j;: 


Kyi = (Haiti + Heriny)s (12.10) 
where 
Hii? = 0, Mh, — *O;Ahy + ATiAdm — AttAlm, (12,11) 
and 
Kyeiy = — Kany = — Kingi = Kiytn. (12.12) 


As Zel'manov shows, the twenty independent components of the 
four-dimensional Riemanntensor R,,,, can be collected in three 
chronometrically invariant tensors expressed interms of F;, Aix, Diz 
and Ki.;. Thus we introduce the following threé-dimensional 
tensors /91/: 
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x#o Rove! yuk = Roe ziti _ pili (12.13) 
ao View’ | 


It is easy to show that the tensors X4, Yue and Zi are chrono- 
metric invariants. Indeed, let Qi? be the components of a world 
(four-dimensional) tensor of rank n in which all the upper indices 
are nonzero while all the lower indices (m in number) are zero. 
Then, performing the transformation (12.1), it is seen that the 
quantities 


pier — Qo: (12.14) 
(go0)™? 


constitute a chronometrically invariant three-dimensional tensor of 
rank n—m. Note incidentally that formula (12.14) can serve as an 
algorithm for the construction of chronometric invariants of the 
form Qg*:? out of world tensors. It is obvious from the definitions 
(12.13) that the tensors X#, Yk and Zi were constructed according 
to this precise rule, and they therefore satisfy the condition of 
chronometric invariance. 

Expressing the components of the world tensor Ry, in terms 
of the chronometrically invariant quantities (12.5) — (12.7), (12.10), 
we obtain the Zel'manov formulas defining the relation between the 
tensors (12.13) and these quantities: 


« * 4 
Xj = "OD y— (Di + Ai) (Dj + An) +°VeFy— FF, (12.15) 
Y ijn = °Vi (Aix + Dix) —°Vi(Aje + Dux) — 24,jF i, (12.16) 
Zin = 2(Ditx Dy; — Aite Any + AigAnr) — Kini. (12,17) 
Here 
Xia we yee at Zt, 4 x4 = — RY (12.18) 
800 V &0 : : 
Xiy= Xz, Yi =— Yur, Yar =9, (12.19) 


and the tensor Z,,;; has the symmetry and antisymmetry properties 
(12,12) of the tensor Ky;;. We see that the twenty independent 
components of the Riemann tensor Rygy can be expressed in terms 
of the six independent components of the tensor X,,, eight of the 
tensor Y;;, and six of the tensor Z,);;. 
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2. Chronometrically invariant definition of 
gravitational-inertial waves 


Thus the Riemann world tensor R,g,3 breaks down into three 
chronometrically invariant three-dimensional tensors (12,13), 
expressible in turn in terms of the gravitational-inertial physical 
characteristics of the three-dimensional space co-moving with the 
chosen reference body, (12.5) — (12.7) and (12.10). We will relate 
the definition of gravitational-inertial waves to the four gravita- 
tional-inertial characteristics (12.5) — (12.7) and (12.10), and also 
to the quantities (12.13) which are expressed in terms of these 
according to formulas (12,15) — (12.17). 

A chronometrically invariant definition of the d'Alembert oper- 
ator in the general theory of relativity, realizing Zel'manov's idea, 
has been formulated in /94/. 

The chronometrically invariant criterion for 
the existence of gravitational-inertial waves is 
that the following three-dimensional chronometric 
invariants — the vectors F,;, the tensors Dyp, Ain, Knit; 
and the scalars composed of these, as well as the 
tensors Xij, Yijn Zrii Expressed in terms of these — 
satisfy equations of the form 


OP =@. (12.20) 


The following notation is introduced here: 
Oa — tara, V2 =bV,"V,, (12.21) 
a 


so that *C] is the chronometrically invariant spatially covariant 
d'Alembert wave operator and a is a scalar function of the coordi- 
nates, It is assumed that @ does not explicitly contain second 
derivatives of the sought function P. The latter will be called 
wave function according to tradition, in the sense of equation (12.20). 
The role of the wave function P will be played by various . 
chronometrically invariant quantities of a three-dimensional- 
tensor nature, Investigating the chronometrically invariant wave 
criterion reduces to analyzing whether it is fulfilled for various 
chronometrically invariant characteristics of the frame of reference 
and gravitational field with respect to this frame, Accordingly it 
is necessary to distinguish between gravitational-inertial waves of 
force F;, deformation D;, and curvature K,,;;. 

In /94/ this criterion was used to analyze a‘number of known 
solutions of the Einstein equations in a vacuum, It was shown that 
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it is met for all known solutions of Petrov type N (solutions of 
Peres /160/, Takeno /153, 163/, Petrov /65/ et al.) but not for 
the "cylindrical wave’ solutions of Einstein-Rosen /187/ and 
Kompaneets /189/, which are not of type NM. Note that the latter 
two solutions (Petrov type I) fail to satisfy any of the generally 
covariant criteria for gravitational waves considered in the pre- 
ceding chapters (criteria of Lichnerowicz, Bel, Pirani, Zel'manov, 
et al.). 

Let us look at wave equations of the form (12.20) for gravitational 
fields corresponding to certain exact solutions of Einstein's 
equations in a vacuum (2.2), 

The Peres solution (9.30) — (9.31) belongs to Petrov type N. 

In order to simplify the calculations we impose the additional 
condition* f; = f,.=0 upon f asa result of which the only field 
equation (9.31) is identically satisfied. 

Using formulas (12.3), we have the following for an arbitrary 
function P ofargument (2 + 2%): 


‘ fot ot 4 
OP = Vi-—2z Ps, Pa Gay Pa 


eae 4 
o*oP = qa ap lP. (1 — 2f) +f 3P.s3l, 


‘ 4 
“Os 0,P = Tapp l Pas (1 = 2f) + 2f 3P.3]- 


From the above, expressing P,, in terms of *@*@P and Pj, in 
terms of *9,*0,P, and noting that 


*V?P = Dik ("0;0,P — Aix'0;P), 
we reduce the ordinary wave equation for the function P(z® + 2°) 
Pys3 — Pio = 0 


to the required form in the given frame of reference 


* In actual fact this condition means that we are confining ourselves to a flat space-time, since 
the nonzero components of the Riemann tensor for the metric (9.30) are given by 


Rsu3 = Rsrio = Rou = Roop =fry1, R323 = R220 = Ro220 = frees 
Rais = R320 = Rois = Roo = fr 


Thus our case will fit inertial waves only. The more general case of gravitational-inertial 
waves in a space-time with R,,/, 0 will be considered later for different exact solutions 
of the Einstein equations. 
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(*V? — *9*9) P = 0. (12.22) 


Since the three-dimensional scalars 


=a 2 = Groh 


for this metric likewise depend on the argument 2x°+ z3 only 
(Ai, = 0 for f, =f, =0), it follows that they satisfy equation (12.22), 
i.e., that theyare solutions of the chronometrically invariant equation 
(12.20) for Q@=0, a=1. 

Let us now obtain the wave equations for the gravitational- 
inertial force vector F‘ and the deformation velocity tensor D,;, in 
the given reference system. The field of the vector F,; is given by 


Fy=0. F=0, F=—_455; 
computing the quantity *CIF;, we obtain the form of the right-hand 
side of equation (12,20) provided that it does not contain 
derivatives higher than the first orderin F, If such derivatives 
are present on the right, it obviously follows that the equation is 
not a wave equation; this is the situation, for example, for the 
Einstein-Rosen and Kompaneets metrics. 

The results of calculations for the Peres metric lead to the 
following chronometrically invariant spatially covariant equations: 


OF = — 3FCV,Fi + 4 FF), (12.23) 
*C) Diy = — 2D;; (3°OD + 2DaD"). (12.24) 


Equations (12,23) — (12.24) are obviously chronometrically in- 
variant equations for a wavelike gravitational-inertial field, the 
nonlinearities on the right representing the sources of the gravita- 
tional-inertial perturbations, 

The Takeno solution (10.5) also belongs to Petrov type N. Upon 
calculating the tensor ofthe angular velocity of absolute rotation (12.6) 
for the metric (10.5), wefindthatit vanishes: A;,=0. But, as 
Zel'manov shows /206/, this is necessary and sufficient to make 
all g); vanish in this region by the transformation (12.1). There- 
fore, assuming from now on that § =0 inthe metric (10.5), we find 
that the scalars 
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; P 
FiFi=—75, 
4 
= ———, (PM ,+ MP 
. amp ( 8 :s) 


satisfy the chronometrically invariant wave equation 
OG = 5h PVGVM, 


where @G is any of the chronometrically invariant scalars FF, and 
Di. 

For the gravitational-inertial force vector the wave equation 
becomes 


‘1 Ft = Fi(— 2F,Fi + 3°0D + 2D;,D*" — D*) + 
+ ayy DEV M (2F)F) — *0D — DD), 


from which it follows, in particular, that the sources of wave per- 
turbations of the gravitational-inertial field depend on the deforma- 
tion of the frame of reference. 

The Petrov solution (9.15) was examined ina different coordi- 
nate system by Bondi, Pirani and Robinson /143/; Synge /172/ 
interpreted it in the language of ''body gravitational waves." 

We will write as follows: 


ds? = dx®*— da'*+ adz®* 4 28 dx*dz3 + y dx3?, (12,25) 
where a, B, y are functions of the argument 2° +2’ related by a 
single differential equation /65/. For this solution, inthe frame of 
reference (12.25), Fi=0, Ai, =0, and the scalar D is given by 
D=--"d In (or — 8) 
and satisfies a scalar wave equation of the type (12.20) for g=1: 
‘OD =D (be *V,D°V,D)'r. 

We turn now to solutions of Petrov type 7,. It appears that the 

wave criterion under consideration — unlike other generally 


covariant criteria — can be satisfied by certain type 7F, metrics 
as well, Thus the metrics 
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ds? = —a'dx®* 4-adz!? + ydz?* + ysh?2*dz5’, 
ds? = — a1dzx°* + adr}? 4 ydx?* 4+ y ch? 2? dz? 
where a and y are functions of the argument 27° +21(a2<0, y< 0) 
satisfying a certain system of differential equations /94/, define 
gravitational fields of type Dp /207/. The three-dimensional 
space ofthe frame of reference is holonomic (A, = 0), and the 


chronometrically invariant scalars D and F;F! satisfy the wave 
equation (12.20) for a=a: 


‘VG — = 0'aG = bE VG, Int 


A solution which merits special discussion is that of Einstein- 
Rosen /187/, which belongs to type I: 


ds? = e%(v-a) (dx9? — dx”) — x1 %e-2a dg?? — et dz’, (12.26) 


where ¥ and @ are functions of z! and 2° satisfying the differential 
equations 


Oa + oO, — Ot,—9 = 0, (12.27) 
Ya = 2" (1)? + (4,0)7], Yo = Axia sao. (12.28) 


It can be shown /94/ that although the "cylindrical waves" equa- 
tion (12.27) admits of chronometrically invariant formulation (12.20) 
for a=1, 


"Cla = bik*Vi0°V;, (a — 7) + "daar — (°02)2, 


neither it not an analogous equation of the type (12.20) with a differ- 
ent right-hand side and different a@ will be satisfied by the chrono- 
metrically invariant scalars D and F,F!( A;,=0), or by any of their 
scalar functions. 

An analogous result holds for the Kompaneets metric /189/, 
which generalizes the Einstein-Rosen metric: 


ds = 0, dx®® — ada!* —y de®* — 28 dr%da® — §da*?, (12,29) 


where a,f,y,6 are functions of z! and 2°, From the field equations 
one can arrive at two equations for "interacting cylindrical waves": 


[2 (PQ — 1)"“P 1},1 — 2 (PQ — 1)""Polso = 9, (12,30) 
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[21(PQ —1)"4Q.tla — (PQ —1)"40, 01,0 = 0, (12.31) 
where the following notation is used: 
P=7(76— 217)", Q = 8 (v5 — 21)", 


As in the case of the Einstein-Rosen metric, the system of equa- 
tions (12.30) — (12.31) admits of chronometrically invariant spati- 
ally-covariant formulation (12.20) for a =1: 


Cy = gas (GA) n+ or}, (12.32) 
118 = gg (BY — 2) @+ 264, (12.33) 


where A = 75 — B*, and a and o are certain functions of «, B, y, 6 
and their first derivatives. However, equations of the (12.20) type 
are likewise not satisfied by D, F;F! or by their scalar functions. 
To close this section we will discuss the question of the general 
relationship between the chronometrically invariant criterion 
(12.20) and the generally covariant Zel'manov criterion (Ch. 7) 
examined in /165/, Writing down the system of 20 equations 


8° Rapys; 00 =0 
in chronometrically invariant spatially covariant form, we arrive, 


after fairly extensive calculations, at the following three systems 
of equations: 


(‘V2 — *a*a) Xi = Ais (six equations), (12,34) 
(1 

(‘V2 — *a°a) Yiit = Av® (eight equations), (12.35) 
(2) 

(*V? — *a'a) Ze — A*s (six equations), (12.36) 
(3) 


where the right-hand sides A4, Auk, At4éi are chronometrically invari- 
Q) (@) (3) 
ant spatialtensors ofthe second, third and fourth rank, respectively, 


and contain no derivatives higher than the first order in the "wave 
functions’ X#, Yuk, Zeii, These equations are given in expanded 
form in Appendix II. 

Thus any space-time V, satisfying the Zel'manov 
generally covariant criterion of gravitational 
waves will also satisfy the chronometrically in- 
variant criterion of gravitational-inertial waves. 
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The role of the wave functions in the corresponding wave equations 
of the type (12.20) is filled by the chronometrically invariant 
tensors (12.13), i.e., by the quantities xX, yuk zeuj, which repre- 
sent the Riemann tensor in the given frame of reference. 


3, Physical conditions of existence of gravitational- 
inertial waves 


Since equations (12.20) are not generally covariant, the waves 
described by these equations are closely related to the physical 
characteristics of the chosen frame of reference and more specifi- 
cally to the chronometrically invariant quantities F,;, A;, and D;, 
(12,5) —(12.7). Let us investigate what part these physical 
quantities play in wave equations of the type (12.20) and how they 
influence the existence of gravitational-inertial waves. It is clear, 
for instance, that the choice of a frame of reference based on the 
form of the quantities Fi, A;, and D;, may lead to restrictions under 
which the tensors X#, yt and Z*'4i become stationary, i.e., time- 
independent. Then the chronometrically invariant d'Alembertian 
(12.21) in equations (12.20) degenerates into a Laplacian, testifying 
to the absence of physical gravitational-inertial waves in the given 
frame of reference, 

Moreover, inthe chosen frame of reference the three-dimensional 
Laplacian °V? entering the expressionfor *C] (12.21) may make the 
"wave function" vanish. This situationmaybe realized, for example, 
whenever the functions P in(12.20)arehomogeneous: *V; P = 0. 

Thus both of the above situations (stationarity and homogeneity 
of the wave function) can be regarded as sufficient (but not gener- 
ally speaking necessary) conditions for the absence of gravitational- 
inertial waves in the specified frame of reference, 

By expressing the identities (7.8) in chronometrically invariant 
notation, one can easily show /208/ that in the Einstein spaces (3.7) the 
chronometrically invariant representatives of the Riemann tensor 
Xi, yt Ziti always satisfy the wave equation (12,20) for a=1 
and a definite choice of the right-hand side of Q. Thus the 
question of the existence of gravitational-inertial waves for the 
quantities x, y#*, 77 in Rinstein spaces reduces merely to in- 
vestigation of the nontriviality (in the above sense) of the left-hand 
sides of these wave equations. 

In /208/ a study was made — not just for Einstein spaces but in 
the general case as well — of the sufficient conditions under which 
a frame of reference does not admit of gravitational-inertial 
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waves, i.e., the left-hand side of equation (12.20) necessarily 
becomes degenerate. A complete classification of frames of refer- 
ence not admitting of gravitational-inertial waves (either due to 
stationarity of the wave function of due to its homogeneity) is given 
in this work for all chronometrically invariant quantities playing 
the part of wave functions, 

Let us quote the results of this study for the case of an arbitrary 
gravitational field in a medium with the energy-momentum tensor 
T,,- We introduce, in accordance with Zel'manov /205/, the 
concepts of chronometrically invariant density, pressure and stress 
tensor of the medium: 


P=TyolGo, Jt=TVV 80, UX =TH. (12.37) 


Let the part of wave functions be played by the chronometrically 
invariant representatives X*), yis* Z'*4 of the Riemann world tensor. 
We will investigate the wave equation (12.20) for these functions in 
frames of reference such that: 

a) all chronometrically invariant mechanical characteristics 
of the frame of reference (12.5) — (12.7) tend to zero, 

b) one of these is nonzero, 

c) two of these are nonzero, and 

d) all three chronometrically invariant mechanical character- 
istics ofthe frame ofreference are nonzero: F'+0, A;, #0, 

D;, ~ 0, and we will determine which frames of reference fail to 
admit of gravitational-inertial waves, due to the homogeneity or 
stationarity of the wave functions, 

Let the following homogeneity conditions be fulfilled:* 


ViFi=0, “Vidi =0, °VjDn=0, °VjKix = 90, 
(12.38) 
“di = 0, ViU ix = 0, “V,J# = 0. 


One can show that in this case all the wave functions are homo- 
geneous, i.e., 


*(1p = — *a*aP 


for any of the tensors (12.13). Thus when the homogeneity condi- 
tions (12.38) are fulfilled, there are no gravitational-inertial 
waves, 


" The conditions (12.38), first pointed out by Zel’manov, differ from the homogeneity conditions 
*d 
which he gives in /205/ in that in (12.38) the equations ri =0, *V;B,, =0 have been 


replaced by the equation *VjU4,—=0, and the equation *Vjq; =0 by the equation *V;J;-=0. 
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Assuming next that the space is inhomogeneous, let us investi- 
gate another sufficient condition for the absence of gravitational- 
inertial waves, namely, stationarity of the wave functions. We 
begin with case (a), where 


F=0,A;,=0, Di, =0, (12.39) 


i.e., the chosen frame of reference fallsfreely, does not rotate and 
is not deformed. * The quantities (12.13) then become 


X#—0, Yi—0, gieu — — Kieu, (12.40) 


It appears that the third condition in (12.39) leads to stationarity of 
the three-dimensional curvature tensor Ki, and consequently 
there are no gravitational-inertial waves in this frame of reference. 
Conditions (12.39), which define the frame of reference, 
make it possible to recover uniquely the general form of the space- 
time metric V,. Indeed, joint fulfillment of the first and second 
conditions of (12.39) means that in the given frame of reference it 
is possible to parametrize the time lines 2° so as to have simul- 
taneously /205/ 


fo = 1, Goi = 0. (12.41) 


The third condition (12.39) will then guarantee stationarity of the 
three-dimensional metric tensor (12.4), According to Cotton's 
result (see /65/, p. 389), the three-dimensional metric 6,, may be 
transformed in this case to a diagonal form. Thus in order for 
there to exist a frame of reference satisfying conditions (12.39), it 
is necessary and sufficient that the given V, be a reducible space 
of the special type 


ds* = dx? + gy,dx'* + goodx®? + gygdzx9?, (12.42) 
Su = Bu (z', a x). : 
It can be proved that the Einstein space (3.7) with metric of the 
type (12.42) is always flat (see /65/, p. 390). 
Consider now the case when the second of conditions (12.39) is 
not fulfilled: 


Fi=0, Ainp#O0, Din =9, (12.43) 


* Here, and henceforth when speaking of the free fall, rotation and deformation of a frame of 
reference, we will be referring to the corresponding motions of the three-dimensional space of 
the given frame of reference with respect to the locally co-moving geodesic system. 
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i.e., theframeofreference rotates while falling freely and without 
being deformed. The first condition in (12.43) permits us to choose 
a parametrization of the time lines z° such that the following 
conditions be fulfilled /205/: 


Zoo = 1, *0go; = 0. 


One can show that in a frame of reference with the properties 
(12.43), all gin are also stationary. From this it follows that the 
metric ga, of the space-time V, is stationary and therefore that the 
quantities (12.13) are stationary too. Thus a frame of reference of 
the type (12.43) will not admit of gravitational-inertial waves. The 
well-known Gédel metric /209/ is an example which fits this case. 

Next, consider the case where the first condition in (12,39) is 
violated: 


Fi +0, Ai, = 0, Di, = 9, (12,44) 


i.e., the frame ofreference accelerates without rotating or becoming 
deformed. The space-time metric can then be reduced to the form 


ds? = oq (x°, x', x*, 5) dx°* + g,,dz1* + goodx®*+ B53 02° a 


Gu = Bu (x', x*, 2°). 


(12.45) 


The quantities (12.13) in this frame of reference assume the form 


Kies > (*ViF) + VIF) — FIFI, (12.46) 
Yi — 0, (12.47) 
Zikij — __ Kirti (12.48) 


Due to the third condition of (12.44) Zz} is stationary and thus 
the problem reduces to investigating the wave properties of X#, 

Using the chronometrically invariant form of the field equations, 
we can write X4# as 


Xo = (Ki + Abi) + * (obi + 204 — DDK, (12.49) 


where A is the cosmological constant, 

From this it follows that waves of X are impossible ina vacuum 
owing to the stationarity of this ''wave function,'' even though as a 
rule the four-dimensional metric gg is then nonstationary, 

In the general case (7,, *«0) it follows from the field equations 
and from the conservation laws that the mass density p is 
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independent of time in the frame of reference under consideration. 
But the stress tensor Us is generally nonstationary. Thus the non- 
stationarity of X‘/ is due in this case to the nonstationarity of the 
stress tensor YJ, and the question of the absence of waves ina 
medium needs to be analyzed in detail for every choice of T,,in 
the field equations. 

Let us first require that the given frame of reference be co- 
moving with the medium. In this case the stress tensor may be 
expressed as 


U4 = phii — Bu = pedi — all, (12.50) 


where £* is the first viscosity, which develops in anisotropic 
deformation, ais the second viscosity, which develops in iso- 
tropic deformation, p is the true pressure, and pq) is the equi- 
librium pressure, determined from the equation of state. Since 
the viscosity of the medium is not felt in the co-moving frame of 
reference in the absence of deformation, it follows that U¥ = po¥ , 
and, therefore, that p = po). 

If the medium is baroclinic, i.e., if po = po(p,t), where + is 
the absolute temperature, then 


"OX" = ——> “ae bh, (12.51) 
and thus the X are generally nonstationary. In the case of a baro- 
tropic medium, for which pq = po)(p), the frame of reference under 
consideration does not admit of gravitational-inertial waves in 
view of the stationarity of the stress tensor. We will show that in 
a barotropic medium the four-dimensional metric tensor Sap iS 
also stationary. 

Under the above assumptions the equations for the relativistic 
conservations laws become 


‘ip <0, (12.52) 


1 
Sp UP = — ln V Boo. (12.53) 


Here equation (12.53) isthe chronometrically invariant analog of 
the condition of equilibrium in hydrodynamics extended to the case 


of the gravitational field /210/. Equations (12,52) — (12,53) together 
lead to the following expression for gy: 
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Zoo = exp 21T (x®) +R (z')l, (12,54) 


where R and 7 are arbitrary functions of their arguments. 
With the help of the coordinate transformation 


dz? = [exp T (2°)|dx®, dz! =dz' (12.55) 
the metric can be reduced to the stationary form 


ds? = exp [2R (z')j dg°? + fu dz? + Bia dz?” + Bas az * 
S _ (12.56) 
Gu = gu (2, 2, F°). 


(An example of such a frame is that in which the Schwarzschild 
metric is usually written). The main conclusion to emerge from 
the case under consideration is that if the medium is barotropic, 
then the frame of reference co-moving with it will not admit of the 
existence of gravitational-inertial waves. 

In discussing the case of an arbitrary energy-momentum tensor 
our only assumption was that the frame of reference is co-moving 
with the medium. However, certain variants of the energy- 
momentum tensor for which this assumption does not hold are of 
interest. The first such example is the energy-momentum tensor 
of an ideal fluid: 


Tag = (0 +P) Uglg — PBag- 


It follows from the field equations that the corresponding frame of 
reference is co-moving with the mass, i.e., J'=0. From the 
chronometrically invariant formulation of the field equations one 
can show that the frame of reference must also be co-moving with 
the medium in this case. This means that in an ideal fluid, as in 
the previous instance, gravitational-inertial waves cannot exist if 
the medium is barotropic. 

Consider another example — the energy-momentum tensor of 
dissipative systems:* 


To, = (p + p) Ugls + P8ap + Ta g- (12,57) 
Inthe frame of reference corresponding to the notation of (12.57) 


viscosity does not develop in view of the third condition of (12.44), 


* The following considerations are valid assuming that the dissipative processes (viscosity and 
thermal conduction) are not too strong /210/. 
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and since the frame of reference is co-moving with the mass there is 
no flux of heat. In this frame of reference, therefore, T,., fora 
dissipative medium reduces to the 7,, for an ideal fluid discussed 
above. 

The third example is that of the energy-momentum tensor of 
the electromagnetic field (2.30). Let 


aS + tbe F py 


be the tensor dual to the Maxwell tensor. We set 


Soa (ea (12.58) 
V 800 "Veo 


where d‘ and h' are the chronometrically invariant stress vectors 
of the electric and magnetic fields, respectively. One can show 
/211/ that these vectors are related to the chronometrically 
invariant representatives of the energy-momentum tensor of the 
electromagnetic field by the following relations: 


p=L +a), ht bah”, d= Dad", (12.59) 
Ji a nid las (12.60) 
Us — ppd — (nth + ald’). (12.61) 


In order for the electromagnetic field F,, to be isotropic, i.e., in 
order for it to satisfy the relations 


FF =0,  Fy*F* =0, 
or, in chronometrically invariant form, 
b=d, b"d, =9, 
it is necessary and sufficient that the condition 
J=op (J = b*JJ,) (12.62) 
be satisfied /211/. It follows from (12.62) that a frame of reference 
co-moving withthe mass," in which j'—0, cannot be realized in 


an isotropic electromagnetic field; from now on, therefore, we 
will consider the electromagnetic field to be anisotropic. 
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The tensor X“ for the electromagnetic field becomes 
xi® = (K* 4 Abily 4 ai, (12.63) 
whence it follows that 
*axi* — —1*0 (hin* + did*). (12.64) 


Let us determine for which anisotropic electromagnetic fields 
the tensor X* is stationary. In chronometrically invariant form 
the Maxwell equations are given by 


gimn ene ne Fn) (Rey V>) es, *a(d' Vd), 
gine (°Vin — Fy) (dn V>) =—*d (hi Vd), 


*v 07 = 2h"Q,, 
“Vink™ = 2d"Qn. 


(12.65) 


From the condition that j'=0 follows the collinearity of the vectors 
h; and d;: 


ad = ypi. (12,66) 


Setting *ad‘ = *ah' = 0 and taking (12.66) into account, we obtain, 
from equations (12.65), 


*Om Wn _— *0,¥hm3 (12.67) 
then 


h, = 6 (z’) “OnXs d, = % (z') S (x’) *Onk- ( 12, 68) 


But every vector /; proportional to the gradient satisfies an 
equation of the form 


eV dy 20: (12.69) 


Thus in an anisotropic electromagnetic field gravitational-inertial 
waves X* are not present if the chronometrically invariant stress 
vectors of the electric and magnetic fields satisfy condition (12.69). 
Lastly, let us consider the case where the third condition of 
(12.39) is not satisfied, i.e., where the frame of reference falls 
freely and does not rotate but undergoes deformation; sucha 
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system is termed semigeodesic (synchronous): 
Fix=0, Ay, =0, Diz #0. (12.70) 


According to (12.15) — (12.17), in this frame of reference the 
three-dimensional tensors (12.13) become 


xX — — DD's + Di Dik + (K4 + Abi) + 


yidt _ epi pit _ evi pie (12.72) 
Ziklj _ pik pi _ pil pri __ Kiki (12.73) 


and are generally nonstationary. 

Let us turn now to case (c), in which two of the conditions 
(12.39) are not met. Let the frame of reference accelerate and 
rotate without becoming deformed: 


Fi-~0, Ai~#0, Di, = 0. (12,74) 
The tensors (12.13) then become 


X4 = 345.AM + (K+ AdY) + 4 (pb + 20% — Ub), (12.75) 
yi _ svi gt eytgie 4 9 dipe (12.76) 
Zatti. Ait gd gitgts 4 oat gel Kins, (12.77) 


Making use of the identity 
"OA + *OF a == 0 (12,78) 


and recalling that in anundeformedframeofreference *0A;, = *dA**, 
we conclude, from (12.75), that the nonstationarity of X# is due to 
the rotational character of the field F' and to the nonstationarity 
of the chronometrically invariant representatives of the energy- 
momentum tensor. Similarly, from (12.76) and (12.77) it is clear 
that the nonstationarity of Yi is due to the rotational character 
and nonstationarity of the field Fi, and that the nonstationarity of 
Zi is due entirely to the rotational character of the field F', 
Thus for the case of rotational fields (*V,,F;,+ 0) a frame of 
reference of the type (12.74) does not preclude the existence of 
gravitational-inertial waves of all three chronometrically invariant 
representatives of the Riemann tensor. 

On the other hand, if the field Fi isirrotational, *V,,F;,=0, then 
it follows from (12.76) that Zjis stationary. The tensor Yuk is 
generally nonstationary in this case, since 
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*oVisk — QA art 0, (12.79) 


The tensor X is nonstationary only for 7.,;0. In the latter case 

the equations for the conservation laws in the frame of reference 

(12.74) lead to the conditions j#0, *@p #0 — unlike case (12.39). 
The condition that jé--0 means that the given frame of reference 

can exist not only in an anisotropic, but also in an isotropic electro- 

magnetic field. Let it be co-moving with the medium. Then 

Uik = pb in formula (12.50) and 


*axti = + ("dp — *dp) b¥, (12.80) 


In this case the function X will be nonstationary in a baroclinic 

as well as ina barotropic medium, unlike the third case considered. 
The special case where the barotropic medium is characterized by 
the equation of state p = p (medium of superhigh density, e.g., near 
a singularity) constitutes an exception; here *@X vanishes. Thus 
aframe ofreference characterizedby acceleration and nonstationary 
rotation does not preclude the existence of gravitational-inertial 
waves X#, Yik, Ziti, A frame ofreference characterized byacceler- 
ation and stationary rotation does not admit of the existence of 
waves X4 (in vacuum or in a medium described by the equation 
p=gp), nor does it admit of the existence of waves Ziti, 

For aframe of reference satisfying the demands (12.44) we find 
that in the case under consideration the functions X#, Yi, Zitli have 
precisely the same form (12.71) — (12.73) as in a semigeodesic 
frame of reference and are, asa rule, nonstationary. 

In the case where 


F;=0, AiO, Din 0, 
the functions X#, yi, Zi*i are given by 


xtk — — 34¢A' — DD* — Dip 4 (K™ + Abi) + 
++ (pot + 2U* — vot), (12.81) 


yisk = *V (Die ae A‘*) i +yi (Da* te Ai*), (1 2.82) 
VAL Di* ps — Di pki _ Kiki Aik Ali ALAR 4 9 Aid gel, (12.83) 
Since F, =0, it follows from the relation 


"0A, = Bs (Dia at bi™ DF) Ain (12.84) 
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that the nonstationarity of the "wave functions’ is due entirely 
to the deformation oftheframe ofreference. Unlike A‘*, the tensor 
A, is stationary in this case. 

Lastly, looking at the case where none of the three character- 
istics of the frame of reference vanishes, 


FF; 0, Aun #0, Da 9, 


we see that all three ‘wave functions" (12.15) — (12.17) are generally 
speaking nonstationary. 

We have been investigating the conditions sufficient for the 
absence of gravitational-inertial waves characterizable by chrono- 
metricallyinvariant representatives" of the Riemann tensor, i.e., 
homogeneity and stationarity of the quantities (12.13). We note 
that the homogeneity of all of these simultaneously is caused by 
the homogeneity of space. 

Besides the waves X yis* and z**¥ studied here, it is also inter- 
esting to analyze gravitational-inertial waves for which it is the 
chronometrically invariant characteristics of the frame of reference 
itself that serve as wave functions: waves of deformation D,,, 
waves of rotation Ay, waves of gravitational-inertial force F,, 
and waves of curvature Ky,,;. The results of this analysis are as 
follows /208/: 

Waves of rotation A, are present for F,; = 0 and for any 
D,, and absent for F,; #0 if Dy +0 or (for Dy =0) *VyFy #0 
(rotational gravitational-inertial field). 

Waves of acceleration (force) F,; are generally 
present (i.e., the field F, is not stationary) for any A, and D,. 

Waves of deformation D, are present only if *V,D,, +0 
(inhomogeneity of the deformation) and for *9*ab, 0, irrespective 
of the properties of A; and F,. 

Waves of curvature rae are present in any detennied 
frame of reference ( Dy, + 0) 

This exhausts the list of gravitational-inertial waves in arbi- 
trary frames of reference. Further research might profitably 
concern the effect of gravitational-inertial waves on concrete 
physical systems. In particular, research on the effect of 
gravitational-inertial waves on systems of test bodies, as men- 
tioned by a number of authors /95, 212 — 214/, could bring to light 
possible ways of detecting gravitational waves in the laboratory. 
The chronometrically invariant approach could play an important 
role here as a means of describing the observables — physical 
quantities, measured by laboratory means. Accordingly it will be 
useful to discuss the problem of the experimental detection of 
gravitational waves. 


Chapter 13 


THE PROBLEM OF GRAVITATIONAL WAVES 
AND PHYSICAL EXPERIMENT 


1. Geodesic deviation of test particles 


The expression of the field functions in the language of observ- 
able quantities accessible to physical measurement is bound to be 
important in any attempt to compare the results of the theory of 
gravitational radiation with the experimental data. Such 
observables have already been considered in Chapter 11, namely 
the invariant geometric objects characterizing the Riemann tensor 
as a field function within the framework of a tetrad formalism. 

In another variant, the chronometrically invariant approach of 
Chapter 12, we introduced the chronometrically invariant '"com- 
ponents’ of the Riemann tensor. We will now consider the problem 
of the physical basis for the experimental detection of gravitational 
radiation. 

The most convenient way of observing the quantities which 
characterize gravitational waves is to record the motion of test 
particles. Let there be given a set of test particles moving along 
geodesic lines -zx*(s), where s is the arc length along the geodesic. 
Let zx* (s, v) be a one-parameter family of such curves, a change in 
the parameter vy corresponding to passage from one geodesic to 
another. We introduce the two vectors 


of which the first is the tangent to the geodesic and the second is the 
infinitesimal displacement for one particle relative to another, 
5 : 


n be the covariant differ- 


whichis orthogonal tothe tangent. Let 


ential of displacement along the geodesic, divided by the element ds. 
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The quantity 


Dy* 


ds? 


is called the geodesic deviation and is interpreted physically 
as the measure of the relative acceleration of two infinitely close 
particles moving along adjacent geodesics. The role of this 
quantity in the theory of gravitation is obvious from the well-known 
equation of geodesic deviation: 


Din" 4 RE gusut x¥=0 (13.1) 
as? B q 


(see, for instance, /60, 172/). 

Equation (13.1) shows that the relative acceleration of two near- 
by particles moving without being influenced by external (nongravi- 
tational) forces is completely determined by the physical components 
of the curvature tensor. Consequently if we perturb the components 
of the Riemann tensor we modify the relative acceleration of the two 
particles. Conversely, by observing the relative acceleration of 
the particles we can measure the perturbations of the Riemann 
tensor due to incoming gravitational radiation. 

Let us suppose that the test particles are connected by a spring. 
The gravitational waves can then be recorded from the oscillations 
of the spring. The characteristic frequency of the system may co- 
incide with one of the spectral harmonics of the gravitational waves; 
it will then be possible to observe even the weakest gravitational 
waves from the resonance. 

Now let the vector y* connect not two adjacent test particles 
but rather two infinitely close points of a piezoelectric crystal. 

The deformation produced in the crystal by a gravitational wave. 
incident upon it will induce an electric field inside it. The integral 
of this field strength is a potential difference the measurement of 
which gives the physical components of the curvature tensor induced 
by the gravitational radiation /95/. If the integral of equation 
(13.1) is treated as the Fourier transform of Rj, to nz, then 
measuring the geodesic deviation reduces to measuring a single 
Fourier component (namely the resonance component) of the entire 
spectrum of gravitational radiation. It is therefore of interest to 
determine which frequencies are possible in the gravitational 
radiation produced by various sources. 
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2. Possible sources of gravitational waves 


In view of the difficulties encountered in generating gravitational 
waves in the laboratory it has proved advisable in experiments to 
look for gravitational waves of cosmic origin. The problem then 
arises of arriving at some sort of theoretical evaluation of the 
possible energy-mass loss due to gravitational emission for differ- 
ent sources. Estimates of this kind have been obtained for binary 
stars (Cooperstock /33/, Forward and Berman /215/), for 
collapsed stars (Zel'dovich-Novikov /34/), quasars (Cooperstock 
/32/), pulsars (Weber /31/, Shklovskii /216/), neutron stars 
describing nonspherical pulsations (Thorne /30/), and also the 
metagalaxy as a whole (Wheeler /36/). Let us briefly review the 
available results. 

The first cosmic objects to be considered as a possible source 
of gravitational radiation were the binary stars /15/. The 
presence of a quadrupole moment in these systems provides 
grounds for believing that they lose energy, and thus the param- 
eters of the stellar orbits are changing, These changes can be 
recorded by astronomical observations. Using the Landau- 
Lifshits energy pseudotensor to evaluate the energy of quadrupole 
gravitational radiation of a binary system in the linear approxima- 
tion, one can show that the energy loss E by a system of two 
bodies of mass m, and m, travelling along circular orbits around 
a common center of inertia at a distance r one from the other, is 
given by 


dE 32k { mm \24 4 
— = a) 8 


where k is the Newtonian gravitational constant and o is the 
circular frequency of revolution. From here the rate at which 
the bodies converge due to radiative energy loss is easily 
determined: 


6448 mym2 (m1 + ma) 
= ors 


v= r = 

Subsequently a formula was also obtained for the mean flux of 
energy emitted by a pair of bodies of mass m, and m, revolving 
along elliptical orbits of semimajor axis a and eccentricity e: 


Tr 2,2 

1 (dE 32 WAmymey (rm + ma) ( (9-25 1 so 

a\ ata ee) + xe}, 
0 


13. GRAVITATION AL WAVES AND PHYSICAL EXPERIMENT 139 


where YT is the period of revolution. Expressions were also ob- 
tained /217/ for the energy loss of binary systems and their gravi- 
tational radiation spectrum in the case of hyperbolic motion, free 
fall and other types of motion. 

As in the case of the problem of energy, research has been 
carried out into the problem of the loss of gravitational momentum 
by binary stars due to their gravitational emission (Cooperstock 
/33/). Within the limits of the linearized theory of gravitation it 
was shown that binary stars may lose gravitational momentum 
(defined via the energy-momentum pseudotensor) with an intensity 
of the same order as energy. A formula was obtained expressing 
explicitly the functional dependence of the total momentum flux on 
the difference in phase of the components of the binary system. 

According to computations by Forward and Berman /215/, for 
binary neutron stars the theoretical maximum emission rate 
is P~6- 104 watt, and it does not depend on the sum of the masses 
of the pair components. 

Other than binary stars cosmic sources of gravitational radiation 
may consist of cosmic bodies (planets, asteroids and so forth) 
falling on collapsed stars. It is easy to compute (in the 
linear approximation) that if such a source lies at distance of 
~ 500 Mpc from the Earth, then for, say, m=m@ and M=10? mo 
(M being the mass of the collapsed star and m the mass of the 
body falling on it in a straight line or along a spiral), near the 
Earth the flux of power emitted may amount to 0.7 erg/cm? sec 
/34/. However, it is not known how frequent such processes are in 
the metagalaxy. Basing themselves on rigorous solutions of the 
field equations, Forward and Berman /215/ propose to refine the 
estimate obtained by Zel'dovich and Novikov /34/ for the gravita- 
tional power emitted by collapsed stars from the linearized theory 
of gravitation. 

Thorne concludes /30/ that neutron stars which pulsate 
in nonspherical manner can give rise to considerably 
stronger gravitational emission (compared with binary stars). In 
contradistinction with the approach usually employed to the prob- 
lem of the gravitational radiation of binary stars, in which the 
stars are represented as material points and the calculations are 
carried out in the linear approximation, in Thorne's work use is 
made of concrete models of neutron stars, the calculations being 
carried out by an approximation method in which the Schwarzschild 
solution is the point of departure. 

On the average the period of pulsation 7 of a neutron star lies 
within the range 10-* —10-? sec. If we assume that the pulsation 
energy is entirely converted into the energy of gravitational radiation, 
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then given certain logical assumptions regarding oscillation ampli- 
tude, a neutron star may, according to Thorne, emit energy of the 
order of 10°! erg in one period of pulsation. These estimates 
agree as to order of magnitude with similar ones by Wheeler /238/, 
although the latter are valid only in the linear approximation. 
Thorne and his colleagues have carried out a general study of 
the nonradial pulsations of stellar models in the Einstein theory 
of gravitation /218 — 222/. We will briefly describe their method. 
Let the unperturbed line element 


(ds?)) = e’dt® — erdr? — r? (dO? 4+. sin? 0 do?) (13.2) 
describe a spherically symmetric equilibrium configuration 


Vv=0(N, N= AC), 


The gravitational field of a perturbed configuration (a pulsating 
and, generally, revolving star) will then be given by the line 
element 


dst = (ds*)y + hap de® da?, 


where h,, is the perturbation of the metric tensor relative to 

the stationary metric (13.2). The displacement & of an element 
of ideal fluid from the equilibrium position, deviations of 
density and pressure from equilibrium values inside the fluid ( 6p 
and 6p, respectively) as well as perturbations hg of the metric, 
are expanded in spherical harmonics /218/. The perturbed metric 
is written as follows: 


ds?= eX(1 + Hy Yh) dt? +2H,Y dt dr—e(1 —H,Y',) dr? — 
—r?2(1 — KY!,) (0? + sin? 0 dg’), (13,3) 
where the functions #H,(t,r), H,(t,r), H(t.r) and K characterize 


perturbations of the metric, and Yi, (0, g) is an ordinary spherical 
harmonic: 


2i+4 (L—m)! Th pm : 
Y7, (8,9) = a Tear 1 (cos 8) em, 


Inserting the metric (13.3) in Einstein's equations and retaining 
only terms linear in the perturbations hg, we can solve these 
equations by numerical integration. Machine calculations for 
Specially selected models of neutron stars jf 22%, 222/ have made it 
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possible to determine the frequency spectrum of their gravitational 
radiation, pulsation periods (~ 1073 sec), pulse powers (~ 10° erg 
per sec) andthe decay times of perturbations due to emission 

(~ 1 sec). 

Cooperstock /32/ investigated the transfer of energy by gravita- 
tional radiation from quasi-stellar sources (quasars). 
He uses the Fowler model /223/, according to which a quasar 
consists of a core formed by two collapsed stars and an outer 
quasi-stable shell; the transfer of energy by gravitational 
radiation from the rotating core to the shell leads to "polar 
explosion," giving rise to an intensive movement of stellar material 
in the direction of the axis of rotation of the core. If the shell 
describes radial oscillations, then the power of resonance gravita- 
tional radiationof the core calculated with the Landau-Lifshits 
pseudotensor will exhibit an angular dependence of the form 
~ sin? 96 (the angle 6 gives the direction of radiation relative 
to the core's axis of rotation). Thus radiation is absent in the 
direction of the rotational axis and reaches a maximum in the 
equatorial plane. In cases where the shell is also rotating, 
the power emitted will also depend on the coincidence or 
noncoincidence of the directions of rotation of the core and 
the shell. 

Estimates have been made of the intensity of gravitational 
radiation of cosmological origin. Since gravitons are 
absorbed only very weakly by matter, once interstellar gravita- 
tional radiation had appeared (e.g., in the primordial ''explosion" 
of the Universe) it might have persisted to the present day, its 
total intensity in the metagalaxy depending on the rate of expansion. 
According to Wheeler's calculations /224/, the present rate of 
expansion of the metagalaxy gives an energy density of gravitational 
radiation for the metagalaxy of 107/19 =78 g/cm, i.e., a power 
flux of 10° erg/cm*sec, with a period of 7 ~ 108 years. 

The intensity of the gravitational bremsstrahlung of 
the Sun due to chaotic thermal (nonrelativistic) atomic collisions 
—a quantity not previously investigated — has been computed by 
Carmeli /35/. His computation is based on application of the 
classical method of Fourier integrals to the Landau-Lifschits 
expression for the intensity of gravitational radiation in the weak- 
field approximation, Carmeli obtained a sensational result accord- 
ing to which the power of the Sun's gravitational bremsstrahlung is 
of the order of P~6- 10 erg/sec, four orders higher than the 
power of gravitational radiation due to the total quadrupole moment 
of the planets of the solar system. Carmeli also obtained the 
frequency spectrum of gravitational bremsstrahluhg of arbitrary 
systems of interacting particles. 
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In addition to gravitational bremsstrahlung, mention must also 
be made of the possibility of gravitational emission by the Sun due 
to nuclear explosions as well as to the thermal motions of atoms 
generating a relative quadrupole moment. However, it is 
very difficult to determine the scale of nuclear explosions on the 
Sun (Braginskii and Rudenko give a few estimates of their power 
/158/). 

Today among the strongest sources of cosmic gravitational 
radiation are believed to be pulsars, including one possible 
model, that of neutron stars pulsating witha period of 10-*— 10-3 sec. 
If we assume that the energy of pulsation is entirely converted 
into energy of gravitational radiation, then — as we noted earlier — 
within a single pulsation period a neutron star may emit energy of 
the order of 10°! erg or about 0.1% of the rest mass of the star 
itself /30/.* This result was used by Weber /31/ to evaluate the 
order of magnitude of the Riemann tensor components 


XA Rs. 


generated by the gravitational radiation of pulsars. It was found that 
for a pulsar mass of M =10°8 g, the lower bound of values of 2;%? 
is of the order of 


Rij > 5.10- cm 2, 


3. Devices for the laboratory detection of gravitational waves 


The gravitational radiation field has been recorded experi- 
mentally with the aid of the quadrupole harmonic mass -detector, 
the first example of which was built in 1964 at the University of 
Maryland by Weber, Zipoy and Forward /212/. The sensitive 
element of the instrument is an aluminum cylinder weighing 
1.5 tons, suspended in a cylindrical vacuum chamber on a metal 
filament. Where it touches the filament the cylinder is sheathed 
in a piezoelectric quartz envelope connected to a sensitive volt- 
meter in the radio receiver system. After certain improvements 
introduced by Sinsky /226/, the sensitivity of the mass detector 


“ For other pulsar models the power emitted may be considerably less. Thus for pulsar NP 0532 
inthe Crab Nebula, within the framework of the model of a rotating neutron stat whose rota- 
tional axis is inclined to the axis of symmetry of its magnetic field, the gravitational power 
emitted was estimated to be of the order of 8x10" — 8x 10 erg/sec /225/. 
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was such as to enable measurement of relative displacements of 
the cylinder ends (e.g., expansions and contractions due to 
incident gravitational radiation) of the order of 10-18cm. An 
experiment for detecting gravitational radiation using two groups 
of n identical, closely placed parallel cylinders was subse - 
quently suggested by Braginskii. If the oscillations in these groups 
are excited in phase, the gravitational radiant power will then be 
roughly four times greater than from a single group. 

In a further refinement of his device, Weber developed a 
system which works on the principle of coincidence of signals of 
identical frequency from two detectors /227/. It consists of two 
mass detectors with a relaxation time of 30 sec timed to the 
frequency w=104* rad/sec and placed at a distance of 2 km from 
each other. During 1967 this system recorded coincidences 
(accurate to At =0.20 sec) of signals with a mean periodicity 
of one per month. In Weber's opinion it is extremely unlikely 
that the recorded coincidences were purely random. Weber 
assumes that the signals he detected were due to cosmic 
gravitational radiation. He demonstrated that to detect gravita- 
tional radiation from pulsars using this device for a measurement time 
of the order of one month, it is sufficient to record effective dis- 
placements §~3-107' cm. The level of modern experimental 
technology may thus be high enough to allow detection of cosmic 
gravitational radiation. 

Increasing the separation of the detectors, one increases the 
sensitivity of the device. Accordingly Weber later employed six 
detectors /228/, one in the Argonne National Laboratory and the 
rest in the laboratory of the University of Maryland. The distance 
between these laboratories amounts to 1,000 km. The detectors 
were adjusted to the expected gravitational radiation frequency of 
collapsed supernovas in our Galaxy — 1,660 Hz. The signal coin- 
cidences recorded by this instrument during a period of several 
months nearly preclude the possibility of attributing them to random 
coincidences. Weber interprets these results as proof of the 
existence of powerful gravitational radiation in the Galaxy. 

Weber later performed a second series of experiments /229/ 
which confirmed his original results. A special feature of these 
experiments was maximum insulation of the apparatus from 
external influences of electromagnetic and seismic nature. 

The possibility of the existence of gravitational radiation in the 
Galaxy in the frequency band around 1,660 Hz impelled a number 
of researchers to seek new sources of gravitational waves in the 
Universe. According to calculations by Greenstein /230/, cosmic 
gravitational radiation at 1,660 Hz may arise fromclose encounters 
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of stars in massive stellar clusters with a mean periodicity of once 
per week. Here the pulsed character of the radiation makes it 
likely that it is due to binary encounters of unstable relativistic 
objects — neutron or collapsed stars. 

One would expect gravitational radiation to have an important 
influence on the evolution not just of individual stars or stellar 
clusters, but of the Galaxy as a whole. Thus the gravitational 
radiation flux of the Galaxy recorded by Weber does not preclude 
our attributing the observable expansion of the Galaxy to the loss 
of its mass ensuing from the emission of gravitational waves 
/231/. However, as many researches show /231 — 235/, this 
interpretation is far from indisputable, and leads to contradictions 
with other astrophysical observations. That the interpretation of 
Weber's experiments is not unique has been demonstrated in works 
by Braginskii, Zel'dovich and Rudenko /236, 237/. 

Seismic methods which make it possible to use the Earth as 
detector may play an important role in the problem of the detection 
of gravitational waves. This possibility seems attractive in that 
the quadrupole moment of the Earth is many orders higher than 
laboratory detectors. The frequency of the Earth's eigenoscilla- 
tions (of the order of 1 millihertz) makes it possible to record the 
resonant harmonics of gravitational radiation from pulsars. How- 
ever, as Weber has shown /31/, this method is limited by the high 
noise temperature of the terrestrial core and thus would require 
measurement of effective displacements 6 ~ 2.107!" cm, which 
lies at the very limit of possibility of modern measuring techniques. 

A more effective approach may be to use individual seismically 
insulated inhomogeneities on the Earth's surface capable of absorb- 
ing gravitational radiation in the frequency band around 1 hz. In 
the linear approximation of Einstein's theory, the problem of the 
reaction of an elastic body to incident gravitational waves was 
investigated by Dyson /258, 259/. He found that the absorption of 
gravitational waves by an elastic body takes place exclusively 
due to inhomogeneities in its shear modulus, such absorption 
being absent ina homogeneousmedium. According to estimates by 
Dyson, the intensity of seismic signals due to gravitational 
radiation from probably theoretical models of pulsars at 1 Hz 
frequency is five orders lower than the noise level. However, the 
possibility of seismic recording of gravitational radiation of pulsars 
cannot be regarded as definitively closed. In particular, DeSabbata 
/239/ has recently suggested using local inhomogeneities on the 
Moon's surface — ''mascons'' — to detect gravitational waves from 
pulsars at the same frequency (1 Hz; see also /240/). 

Mironovskii /241/ suggested using a variant of the Weber mass 
detector to detect gravitational radiation from binary stars. The 
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receiver is a friction-free torsion pendulum of period 7). The 
formula expressing the gravitational power emitted by a system of 
two material points moving along circular orbits about a common 
center of gravity gives the following expression for the spectral 
density of gravitational radiation of binary stars: 


9 (T) = 2Nf (27) g (7). 


Here N is the number of stars of a given type in the Galaxy and 

T = 2n/m is the period of the gravitational wave (T= 1/2, where 

t is the period of revolution of the star's components about the 
common center of gravity; due to the equivalence of the components, 
the frequency of the gravitational radiation is doubled). 

In recording the emission of binary stars (fairly widespread in 
the Galaxy), the most suitable stars belong to the WUM class, 
totalling N ~10® inthe Galaxy. For stars in this class the 
function p (7) has a sharp maximum for 7 = 04,15 (d standing for 
"day''). The gravitational radiation spectrum corresponds to 
periods of revolution of 0,1 — 0.5, 

The energy emitted by close pairs of binary stars of 
class WUM can be evaluated using the Landau-Lifshits pseudotensor, 
Knowing the function p (T) and integrating numerically over periods 
with N =10°, Mironovskii obtained 10°® erg/sec for the total 
gravitational power emitted by this type in the Galaxy, which is 
only five orders less than the corresponding electromagnetic 
radiation. Within the Solar System the flux should then be about 
10-7 erg/sec cm’, 

In order to explore the possibility of recording this radiation, 
Mironovskii studied the equation of motion of a torsion pendulum. 
If a pendulum in equilibrium is oriented along the x-axis of an 
orhtogonal geodesic coordinate system, then the equation of geo- 
desic deviation for motion in the oscillation plane zy of a point of 
the pendulum lying at the distance |] from the rotational axis is 
given by 


yo oy == c*l >) Reo19 (7; ;) Sin (22;t + Q;). (13.4) 


Here Rj, is the Riemann tensor component corresponding to 
the radiation field of the i-th star, Q; is its angular velocity of 
revolution, and 7r; and m,; are vectors characterizing the position 
of the star and the orientation of its orbit in space, respectively. 

Solving equation (13.4) and averaging over the distribution of 
the stars in space, orientations of the orbits and receiving device 
and distributions of phases and frequencies, we obtain /241/: 
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(92) = 10-2979" VF (2T,) t [cm], (13.5) 


where /(t) is the probability density for distribution of WUM stars 
over periods of revolution, YT, is the period of eigenoscillations 
of the pendulum and ¢ is the observation time. 

Formula (13.5) gives the mean deviation g of the pendulum from 
the x-axis and is independent of the choice of energy pseudotensor. 
Mironovskii gives the following estimate for two WUM stars near us: 


(72) = 6-107* cm. 


As we know from the works of Weber and Braginskii, it is now 
possible in practice to measure periodic mechanical displacements 
up to an order to 10-18 cm. One may therefore hope to see 
Mironovskii's suggestion realized experimentally. 

Besides the principle of deviation of geodesics for test particles, 
used in experiments with the quadrupole mass detector, 
researchers have recently turned to another principle by means of 
which the interaction of the electromagnetic and gravitational fields 
can be exploited for the experimental detection of gravitational 
waves. 

As we remarked in previous chapters, there exist a number of 
rigorous solutions of the gravitational field equations describing the 
propagation of gravitational and electromagnetic waves along the 
same trajectories. The existence of such gravitational fields 
follows from the general solution of the Cauchy problem for Ein- 
stein's equations, according to which the characteristic hypersur- 
faces of the equations of gravitation and electromagnetism (wave- 
front surfaces), and their bicharacteristics (trajectories of wave 
propagation) as well, coincide. Asa result a field of electro- 
magnetic radiation can induce gravitational waves. It is natural 
to expect that the converse effect — stimulation of dynamic electro- 
magnetic fields by gravitational waves — is also present (e.g., as 
a result of the action of gravitational radiation on the field of a 
system of charged bodies). 

This effect was studied by Heintzman /242/, who examined the 
motion of test particles in a gravitational plane wave in the weak- 
field approximation. As the velocity of an uncharged test 
particle is not altered by interaction with a passing gravitational 
wave, it follows that a system of electrically neutral test particles 
cannot absorb gravitational radiant energy. Hence the idea 
suggests itself of using a system of charged particles and recording 
the gravitational radiation according to the variations of the 
electromagnetic energy of the detector system, Heintzman 
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computed the amount of gravitational wave energy absorbed not only 
for the case of weak fields in the linear approximation, but also for 
two exact solutions of the field equations: the metrics of Einstein- 
Rosen cylindrical waves and Bondi plane waves. He found that in 
allthree cases the electromagnetic energy of the system of 
charged test particles is altered by the absorption of energy of 
gravitational waves: in the first and third case (plane waves in the 
linear approximation and Bondi plane waves) the amounts of 
absorbed energy are finite, while inthe second case (Einstein- Rosen 
cylindrical waves) the detector can receive an unlimited amount of 
gravitational wave energy. These results are sufficient grounds 
for using a system of charged test particles as laboratory detector 
of gravitational waves. 

A slightly different method of detection was considered by 
Vodyanitskii and Dimanshtein /465/. They obtained a solution 
of the Einstein-Maxwell system of equations for the weak field in 
the form of a plane monochromatic wave. The following formula 
was put forward for the power of the gravitational: signal (received 
by an electrical antenna with amplifier): 


4g (Fo)® (2 + @) (27))? | 
Px pS — aa 


here S is the effective surface of the antenna, 2 and o are the 
frequencies of the gravitational and electromagnetic waves, 
respectively, propagating in opposite directions along the Oz axis, 
and Ej, and hj are the amplitudes of the corresponding electro- 
magnetic and gravitational field components. 

A method for the experimental detection of gravitational waves 
using dynamical electromagnetic fields has also been considered. 
It is based on observation of the fluctuations produced by gravita- 
tional waves in electromagnetic radiation, 

A theory of perturbations (fluctuations) of the electromagnetic 
field due to the gravitational field was constructed in the linear 
approximation by Cooperstock /243/. He solved the equations for 
perturbations of the electromagnetic and gravitational fields in the 
case where a plane-polarized monochromatic electromagnetic 
wave of high frequency propagates between two ideally conducting 
parallel walls, interacting with a gravitational plane wave of low 
frequency. The fluctuations of the field were computed for two 
mutual orientations of the direction of wave propagation: where the . 
directions coincided ("longitudinal orientation''), and where they 
were mutually orthogonal ("transverse orientation"), Since it may 
be that macroscopic systems generate low-frequency gravitational 
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radiation, Cooperstock suggested that this radiation might be 
detected by observing fluctuations in the intensity of electro- 
magnetic radiation of celestial bodies. 

This problem was subjected to detailed investigation by Winter- 
berg /244/. He demonstrated that in the linear approximation the 
problem of the intensity fluctuations of a light signal in a medium 
with statistically distributed gravitational waves is equivalent to 
the problem of intensity fluctuations of a light signal in a medium 
with statistical distribution of inhomogeneities. On the average 
this effect is described by the formula 


ay-22 yay 18.0) 


To no 


where J, is the signal intensity averaged over time, AJ the devia- 
tion from the mean value /,, n, the mean value of the refractive 
index of the medium, An the deviation from nm, x the distance 
between the source and the observer, and 7 the characteristic 
dimension of inhomogeneities in density. Defining the effective 
value n(Q, g) of the refractive index of otherwise empty space 
filled with gravitational waves by means of the relation 


one can obtain the effective An and n, for the case of plane waves 
in a linear approximation, Then averaging with use of formula 
(13.6) gives us the expression 


1.6/4 \(Z)”. (13.7) 
Here / is the characteristic length of the gravitational wave, 
h=det |i hy |, 2» being a Small correction to the pseudo- Euclidean 
metric. 

The values of |h | and / were estimated for three types of 
Sources: binary stars, quasars and the entire Universe. (Inthe 
latter case, the author had in mind the relict gravitational 
radiation, under the assumption that at an early stage in the evolu- 
tion of the Universe it was in thermodynamic equilibrium with the 
relict electromagnetic radiation at a temperature of 3°K). For 
binary stars the coefficient s =1.6 |h|/-‘ in formula (13.7) is 
found to be k =5.9-107*, which for a distance z =10"%-cm yields 
fluctuations of radiation intensity of the order of 5.9- 107°, i.e., 
values which can be recorded by sensitive scintillation counters 
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outside the Earth's atmosphere. The possible gravitational 
emission of quasars was estimated to be k= 2.3-10-*, which for 
z= 1077 em yields values easily measured even under terrestrial 
conditions: A/J/I,~7.3,. Lastly, for gravitational waves of cosmo- 
logical origin fluctuations amount to ~ 0.5 for a star 10? light years 
away. However, the latter value is practically impossible to ob- 
serve since in observations of this kind the stellar diameter would 
have to exceed the length of the gravitational wave. 

Nontheless from the standpoint of principle Winterberg's method 
has not been given sufficient theoretical grounding. As Zipoy and 
Bertotti point out /245/, extrapolation of the formula (13.6) given 
by Scheffler /246/ for electromagnetism to the case of gravitational 
fluctuations is not admissible in the case of a strong gravitational 
field. Stellar scintillation estimates may therefore refer to purely 
coordinate effects, i.e., nonphysical and thus unobservable effects. 

We also note that the methods described above for direct mea- 
surement of gravitational radiation intensity with quadrupole mass 
detectors entirely disregard the possible quantum effects of the 
interaction of gravitational waves with the crystal detector, These 
effects have been investigated by some authors /247 — 254/. Their 
researches tend to confirm that it is possible to detect and generate 
gravitational waves by the techniques of quantum electronics. 

The plans for a resonant receiver of gravitational waves using 
the interaction of gravitational waves with the atomic structure of 
matter have been worked out in principle by Lavrent'ev /252/, and 
also by Kopvillem and Nagibarov /249,.250/. Energy is first 
stored in the receiver of directed gravitational radiation by optical 
stimulation, and the gravitational ray serves to create optimum 
conditions for release of this energy in a definite direction in the 
form of a coherent electromagnetic ray. 

Kopvillem and Nagibarov /248, 255, 256/ also investigated the 
possibility of generating directed gravitational radiation in the 
laboratory by exciting coherent periodic oscillations of mass 
quadrupoles in the electron shell of atoms with the help of lasers. 
The general idea of these works consists of creating a special, 
so-called ''superemissive’ state of matter, the stimulation of 
which by sequences of brief, powerful pulses might be used to 
generate gravitational rays. A superemissive state of matter can 
be induced by coherent fluxes not only of photons but also of other. 
elementary particles (electrons, neutrons, protons and so forth), 
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4. Relationship between the theoretical and experimental 
aspects of the problem of gravitational waves 


The foregoing methods for experimental detection of gravita- 
tional waves, aS we saw, are not free of flaws as far as theoretical 
grounding is concerned. Thus in Weber's calculations the equation 
of geodesic deviationis not related to physical observables, while the 
very concept of observables in the general theory of relativity is 
still far from unambiguous. Hence the necessity for choosing a 
certain preferred frame of reference, and this choice is never 
unique. 

Others among the methods considered are based on the linear- 
ized theory of gravitation, whose flaws we discussed earlier; on 
this argument rests, in particular, the critique of Winterberg's 
method given by Zipoy and Bertotti /245/. 

Let us notea few basic traits of the relationship between the 
theory of gravitational waves and experimental physics. 

The design of any experiment is necessarily based on definite 
premises without which interpretation of the experimental data 
would be impossible. Thus the process of measurement (compar- 
ing length with a standard) which underlies experiments is basedon 
a definite specification of a method for identifying objects of one type. 
From this point of view the principle of relativity which underlies 
the physical theory precedes experiment (see /257/). The theoret- 
ical formulation of the concept of gravitational waves may there- 
fore be regarded as the basic premise for the interpretation of 
experimental results. 

In this respect a significant drawback of the conception of gravi- 
tational waves considered earlier is their internal incom- 
pleteness, resulting in multiplicity and inconsistency of the 
criteria, The latter is evidence that the level of the theoretical 
concepts themselves is not high enough in this field. 

Thus with regard to the best known of these criteria — the 
Lichnerowicz criterion — the very formulation contains arbitrary 
elements and inconsistencies with the initial premises. Accord- 
ing to the basic premise, a gravitational wave is characterized 
by discontinuities in derivatives of the type gij,.9 on the character- 
istic hypersurface in that coordinate system in which this hyper- 
surface is expressed by equation (2.15). Yet one can point to gravi- 
tational fields satisfying the Lichnerowicz criterion for which the 
metric components gag are smooth functions of the coordinates, 
their derivatives nowhere exhibiting Hadamard discontinuities, An 
example is the Petrov metric 


ds? = 2dax°dx! — sh? 2° dx?” — sin? x° da”, 


18. GRAVITATIONAL WAVES AND PHYSICAL EXPERIMENT 151 
for which 
B 
Oa = 812a8 = Bia» 


i.e., p=2°+const. It is obvious that on the "wave front" surface 
(x° = const) the curvature tensor components 


Roses oe ae sh? rae Rosos = sin? Fad 


cannot have Hadamard discontinuities. 

Other covariant criteria for gravitational waves based on the 
concept of Hadamard discontinuity in the Riemann tensor com- 
ponents — discussed in previous chapters — suffer from the same 
flaw. 

In conclusion we can state that, from the theoretical as well as 
from the experimental point of view, the problem of gravitational 
waves is still very far from solved. But it remains one of the most 
pressing and fundamentally important problems not just of gravita- 
tion theory, but of modern physics as a whole, a problem to which 
researchers are devoting ever more attention. It is therefore 
reasonable to expect that it will find an elegant solution and 
assume its proper place in the rigorous picture of Einstein's 
theory of gravitation. 


Appendix I 


We will prove certain theorems used in the text. Let Rapes be the Riemann tensor 
of a space V4, antisymmetric in each of the index pairs eB and yd. We introduce the 
two tensors conjugate to it: 


1 Gee » 1 oe 
*Rasys ss ~T Najes Ree Rapes =D Nybe0 Reap" (1.1) 


It follows in an obvious manner from the definition (1.1) that 


¢ 


* Ro ays = RB gap: (1.2) 


Theorem 1. The following relations hold in Einstein spaces: 


*R ays = *R sap: (1.3) 
Proof. Due to (1.2) condition (1.3) is equivalent to 
wt 
*Rooys = Rapys: (1.4) 


Thus to prove the theorem it is sufficient to prove (1.4). In order to do so we will 
make use of the relations (/172/) 


i» 
nour Noah =— Oey? (1.5) 
1 Necaa = — 2803" (he) 
ie! Bee = 66k, (1 7) 


where the tensor OF is a generalized Kronecker symbol obeying the following 
tules: if u, v, a... are all different and A, f, y,... are obtained from them by 
a certain permutation, then it is equal to + 1 depending on whether the permutation 


ae... is even or odd; in the remaining cases it is zero. Multiplying both sides of 


1 


*R vk =a. Nooap Bae 


153 
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by —- - nero, we obtain, owing to relations (1.6) 


re | = 


1 6 € nap. Yee 
— are be ee = ane Sein a RY a 
Contracting in the indices e and v, we obtain from this 


1 Vv 
R= —37 1 lea Ln (1.8) 


Multiplying (1.8) by Txpy? Placing 4 by 6 and using (1.5), we have 
1 SYST 
Naasy RE = ie Oopy F oovs — 
1 


a FR avs a *Rovas se * Rags o *R pons = *R a ios C2 *R p08) = 3* Fr aayys- 


Consequently, in the Einstein spaces (3.7) 


3* Rio ovjs = MyapyOh = *Msapy « m8) 
Rewriting (1.9) as | 
3*Raasyy = sap 
and adding, we obtain 
2* Raps = *Ravas + *Aasay + *Rasyp + *Rypas + 2 MNsapy - (1.10) 
We interchange the indices q and y, as wellas B and 4, in (1.10): 
i a LT + Rega 4 “Raat sas Lar + 2x14. (1.11) 


Comparing (1.10) and (1.11), we see that relations (1.4) hold, which proves the theorem. 
In conformity with this theorem, it is customary to set the following in Einstein 
spaces: 


Rae 


* — 
Rays ~~ aby apys°’ 


(1.12) 
Lemma. Inspaces V, the equations 


I* Roo, = 0 (1.13) 
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are equivalent to equations 


LR aor +l Roays + Ig Ry avs = 0. (1.14) 


Proof. Let there exist a vector /* in V, satisfying equations (1.13). Multiplying 
equations (1.13) by ty,» We have, owing to relations (1.5), 


i 4 { 
Nap PAE, Tan Nuapy colt == a. Oopy LR cake a (1.15) 


= — (1, Rovar zi Lg Ryart sa Lape) = 0, 


Ve 7p 


i.e., 1* satisfies equations (1.14) as well. 

Conversely, let there exist a vector /* satisfying equations (1.14); we will show that 
it also satisfies equations (1.13). Performing the calculations of (1.15) in reverse order, 
we obtain relations equivalent to the initial condition: 

Magy ee, = 0, 


Let us multiply these equations by n&*8%. Owing to the identities (1.7) we obtain 
the following result; 


€ Vee 
by LeREYT = 0, 
whence 


L*RE5, = 0 


? 


i.e., [* satisfies condition (1.13), which proves the lemma. 

Note that in the formulation of the lemma it is essential that 2* be contracted with 
one of the first two indices of R,,/,, in view of the inequivalence of the pairs af 
and y6, which follows from definition (1.1). In the case of Einstein spaces, on the 
other hand, owing to Theorem 1 the index pairs are equivalent and the conditions 
ig "Roa = 0 and I°*R_,,,=0 become equivalent. 

Theorem 2. In the Einstein spaces *7,; (3.7) the equations 


Rays = 0 (1.16) 
are equivalent to the equations 
LRasys + laRanys +p Paars = 9: (1.17) 


Proof. Let there exist a vector * in *7, satisfying equations (1.16); we will 
show that it satisfies equations (1.17) as well. Multiplying the initial equations 


1 ; 
PR =0 by F "Nweo? WE obtain 


eRe =. (1.18) 
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Owing to Theorem 1 the tensor *R,,,/, is symmetric with respect to the index 
pairs pv and yé, and it follows from (1.19) that 


PR. =O. (1.19) 


But, according to the lemma, condition (1.19) is equivalent to condition (1.17), which 
proves the first half of Theorem 2. 

Conversely, let a certain vector /* in *7, satisfy equations (1.17). According to 
the lemmia, it satisfies equations (1.13); but, owing to Theorem 1, these equations 
are equivalent to the following: 


Thee? pax m= 0, (1.20) 
° > © S.. . P 2 ‘ 
Multiplying (1.20) by n¥°": and using the identity (1.6), we obtain 
beck 


e3a3 1% = 2Rerapl” =— 0, 


i.e., the vector J% also satisfies equations (1.16). Thus Theorem 2 is proved. 
Theorem 3. The following identity holds in the Einstein spaces (3.7): 


sage | Seer a R35 yso- a 2 (Ryu tay = Ree gee + XR avs) =0. (1-21) 


Proof; Differentiating covariantly the Bianchi identies 


R a py3:0 a Ra aseyw = Ra spxs8 =0 (1-22) 
and multiplying the result by the tensor g°°, we obtain 
BR ays, Pa Ropero =: Ra aris0 =0. (1.23) 


Next, applying the Ricci differential identities (see /9,’) to the Riemann tensor 
cote Perey Sack sock 
27 38,[¢0] = Rays Rosa. Ronse sop. + Rane? poy. TR yaya Roos. (1.24) 
aby; o8 ? 


respectively, as well as terms quadratic in the Riemann tensor (and not containing 
derivatives of Ryoys)* Allowing further for the identity 


we express the second and third terms in (1.23) in terms of Rj;,°.,, and R3,° 


Ry py sa _ QR fa:B}s ; (1.25) 
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which follows from the Bianchi identities (1.22), we see that in the Einstein 
spaces (3.7) 


Raby vo =O (1.26) 


Using relations (1.26), we cancel terms of the form arand R yay; gs in the 


equalities obtained and arrive at relations (1.21). Theorem 3 has been proved. 


Appendix II 


We will give the wave equations (12,34) — (12.36) in expanded form for the 
chronometrically invariant components XxX, Yi3* and Z#Jk! of the Riemann tensor. 
To this end let us use the Zel'manov generally covariant definition of gravitational 
waves: 


afy8ie 0. (II.1) 


PR 
The system of twenty equations (II.1) is obviously equivalent to the following 
three systems of equations, the first consisting of six equations, the second of eight 
and the third of six again: 


BR. 06 a a + Hee = 0, (II.2) 

a ee = oh, af RE =0, (11.3) 
sepkijt _. pkisl;o« kaglme __ 

g igs = RN 0 + RN) ss =U. (11.4) 


We will write each of the three systems in chronometrically invariant form. 
To do this we will use definitions (12.13) as well as the formulas for chronometrically 
invariant and spatially-covariant differentiation (see Ch.11). 
Describing equations (11.2) term by term and expressing each term in terms of the 
corresponding chronometric invariants, we obtain the following equations: 
i (i 
Cyn — tata X47 + 2¢¢y, — FY [Ya (Dh + AM] — 
Nk ij ij ; 3 nj ; j ni 
—[F™y XY 4 peax 4 (DE + Atytax™ +4 (Di 4 Adytax™) 4 
+ (9 + D) (2F"¥ 0) — x" (i 4 AZ) — x” (E+ Ady] + 
+ 2(D™™ 4 AP™ ty VD 4 oF, aye) — ym (De + AM) + 
POP MD te Aa Ore bE, ARO Ea 
opr ae (OL a As yaox (oe de Aye 
Ae i Speake kl kL l 
—2X"™ (DE + Ast) (D2, + Af) + 2X7 DD" + Ay A FP) + 
ivi mre l l ‘ 
+ oF Gx 1 97 G) (Rh R™ (pe 4 AM (De 4 A] = 0. (11.5) 
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Here, F!, Di}, Aqj are the vector of gravitational-inertial force, the deformation 
velocity tensor and the tensor of angular velocity of rotation of the three -dimensional 
space of the given reference system = with respect to a locally co-moving geodesic 
system %, /205/ (these quantities were introduced in Chapter 12). Despite the 
unwieldy aspect of these equations, each of the terms that enter into them has a fairly 
clear physical interpretation. The first term evidently represents the result of the 
action of the wave operator (12.21) on a chronometrically invariant wave function. 
Let us look at the physical meaning of other chronometrically invariant operators in 
equations (II.5). Thus the relativistic operator 


*V; a F, ’ 


contracted with a certain chronometrically invariant three-dimensional vector ti, 
expresses the "physical divergence” of the vector ¢, the difference between it and 
the mathematical divergence *V;t' being due to the fact that at different points of 
the boundary of an element of three -dimensional volume, the value dt of the 
chronometrically invariant time interval 


varies for the same value of di. The relativistic operator of "physical differentiation" 
with respect to the time *@ + D differs from the mathematical (chronometrically 
invariant) operator of differentiation with respect of the time *9 in that it allows for 
deformation with time of the spatial coordinate network in which the function to be 
differentiated is specified. 

Analogously, expressing equations (II.3) in chronometrically invariant form, we 
obtain the wave equations for the chronometrically invariant tensor yk ; 


Or ety, — *080) ¥* + a + D) (Y™ (DI + Ast) — 
YP De 4 Ay ¥™E (DE + Ast) + QRtixaik — ptgrtiky) — 
— (y, =F) (2 (Dt AM) — x8 (DT + A) + 
+ Zyit (D™ + A) + (Age, — DB) tty, — Fr) oe + 
+ (D+ AM yey KM (DME aMyey xy petty ye 
+ Dey + (Di + Ais ay" (Di 4 Adyoy™® 4 
+ (DE + Asha” 4 aRlitgxit 4 pMag-iikyy 4 
ie he omy sek Loph yeh cat Aa) BI ao pa 
+ 2(ay""lird 4 20 pi 4 act y — ZF ot 4 ai yp 
+ 2" (DE + ARI} +2 2PEX DEL 42 4 


4X" cr (Di + sty — Ft + Aziyy + 
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(1.6) 


Lastly, equations (11.4) yield the following system of chronometrically invariant 


wave equations for the tensor Zkil; 


on™” Vn — *9*9) grits. (*0 + D) p2yilk el) su ay kilt wily 
. nik . i . , . 
— (Di 4 Asty ZA 
—(*y, ome F) (pb a yi + (p™ 4 A") yr - 
—(D yn L Amy yuk _ (p™ Pe A”) yeti cs 
é nl rik nk il 
+ (D" 4 A") y YS —(D" ee y¥y ¥ a 
ie n 
- (p™ + A™) yy" me we a at ae ae -_ 
a! [ety gt. paz" > Feagy 2! 7 
n 
— Flegy tk + Fisay* 4 __ plagy*ti 4 
+ (Dk + A;k) egz? (DE + At!) aggzriik ef (Di 4 Asi) aggzinil = 
— (Di + assy saghh'Y + 
Te 
+4 (riplky li 4 Fi pl ykliy te Fo (yue ip 4 A™) _ 
ep gud iu wae 
(o™ 4 A” _. yep" 4 A” ee 
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+2400" + A) (x (Dh 4 ant) — x4 (DEE ashy] + 
kj . ' 
+ (D™ + A") (x9 (DEF As — x DL 4 Ady) + 
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+Y 


+ (Di + As) aynillek) — plykis) 
of (pk + A‘*) 2yMlipa abi F lyi m i 
+ (Di, + A; pty (2ynkl pil Ff yin) Sh 
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4 2°™ (Di 4 AS )| — (DE 4 Ady (2h (DI, + Ast) + 
ik , i +i jie z 
+ ZP™ (Di, + Amd + (Dm + Am.) (20 Wi Ag) — 
ioe 
— 2P (DA A AR) == 0. (11.7) 


The physical interpretation of equations (II.6) and (I1.7) is analogous to the above. 
Equations (11.5) — (11.7) constitute the complete system of chronometrically invariant 
equations describing gravitational-inertial waves in the specific frame of reference. 
Since these equations are none other than the generally covariant equations (II.1) 
written down in an arbitrary (fixed) frame of reference, it follows that they are 
satisfied for any empty space-time of Petrov type W (in the given frame of reference). 
Thus equations (1J.5) — (11.7) can serve as chronometrically invariant characteristic 
of type WM gravitational fields in a vacuum. 
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